Assessing the influence of sample tagging and library preparation on DNA metabarcoding

Metabarcoding is increasingly used to assess species diversity by high‐throughput sequencing where millions of sequences can be generated in parallel and multiple samples can be analysed in one sequencing run. Generating amplified fragments with a unique sequence identifier ('tag') for each sample is crucial, as it allows assigning sequences to the original samples. The tagging through so‐called fusion primers is a fast and cheap alternative to commercially produced ligation‐based kits. However, little is known about potential bias and inconsistencies introduced by the long nucleotide tail attached to those primers, which could lead to deficient reports of community composition in metabarcoding studies. We therefore tested the consistency and taxa detection efficiency of fusion primers in (1) a one‐step and (2) two‐step PCR protocol as well as (3) a commercially manufactured Illumina kit using mock communities of known composition. The Illumina kit delivered the most consistent results and detected the highest number of taxa. However, success of the two‐step PCR approach was only marginally lower compared to the kit with the additional advantage of a much more competitive price per library. While most taxa were also detected with the one‐step PCR approach, the consistency between replicates including read abundance was substantially lower. Our results highlight that method choice depends on the precision needed for analysis as well as on economic considerations and recommend the Illumina kit to obtain most accurate results and the two‐step PCR approach as a much cheaper yet very robust alternative.     

Highly-multiplexed barcode sequencing: an efficient method for parallel analysis of pooled samples

Next-generation sequencing has proven an extremely effective technology for molecular counting applications where the number of sequence reads provides a digital readout for RNA-seq, ChIP-seq, Tn-seq and other applications. The extremely large number of sequence reads that can be obtained per run permits the analysis of increasingly complex samples. For lower complexity samples, however, a point of diminishing returns is reached when the number of counts per sequence results in oversampling with no increase in data quality. A solution to making next-generation sequencing as efficient and affordable as possible involves assaying multiple samples in a single run. Here, we report the successful 96-plexing of complex pools of DNA barcoded yeast mutants and show that such ‘Bar-seq’ assessment of these samples is comparable with data provided by barcode microarrays, the current benchmark for this application. The cost reduction and increased throughput permitted by highly multiplexed sequencing will greatly expand the scope of chemogenomics assays and, equally importantly, the approach is suitable for other sequence counting applications that could benefit from massive parallelization.     

HPLC purification of polymerase chain reaction products for direct sequencing

Same day PCR amplification and sequencing is desired in situations where one needs to sequence a number of PCR products. The rapid, high-yield purification of PCR products via the use of high performance, anion-exchange chromatography yields sequencing results comparable to those obtained from techniques requiring subcloning of the PCR product. This can be achieved by standard dideoxynucleotide sequencing technology without the need to prepare prelabeled primers and additional internal primers or to gel purify the PCR product. In addition, this chromatographic technique offers the potential of isolating several PCR products from the same amplification mixture.     

Comparison of Normalization Methods for Construction of Large,Multiplex Amplicon Pools for Next-Generation Sequencing

Constructing mixtures of tagged or bar-coded DNAs for sequencing is an important requirement for the efficient use of next-generation sequencers in applications where limited sequence data are required per sample. There are many applications in which next-generation sequencing can be used effectively to sequence large mixed samples; an example is the characterization of microbial communities where ≤1,000 sequences per samples are adequate to address research questions. Thus, it is possible to examine hundreds to thousands of samples per run on massively parallel next-generation sequencers. However, the cost savings for efficient utilization of sequence capacity is realized only if the production and management costs associated with construction of multiplex pools are also scalable. One critical step in multiplex pool construction is the normalization process, whereby equimolar amounts of each amplicon are mixed. Here we compare three approaches (spectroscopy, size-restricted spectroscopy, and quantitative binding) for normalization of large, multiplex amplicon pools for performance and efficiency. We found that the quantitative binding approach was superior and represents an efficient scalable process for construction of very large, multiplex pools with hundreds and perhaps thousands of individual amplicons included. We demonstrate the increased sequence diversity identified with higher throughput. Massively parallel sequencing can dramatically accelerate microbial ecology studies by allowing appropriate replication of sequence acquisition to account for temporal and spatial variations. Further, population studies to examine genetic variation, which require even lower levels of sequencing, should be possible where thousands of individual bar-coded amplicons are examined in parallel.Emergent technologies that generate DNA sequence data are designed primarily to perform resequencing projects at reasonable cost. The result is a substantial decrease in per base costs from traditional methods. However, these next-generation platforms do not readily accommodate projects that require obtaining moderate amounts of sequence from large numbers of samples. These platforms also have per run costs that are significant and generally preclude large numbers of single-sample, nonmultiplexed runs. One example of research that is not readily supported is rRNA-directed metagenomics study of some human clinical samples or environmental rRNA analysis of samples from communities with low community diversity that require only thousands of sequences. Thus, strategies to utilize next-generation DNA sequencers efficiently for applications that require lower throughput are critical to capitalize on the efficiency and cost benefits of next-generation sequencing platforms.Directed metagenomics based on amplification of rRNA genes is an important tool to characterize microbial communities in various environmental and clinical settings. In diverse environmental samples, large numbers of sequences are required to fully characterize the microbial communities (15). However, a lower number of sequences is generally adequate to answer specific research questions. In addition, the levels of diversity in human clinical samples are usually lower than what is observed in environmental samples (for example, see reference 7).The Roche 454 genome sequencer system FLX pyrosequencer (which we will refer to as 454 FLX hereafter) is the most useful platform for rRNA-directed metagenomics because it currently provides the longest read lengths of any next-generation sequencing platform (1, 14). Computational analysis has shown that the 250-nucleotide read length (available from the 454 FLX-LR chemistry) is adequate for identification of bacteria if the amplified region is properly positioned within variable regions of the small-subunit rRNA (SSU-rRNA) gene (9, 10).In this study, we used the 454 FLX-LR genome sequencing platform and chemistry, which provides >400,000 sequences of ∼250 bp per run. After we conducted this study, a new reagent set (454 FLX-XLR titanium chemistry) was released, which further increases reads to >1,000,000 and read lengths to >400 bp (Roche). The 454 FLX platform dramatically reduces per base costs of obtaining sequence, and physical separation into between 2 and 16 lanes is available; this physical separation on the plate reduces sequencing output overall, up to 40% comparing 2 lanes versus 16 lanes. For applications where modest sequencing depth (∼1,000 sequences per sample) is adequate to address research questions, physical separation does not allow adequate sample multiplexing because even a 1/16 454 FLX-LR plate run is expected to produce ∼15,000 reads. Further, the utility of the platform as a screening tool at 16-plex is limited by cost per run.A solution to make next-generation sequencing economical for projects such as rRNA-directed metagenomics is to use bar-coded primers to multiplex amplicon pools so they can be sequenced together and computationally separated afterward (6). To successfully accomplish this strategy, precise normalization of the DNA concentrations of the individual amplicons in the multiplex pools is essential for effective multiplex sequencing when large numbers of pooled samples are sequenced in parallel. There are several potential methods available for normalizing concentrations of amplicons included in multiplex pools, but the relative and absolute performance of each approach has not been compared.In this study, we present a direct quantitative comparison of three available methods for amplicon pool normalization for downstream next-generation sequencing. The central goal of the study was to identify the most effective method for normalizing multiplex pools containing >100 individual amplicons. We evaluated each pooling approach by 454 sequencing and compared the observed frequencies of sequences from different pooled bar-coded amplicons. From these data, we determined the efficacy of each method based on the following factors: (i) how well normalized the sequences within the pool were, (ii) the proportion of samples failing to meet a minimum threshold of sequences per sample, and (iii) the overall efficiency (speed and labor required) of the process to multiplex samples.     

在单个PCR管内快捷完成定点突变

采用大引物方法,利用质粒多克隆位点两侧的普通测序引物作为旁侧引物,在单个PCR管内,经2个步骤共34个循环进行定点突变. 该方法通过优化模板和引物的量达到降低PCR循环次数, 通过加入10个在68℃复性条件下的PCR循环达到增加突变效率而无需胶纯化.本方法达到平均62%的突变效率,而且全长扩增产物的产率很高.     

Microfluidic PCR-based target enrichment: A case study in two rapid radiations of Commiphora (Burseraceae) from Madagascar

Developing effective and cost-efficient multilocus nuclear datasets for angiosperm species is a continuing challenge to the systematics community. Here we describe the development and validation of a novel set of 91 nuclear markers for PCR-based target enrichment. Using microfluidic PCR and Illumina MiSeq, we generated nuclear, subgenomic libraries for 96 species simultaneously and sequenced them for a total cost of ca. $6000 USD. Approximately half of these costs include reusable reagents (primers, barcodes, and custom sequencing primers) and taxon sampling could be increased by an order of magnitude to maximize sequencing depth efficiency. The principle benefit of microfluidic PCR over alternative target enrichment strategies is that it bypasses costly library preparation. After sequencing, we evaluated the ability of the loci to resolve species level relationships within two recently radiated lineages of endemic Madagascan Commiphora Jacq. (Burseraceae) species. Our results demonstrate that (i) effective nuclear markers can be designed for non-model angiosperm taxa from these publicly available datasets; (ii) that microfluidic PCR amplification followed by high throughput sequencing can produce highly complete taxon by locus sequence data matrices with minimal resource investment; and (iii) that these numerous nuclear phylogenomic markers can improve our understanding of phylogenetic relationships withinCommiphora. We provide a synopsis of ongoing activities to enhance this microfluidic PCR-based target enrichment strategy through broader primer assays, multiplexing, and increased efficiency of sequencing depth.     

Large scale sequencing projects using rapidly prepared double-stranded plasmid DNA.

We have developed a simple rapid plasmid DNA mini-preparation method which yields DNA of sufficient quality to be used in large scale sequencing projects. The method, which is a modification of the alkaline method of Birnboim and Doly (1979), requires less than two hours. We have eliminated the use of organic extractions, RNase digestion and alkaline denaturation of the DNA for annealing of the primer. The proportion of supercoiled plasmid DNA obtained is close to 100%. Greater than 80% of the clones yield at least 500 bp of sequence information per primer. The sequencing reactions from these double-stranded templates can be done on both strands using the universal and reverse sequence primers with the usual two reactions per primer, one to read close to the primer and one to read far from it. Thus, each clone yields at least 1 kb of sequence information. The preparation of the templates and the sequencing reactions can be done in less than three hours so that the sequencing gel can be run the same day.     

Phylum- and class-specific PCR primers for general microbial community analysis

Amplification of a particular DNA fragment from a mixture of organisms by PCR is a common first step in methods of examining microbial community structure. The use of group-specific primers in community DNA profiling applications can provide enhanced sensitivity and phylogenetic detail compared to domain-specific primers. Other uses for group-specific primers include quantitative PCR and library screening. The purpose of the present study was to develop several primer sets targeting commonly occurring and important groups. Primers specific for the 16S ribosomal sequences of Alphaproteobacteria, Betaproteobacteria, Bacilli, Actinobacteria, and Planctomycetes and for parts of both the 18S ribosomal sequence and the internal transcribed spacer region of Basidiomycota were examined. Primers were tested by comparison to sequences in the ARB 2003 database, and chosen primers were further tested by cloning and sequencing from soil community DNA. Eighty-five to 100% of the sequences obtained from clone libraries were found to be placed with the groups intended as targets, demonstrating the specificity of the primers under field conditions. It will be important to reevaluate primers over time because of the continual growth of sequence databases and revision of microbial taxonomy.     

Target region amplification polymorphism: A novel marker technique for plant genotyping

The advent of large-scale DNA sequencing technology has generated a tremendous amount of sequence information for many important organisms. We have developed a rapid and efficient PCR-based technique, which uses bioinformatics tools and expressed sequence tag (EST) database information to generate polymorphic markers around targeted candidate gene sequences. This target region amplification polymorphism (TRAP) technique uses 2 primers of 18 nucleotides to generate markers. One of the primers, the fixed primer, is designed from the targeted EST sequence in the database; the second primer, the arbitrary primer, is an arbitrary sequence with either an AT-or GC-rich core to anneal with an intron or exon, respectively. PCR amplification is run for the first 5 cycles with an annealing temperature of 35°C, followed by 35 cycles with an annealing temperature of 50°C. For different plant species, each PCR reaction can generate as many as 50 scorable fragments with sizes ranging from 50–900 bp when separated on a 6.5% polyacrylamide sequencing gel. The TRAP technique should be useful in genotyping germplasm collections and in tagging genes governing desirable agronomic traits of crop plants.     

A Next Generation Semiconductor Based Sequencing Approach for the Identification of Meat Species in DNA Mixtures

The identification of the species of origin of meat and meat products is an important issue to prevent and detect frauds that might have economic, ethical and health implications. In this paper we evaluated the potential of the next generation semiconductor based sequencing technology (Ion Torrent Personal Genome Machine) for the identification of DNA from meat species (pig, horse, cattle, sheep, rabbit, chicken, turkey, pheasant, duck, goose and pigeon) as well as from human and rat in DNA mixtures through the sequencing of PCR products obtained from different couples of universal primers that amplify 12S and 16S rRNA mitochondrial DNA genes. Six libraries were produced including PCR products obtained separately from 13 species or from DNA mixtures containing DNA from all species or only avian or only mammalian species at equimolar concentration or at 1:10 or 1:50 ratios for pig and horse DNA. Sequencing obtained a total of 33,294,511 called nucleotides of which 29,109,688 with Q20 (87.43%) in a total of 215,944 reads. Different alignment algorithms were used to assign the species based on sequence data. Error rate calculated after confirmation of the obtained sequences by Sanger sequencing ranged from 0.0003 to 0.02 for the different species. Correlation about the number of reads per species between different libraries was high for mammalian species (0.97) and lower for avian species (0.70). PCR competition limited the efficiency of amplification and sequencing for avian species for some primer pairs. Detection of low level of pig and horse DNA was possible with reads obtained from different primer pairs. The sequencing of the products obtained from different universal PCR primers could be a useful strategy to overcome potential problems of amplification. Based on these results, the Ion Torrent technology can be applied for the identification of meat species in DNA mixtures.     

The real-time polymerase chain reaction

The scientific, medical, and diagnostic communities have been presented the most powerful tool for quantitative nucleic acids analysis: real-time PCR [Bustin, S.A., 2004. A-Z of Quantitative PCR. IUL Press, San Diego, CA]. This new technique is a refinement of the original Polymerase Chain Reaction (PCR) developed by Kary Mullis and coworkers in the mid 80:ies [Saiki, R.K., et al., 1985. Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia, Science 230, 1350], for which Kary Mullis was awarded the 1993 year's Nobel prize in Chemistry. By PCR essentially any nucleic acid sequence present in a complex sample can be amplified in a cyclic process to generate a large number of identical copies that can readily be analyzed. This made it possible, for example, to manipulate DNA for cloning purposes, genetic engineering, and sequencing. But as an analytical technique the original PCR method had some serious limitations. By first amplifying the DNA sequence and then analyzing the product, quantification was exceedingly difficult since the PCR gave rise to essentially the same amount of product independently of the initial amount of DNA template molecules that were present. This limitation was resolved in 1992 by the development of real-time PCR by Higuchi et al. [Higuchi, R., Dollinger, G., Walsh, P.S., Griffith, R., 1992. Simultaneous amplification and detection of specific DNA-sequences. Bio-Technology 10(4), 413-417]. In real-time PCR the amount of product formed is monitored during the course of the reaction by monitoring the fluorescence of dyes or probes introduced into the reaction that is proportional to the amount of product formed, and the number of amplification cycles required to obtain a particular amount of DNA molecules is registered. Assuming a certain amplification efficiency, which typically is close to a doubling of the number of molecules per amplification cycle, it is possible to calculate the number of DNA molecules of the amplified sequence that were initially present in the sample. With the highly efficient detection chemistries, sensitive instrumentation, and optimized assays that are available today the number of DNA molecules of a particular sequence in a complex sample can be determined with unprecedented accuracy and sensitivity sufficient to detect a single molecule. Typical uses of real-time PCR include pathogen detection, gene expression analysis, single nucleotide polymorphism (SNP) analysis, analysis of chromosome aberrations, and most recently also protein detection by real-time immuno PCR.     

Microbiome profiling by illumina sequencing of combinatorial sequence-tagged PCR products

We developed a low-cost, high-throughput microbiome profiling method that uses combinatorial sequence tags attached to PCR primers that amplify the rRNA V6 region. Amplified PCR products are sequenced using an Illumina paired-end protocol to generate millions of overlapping reads. Combinatorial sequence tagging can be used to examine hundreds of samples with far fewer primers than is required when sequence tags are incorporated at only a single end. The number of reads generated permitted saturating or near-saturating analysis of samples of the vaginal microbiome. The large number of reads allowed an in-depth analysis of errors, and we found that PCR-induced errors composed the vast majority of non-organism derived species variants, an observation that has significant implications for sequence clustering of similar high-throughput data. We show that the short reads are sufficient to assign organisms to the genus or species level in most cases. We suggest that this method will be useful for the deep sequencing of any short nucleotide region that is taxonomically informative; these include the V3, V5 regions of the bacterial 16S rRNA genes and the eukaryotic V9 region that is gaining popularity for sampling protist diversity.     

Optimization strategies for the polymerase chain reaction

The GeneAmp polymerase chain reaction (PCR) process has now become a key procedure in molecular biology research laboratories. The PCR technique is an in vitro method in which genomic or cloned target sequences are specifically enzymatically amplified as directed by a pair of oligonucleotide primers. This technique has been quite robust in the hands of the majority of researchers and is extremely flexible, as evidenced by the increasing number of related PCR formats (i.e., inverse PCR, anchored PCR, asymmetric PCR, labeled primer PCR and RNA-PCR). Today's applications include direct sequencing, genomic cloning, DNA typing, detection of infectious microorganisms, site-directed mutagenesis, prenatal genetic disease research, and analysis of allelic sequence variations. Scientists at Cetus and Perkin-Elmer have collaborated for several years to better understand the interacting biochemical and biophysical parameters which affect PCR optimization. Following are many of the current recommendations, offered with the caveat that our understanding of the PCR process is continually evolving.     

Random PCR-based genome sequencing: a non-divide-and-conquer strategy.

We propose a genome sequencing strategy, which is neither divide-and-conquer (clone by clone) nor the shotgun approach. Random PCR-based and PCR relay sequencing constitute the basis of this novel strategy. Most of the genome is sequenced by the former process that requires only a set of non-specific primers and a template DNA. Random PCR-based sequencing reduces redundancy in sequencing by exploiting known sequence information. The number of primers required for random PCR was significantly diminished by using a combination of primers. The former process can be partially replaced by the shotgun method, if necessary. The gap-filling process can be effectively performed by way of PCR relay. The feasibility of this strategy was demonstrated using the Escherichia coli genome. This strategy enhances the global effort towards genome sequencing by being available through the Internet and by allowing the use of preexisting sequence data.     

Direct PCR sequencing of dystrophin polymorphic CACA alleles after purification to remove shadow bands.

A method is described that allows the sequencing of polymerase chain reaction (PCR) products containing CACA repeats. The method was tested using a DNA polymorphism that exists at the 3' end of the dystrophin gene. This polymorphism consists of a variation in the length of a CACA dinucleotide repeat. Four alleles from a total of 16 individuals were sequenced at this locus after the DNA sequence had been amplified by the PCR. Five examples of each of the common alleles were sequenced. For each allele all five sequences were the same. The only example of a rare allele was also sequenced. The PCR products of DNA sequences containing dinucleotide repeats consist of a number of bands differing by 2 bp below the most intense main band. Previously, direct sequencing of the PCR products lead to ambiguities and smearing at and above the CACA repeat. In this paper, the main PCR band was cut out of a sequencing gel and directly sequenced to give a clear DNA sequence. Our results indicate that for a particular allele, all individuals had exactly the same DNA sequence. This implies that with the appropriate choice of oligonucleotide primers, polymorphisms could be detected without electrophoresis.     

Phylum- and Class-Specific PCR Primers for General Microbial Community Analysis

Amplification of a particular DNA fragment from a mixture of organisms by PCR is a common first step in methods of examining microbial community structure. The use of group-specific primers in community DNA profiling applications can provide enhanced sensitivity and phylogenetic detail compared to domain-specific primers. Other uses for group-specific primers include quantitative PCR and library screening. The purpose of the present study was to develop several primer sets targeting commonly occurring and important groups. Primers specific for the 16S ribosomal sequences of Alphaproteobacteria, Betaproteobacteria, Bacilli, Actinobacteria, and Planctomycetes and for parts of both the 18S ribosomal sequence and the internal transcribed spacer region of Basidiomycota were examined. Primers were tested by comparison to sequences in the ARB 2003 database, and chosen primers were further tested by cloning and sequencing from soil community DNA. Eighty-five to 100% of the sequences obtained from clone libraries were found to be placed with the groups intended as targets, demonstrating the specificity of the primers under field conditions. It will be important to reevaluate primers over time because of the continual growth of sequence databases and revision of microbial taxonomy.     

Multiplex restriction amplicon sequencing: a novel next‐generation sequencing‐based marker platform for high‐throughput genotyping

To enable rapid selection of traits in marker‐assisted breeding, markers must be technically simple, low‐cost, high‐throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3’‐ends, preceded by 6‐10 bases of specific or degenerate nucleotide sequences and then by a unique M13‐tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next‐generation sequencing‐based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.     

Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads

Microbial communities host unparalleled taxonomic diversity. Adequate characterization of environmental and host-associated samples remains a challenge for microbiologists, despite the advent of 16S rRNA gene sequencing. In order to increase the depth of sampling for diverse bacterial communities, we developed a method for sequencing and assembling millions of paired-end reads from the 16S rRNA gene (spanning the V3 region; ～200 nucleotides) by using an Illumina genome analyzer. To confirm reproducibility and to identify a suitable computational pipeline for data analysis, sequence libraries were prepared in duplicate for both a defined mixture of DNAs from known cultured bacterial isolates (>1 million postassembly sequences) and an Arctic tundra soil sample (>6 million postassembly sequences). The Illumina 16S rRNA gene libraries represent a substantial increase in number of sequences over all extant next-generation sequencing approaches (e.g., 454 pyrosequencing), while the assembly of paired-end 125-base reads offers a methodological advantage by incorporating an initial quality control step for each 16S rRNA gene sequence. This method incorporates indexed primers to enable the characterization of multiple microbial communities in a single flow cell lane, may be modified readily to target other variable regions or genes, and demonstrates unprecedented and economical access to DNAs from organisms that exist at low relative abundances.     

Genome-Tagged Amplification (GTA): a PCR-based method to prepare sample-tagged amplicons from hundreds of individuals for next generation sequencing

Sampling the sequence of a relatively small fraction of the genome in large numbers of individuals is an important objective for population genetics and association genetics approaches. However, currently available ‘sequence capture’ methods either require expensive instrumentation or have problems dealing with high sample numbers and relatively small target sizes. We have developed Genome-Tagged Amplification (GTA) as a flexible PCR-based method for preparing pools of hundreds of amplicons from hundreds of samples for next generation sequencing. The method involves tagging of genomic DNA with barcode adapters at restriction sites, followed by PCR amplification from flanking DNA. It is freely scalable for both sample number and amplicon number and has no specialized equipment requirement. An optimized protocol is presented which provides a matrix of 96 × 192 combinations of samples x amplicons, corresponding to a complete 454 Titanium run. Initially, we used 454 sequencing; however, GTA could easily be adapted to Illumina sequencing platforms as read lengths have significantly increased in this system.