Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within‐population ITS1 sequence variation in a microparasite infecting Daphnia

• Published in: Molecular ecology resources Vol. 15; no. 6; pp. 1385 - 1395
• Main Authors: González-Tortuero, E., Rusek, J., Petrusek, A., Gießler, S., Lyras, D., Grath, S., Castro-Monzón, F., Wolinska, J.
• Format: Journal Article
• Language: English
• Published: England Blackwell Pub 01-11-2015; Blackwell Publishing Ltd; Wiley Subscription Services, Inc
• Subjects: amplicon sequencing; Animals; Caullerya mesnili; Cladocera; Computational Biology - methods; Daphnia; Daphnia - parasitology; DNA, Ribosomal Spacer - chemistry; DNA, Ribosomal Spacer - genetics; Genetic Variation; High-Throughput Nucleotide Sequencing; Ichthyosporea; ITS1 region; Mesomycetozoea - classification; Mesomycetozoea - genetics; Next Generation Sequencing; pipeline; Sequence Analysis, DNA; Software; pipeline; ITS1 region; amplicon sequencing; Caullerya mesnili; Ichthyosporea; Next Generation Sequencing
• ISSN: 1755-098X; 1755-0998
• DOI: 10.1111/1755-0998.12396

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within‐population variation. Additionally, a public Illumina data set was used to validate the pipeline on community‐level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within‐population structure but also the successful application of the QRS pipeline on Illumina‐generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences.