Length and GC-biases during sequencing library amplification: a comparison of various polymerase-buffer systems with ancient and modern DNA sequencing libraries

High-throughput sequencing technologies frequently necessitate the use of PCR for sequencing library amplification. PCR is a sometimes enigmatic process and is known to introduce biases. Here we perform a simple amplification-sequencing assay using 10 commercially available polymerase-buffer systems to amplify libraries prepared from both modern and ancient DNA. We compare the performance of the polymerases with respect to a previously uncharacterized template length bias, as well as GC-content bias, and find that simply avoiding certain polymerase can dramatically decrease the occurrence of both. For amplification of ancient DNA, we found that some commonly used polymerases strongly bias against amplification of endogenous DNA in favor of GC-rich microbial contamination, in our case reducing the fraction of endogenous sequences to almost half.     

Sources of PCR-induced distortions in high-throughput sequencing data sets

PCR permits the exponential and sequence-specific amplification of DNA, even from minute starting quantities. PCR is a fundamental step in preparing DNA samples for high-throughput sequencing. However, there are errors associated with PCR-mediated amplification. Here we examine the effects of four important sources of error—bias, stochasticity, template switches and polymerase errors—on sequence representation in low-input next-generation sequencing libraries. We designed a pool of diverse PCR amplicons with a defined structure, and then used Illumina sequencing to search for signatures of each process. We further developed quantitative models for each process, and compared predictions of these models to our experimental data. We find that PCR stochasticity is the major force skewing sequence representation after amplification of a pool of unique DNA amplicons. Polymerase errors become very common in later cycles of PCR but have little impact on the overall sequence distribution as they are confined to small copy numbers. PCR template switches are rare and confined to low copy numbers. Our results provide a theoretical basis for removing distortions from high-throughput sequencing data. In addition, our findings on PCR stochasticity will have particular relevance to quantification of results from single cell sequencing, in which sequences are represented by only one or a few molecules.     

Amplification of Nucleic Acids by Polymerase Chain Reaction (PCR) and Other Methods and their Applications

Abstract

The in vitro replication of DNA, principally using the polymerase chain reaction (PCR), permits the amplification of defined sequences of DNA. By exponentially amplifying a target sequence, PCR significantly enhances the probability of detecting target gene sequences in complex mixtures of DNA. It also facilitates the cloning and sequencing of genes. Amplification of DNA by PCR and other newly developed methods has been applied in many areas of biological research, including molecular biology, biotechnology, and medicine, permitting studies that were not possible before. Nucleic acid amplification has added a new and revolutionary dimension to molecular biology. This review examines PCR and other in vitro nucleic acid amplification methodologies—examining the critical parameters and variations and their widespread applications—giving the strengths and limitations of these methodologies.     

Re-usable DNA template for the polymerase chain reaction.

DNA covalently bound to an uncharged nylon membrane was used for consecutive amplifications of several different genes by PCR. Successful PCR amplifications were obtained for membrane-bound genomic and plasmid DNA. Membrane-bound genomic DNA templates were re-used at least 15 times for PCR with specific amplification of the desired gene each time. PCR amplifications of specific sequences of p53, p16, CYP1A1, CYP2D6, GSTM1 and GSTM3 were performed independently on the same strips of uncharged nylon membrane containing genomic DNA. PCR products were subjected to restriction fragment length polymorphism analysis, single-strand conformational polymorphism analysis and/or dideoxy sequencing to confirm PCR-amplified gene sequences. We found that PCR fragments obtained by amplification from bound genomic DNA as template were identical in sequence to those of PCR products obtained from free genomic DNA in solution. PCR was performed using as little as 5 ng genomic or 4 fg plasmid DNA bound to membrane. These results suggest that DNA covalently bound to membrane can be re-used for sample-specific PCR amplifications, providing a potentially unlimited source of DNA for PCR.     

Development of a comprehensive detection method for medicinal and toxic plant species

Pharmacologically active ingredients in plants can cause significant morbidity through their increasingly common use in herbal alternative medicines and dietary supplements. Monitoring consumer products for the presence of toxic plants is encumbered by the lack of rapid and specific assays. To create a sensitive, reliable, fast, and broad-spectrum assay for medicinal or toxic plant species, we tested multiplexed ligation-dependent probe amplification (MLPA), which requires partial genomic DNA sequences from species of plants that are not well represented in currently available genetic databases. Genomic DNA was obtained from 21 species of medicinal and/or toxic plants. The PCR products were amplified from these plants and cloned for sequencing. The MLPA method was successful with DNA samples from many different species. The use of a microarray to facilitate screening of potentially thousands of plants in a single assay also was successful. The combination of the specificity of the MLPA assay with the broad-scale capabilities of microarray technology should make this an especially useful tool in screening in foods and commercial herbal preparations to identify the plant compounds actually present. Other applications could potentially extend to the identification of any plant species in samples for academic botanical studies and for biodefense and forensics applications.     

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.     

A highly specific microarray method for point mutation detection

Improvements of microarray techniques for genotyping purposes have focused on increasing the reliability of this method. Here we report the development of a genotyping method where a microarray was spotted with stemloop probes, especially designed to optimize the hybridization specificity of complementary DNA sequences. This accurate method was used to screen for four common disease-causing mutations involved in a neurological disorder called Charcot-Marie-Tooth disease (CMT). Healthy individuals' and patients' DNA were amplified and labeled by PCR and hybridized on microarray. The spot signal intensities were 81 to 408 times greater for perfect compared with mismatched target sequences, differing by only one nucleotide (discrimination ratio) for healthy individual "homozygous" DNA. On the other hand, "heterozygous" mutant DNA samples gave rise to signal intensity ratios close to 1, as expected. The genotypes obtained by this method were perfectly consistent with those determined by direct PCR sequencing. Cross-hybridization rates were very low, resulting in further multiplexing improvements. In this study, we also demonstrated the feasibility of real-time hybridization detection of labeled synthetic oligonucleotides with concentrations as low as 2.5 nM.     

High-Density Universal 16S rRNA Microarray Analysis Reveals Broader Diversity than Typical Clone Library When Sampling the Environment

Molecular approaches aimed at detection of a broad-range of prokaryotes in the environment routinely rely on classifying heterogeneous 16S rRNA genes amplified by polymerase chain reaction (PCR) using primers with broad specificity. The general method of sampling and categorizing DNA has been to clone then sequence the PCR products. However, the number of clones required to adequately catalog the majority of taxa in a sample is unwieldy. Alternatively, hybridizing target sequences to a universal 16S rRNA gene microarray may provide a more rapid and comprehensive view of prokaryotic community composition. This study investigated the breadth and accuracy of a microarray in detecting diverse 16S rRNA gene sequence types compared to clone-and-sequencing using three environmental samples: urban aerosol, subsurface soil, and subsurface water. PCR products generated from universal 16S rRNA gene-targeted primers were classified by using either the clone-and-sequence method or by hybridization to a novel high-density microarray of 297,851 probes complementary to 842 prokaryotic subfamilies. The three clone libraries comprised 1391 high-quality sequences. Approximately 8% of the clones could not be placed into a known subfamily and were considered novel. The microarray results confirmed the majority of clone-detected subfamilies and additionally demonstrated greater amplicon diversity extending into phyla not observed by the cloning method. Sequences matching operational taxonomic units within the phyla Nitrospira, Planctomycetes, and TM7, which were uniquely detected by the array, were verified with specific primers and subsequent amplicon sequencing. Subfamily richness detected by the array corresponded well with nonparametric richness predictions extrapolated from clone libraries except in the water community where clone-based richness predictions were greatly exceeded. It was concluded that although the microarray is unreliable in identifying novel prokaryotic taxa, it reveals greater diversity in environmental samples than sequencing a typically sized clone library. Furthermore, the microarray allowed samples to be rapidly evaluated with replication, a significant advantage in studies of microbial ecology.     

Environmental whole-genome amplification to access microbial populations in contaminated sediments

Low-biomass samples from nitrate and heavy metal contaminated soils yield DNA amounts that have limited use for direct, native analysis and screening. Multiple displacement amplification (MDA) using phi29 DNA polymerase was used to amplify whole genomes from environmental, contaminated, subsurface sediments. By first amplifying the genomic DNA (gDNA), biodiversity analysis and gDNA library construction of microbes found in contaminated soils were made possible. The MDA method was validated by analyzing amplified genome coverage from approximately five Escherichia coli cells, resulting in 99.2% genome coverage. The method was further validated by confirming overall representative species coverage and also an amplification bias when amplifying from a mix of eight known bacterial strains. We extracted DNA from samples with extremely low cell densities from a U.S. Department of Energy contaminated site. After amplification, small-subunit rRNA analysis revealed relatively even distribution of species across several major phyla. Clone libraries were constructed from the amplified gDNA, and a small subset of clones was used for shotgun sequencing. BLAST analysis of the library clone sequences showed that 64.9% of the sequences had significant similarities to known proteins, and "clusters of orthologous groups" (COG) analysis revealed that more than half of the sequences from each library contained sequence similarity to known proteins. The libraries can be readily screened for native genes or any target of interest. Whole-genome amplification of metagenomic DNA from very minute microbial sources, while introducing an amplification bias, will allow access to genomic information that was not previously accessible. The reported SSU rRNA sequences and library clone end sequences are listed with their respective GenBank accession numbers, DQ 404590 to DQ 404652, DQ 404654 to DQ 404938, and DX 385314 to DX 389173.     

Barcoded primers used in multiplex amplicon pyrosequencing bias amplification

"Barcode-tagged" PCR primers used for multiplex amplicon sequencing generate a thus-far-overlooked amplification bias that produces variable terminal restriction fragment length polymorphism (T-RFLP) and pyrosequencing data from the same environmental DNA template. We propose a simple two-step PCR approach that increases reproducibility and consistently recovers higher genetic diversity in pyrosequencing libraries.     

Polymerase chain reaction and its applications: special emphasis on its role in embryo sexing

The polymerase chain reaction (PCR) has developed into one of the most promising methods for in vitro enzymatic amplification of DNA and has found widespread application in DNA cloning, sequencing and mutagenesis related studies. This innovative technique can selectively amplify a single target DNA molecule a billion-fold in a span of a few hours. Amplification of specific DNA sequences by PCR is useful in identification of sex, novel genes, pathogens and diseases. PCR has facilitated the establishment of evolutionary relationships among species and in revealing structural intricacies of single cells. In this article we review some of the major advances and applications of PCR, especially, its role in embryo sexing.     

Production of complex nucleic acid libraries using highly parallel in situ oligonucleotide synthesis

Generation of complex libraries of defined nucleic acid sequences can greatly aid the functional analysis of protein and gene function. Previously, such studies relied either on individually synthesized oligonucleotides or on cellular nucleic acids as the starting material. As each method has disadvantages, we have developed a rapid and cost-effective alternative for construction of small-fragment DNA libraries of defined sequences. This approach uses in situ microarray DNA synthesis for generation of complex oligonucleotide populations. These populations can be recovered and either used directly or immortalized by cloning. From a single microarray, a library containing thousands of unique sequences can be generated. As an example of the potential applications of this technology, we have tested the approach for the production of plasmids encoding short hairpin RNAs (shRNAs) targeting numerous human and mouse genes. We achieved high-fidelity clone retrieval with a uniform representation of intended library sequences.     

ROLL: a method of preparation of gene-specific oligonucleotide libraries

The selection of nucleic acid sequences capable of specifically and efficiently hybridizing to target sequences is crucial to the success of many applications, including microarrays, PCR and other amplification procedures, antisense inhibition, ribozyme-mediated cleavage, and RNA interference (RNAi). Methods of selection using nucleotide sequence libraries have several advantages over rational approaches using defined sequences. However, the high complexity of completely random (degenerate) libraries and their high toxicity in cell-based assays make their use in many applications impractical. Gene-specific oligonucleotide libraries, which contain all possible sequences of a certain length occurring within a given gene, have much lower complexity and, thus, can significantly simplify and accelerate sequence screening. Here, we describe a new method for the preparation of gene-specific libraries using the ligation of randomized oligonucleotide probes hybridized adjacently on target polynucleotide templates followed by PCR amplification. We call this method random oligonucleotide ligated libraries (ROLL).     

Detection of bacterial pathogens in municipal wastewater using an oligonucleotide microarray and real-time quantitative PCR

As a first step toward building a comprehensive microarray, two low density DNA microarrays were constructed and evaluated for the accurate detection of wastewater pathogens. The first one involved the direct hybridization of wastewater microbial genomic DNA to the functional gene probes while the second involved PCR amplification of 23S ribosomal DNA. The genomic DNA microarray employed 10 functional genes as detection targets. Sensitivity of the microarray was determined to be approximately 1.0 microg of Esherichia coli genomic DNA, or 2 x 10(8) copies of the target gene, and only E. coli DNA was detected with the microarray assay using municipal raw sewage. Sensitivity of the microarray was enhanced approximately by 6 orders of magnitude when the target 23S rRNA gene sequences were PCR amplified with a novel universal primer set and allowed hybridization to 24 species-specific oligonucleotide probes. The minimum detection limit was estimated to be about 100 fg of E. coli genomic DNA or 1.4 x 10(2) copies of the 23S rRNA gene. The PCR amplified DNA microarray successfully detected multiple bacterial pathogens in wastewater. As a parallel study to verify efficiency of the DNA microarray, a real-time quantitative PCR assay was also developed based on the fluorescent TaqMan probes (Applied Biosystems).     

GPMiner: an integrated system for mining combinatorial cis-regulatory elements in mammalian gene group

Background

Massively parallel sequencing technology is revolutionizing approaches to genomic and genetic research. Since its advent, the scale and efficiency of Next-Generation Sequencing (NGS) has rapidly improved. In spite of this success, sequencing genomes or genomic regions with extremely biased base composition is still a great challenge to the currently available NGS platforms. The genomes of some important pathogenic organisms like Plasmodium falciparum (high AT content) and Mycobacterium tuberculosis (high GC content) display extremes of base composition. The standard library preparation procedures that employ PCR amplification have been shown to cause uneven read coverage particularly across AT and GC rich regions, leading to problems in genome assembly and variation analyses. Alternative library-preparation approaches that omit PCR amplification require large quantities of starting material and hence are not suitable for small amounts of DNA/RNA such as those from clinical isolates. We have developed and optimized library-preparation procedures suitable for low quantity starting material and tolerant to extremely high AT content sequences.

Results

We have used our optimized conditions in parallel with standard methods to prepare Illumina sequencing libraries from a non-clinical and a clinical isolate (containing ~53% host contamination). By analyzing and comparing the quality of sequence data generated, we show that our optimized conditions that involve a PCR additive (TMAC), produces amplified libraries with improved coverage of extremely AT-rich regions and reduced bias toward GC neutral templates.

Conclusion

We have developed a robust and optimized Next-Generation Sequencing library amplification method suitable for extremely AT-rich genomes. The new amplification conditions significantly reduce bias and retain the complexity of either extremes of base composition. This development will greatly benefit sequencing clinical samples that often require amplification due to low mass of DNA starting material.     

Endonuclease-mediated long PCR and its application to restriction mapping

The polymerase chain reaction (PCR) is the most widely used technique for the study of DNA. Applications for PCR have been extended significantly by the development of "long" PCR, a technique that makes it possible to amplify DNA fragments up to 40 kb in length. This article describes two novel applications of the long PCR technique, one which simplifies restriction mapping and another which enhances amplification specificity and yield. The same primers used to perform the long PCR amplification can be used as probes to perform restriction mapping of the DNA fragment amplified. Restriction digestion performed prior to long PCR amplification can be used to selectively suppress the amplification of members of families of closely related DNA sequences, thereby making it possible to selectively amplify one of a group of highly homologous sequences. These two complimentary techniques, both involving use of the long PCR paired with restriction digestion, have potential application in any laboratory in which PCR is performed.     

Mung Bean Nuclease Treatment Increases Capture Specificity of Microdroplet-PCR Based Targeted DNA Enrichment

Targeted DNA enrichment coupled with next generation sequencing has been increasingly used for interrogation of select sub-genomic regions at high depth of coverage in a cost effective manner. Specificity measured by on-target efficiency is a key performance metric for target enrichment. Non-specific capture leads to off-target reads, resulting in waste of sequencing throughput on irrelevant regions. Microdroplet-PCR allows simultaneous amplification of up to thousands of regions in the genome and is among the most commonly used strategies for target enrichment. Here we show that carryover of single-stranded template genomic DNA from microdroplet-PCR constitutes a major contributing factor for off-target reads in the resultant libraries. Moreover, treatment of microdroplet-PCR enrichment products with a nuclease specific to single-stranded DNA alleviates off-target load and improves enrichment specificity. We propose that nuclease treatment of enrichment products should be incorporated in the workflow of targeted sequencing using microdroplet-PCR for target capture. These findings may have a broad impact on other PCR based applications for which removal of template DNA is beneficial.