Gene capture and random amplification for quantitative recovery of homologous genes

A comprehensive benchmarking study of protocols and sequencing platforms for 16S rRNA community profiling

Background

In the last 5 years, the rapid pace of innovations and improvements in sequencing technologies has completely changed the landscape of metagenomic and metagenetic experiments. Therefore, it is critical to benchmark the various methodologies for interrogating the composition of microbial communities, so that we can assess their strengths and limitations. The most common phylogenetic marker for microbial community diversity studies is the 16S ribosomal RNA gene and in the last 10 years the field has moved from sequencing a small number of amplicons and samples to more complex studies where thousands of samples and multiple different gene regions are interrogated.

Results

We assembled 2 synthetic communities with an even (EM) and uneven (UM) distribution of archaeal and bacterial strains and species, as metagenomic control material, to assess performance of different experimental strategies. The 2 synthetic communities were used in this study, to highlight the limitations and the advantages of the leading sequencing platforms: MiSeq (Illumina), The Pacific Biosciences RSII, 454 GS-FLX/+ (Roche), and IonTorrent (Life Technologies). We describe an extensive survey based on synthetic communities using 3 experimental designs (fusion primers, universal tailed tag, ligated adaptors) across the 9 hypervariable 16S rDNA regions. We demonstrate that library preparation methodology can affect data interpretation due to different error and chimera rates generated during the procedure. The observed community composition was always biased, to a degree that depended on the platform, sequenced region and primer choice. However, crucially, our analysis suggests that 16S rRNA sequencing is still quantitative, in that relative changes in abundance of taxa between samples can be recovered, despite these biases.

Conclusion

We have assessed a range of experimental conditions across several next generation sequencing platforms using the most up-to-date configurations. We propose that the choice of sequencing platform and experimental design needs to be taken into consideration in the early stage of a project by running a small trial consisting of several hypervariable regions to quantify the discriminatory power of each region. We also suggest that the use of a synthetic community as a positive control would be beneficial to identify the potential biases and procedural drawbacks that may lead to data misinterpretation. The results of this study will serve as a guideline for making decisions on which experimental condition and sequencing platform to consider to achieve the best microbial profiling.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2194-9) contains supplementary material, which is available to authorized users.

Deconstructing the polymerase chain reaction: understanding and correcting bias associated with primer degeneracies and primer-template mismatches

The polymerase chain reaction (PCR) is sensitive to mismatches between primer and template, and mismatches can lead to inefficient amplification of targeted regions of DNA template. In PCRs in which a degenerate primer pool is employed, each primer can behave differently. Therefore, inefficiencies due to different primer melting temperatures within a degenerate primer pool, in addition to mismatches between primer binding sites and primers, can lead to a distortion of the true relative abundance of targets in the original DNA pool. A theoretical analysis indicated that a combination of primer-template and primer-amplicon interactions during PCR cycles 3–12 is potentially responsible for this distortion. To test this hypothesis, we developed a novel amplification strategy, entitled “Polymerase-exonuclease (PEX) PCR”, in which primer-template interactions and primer-amplicon interactions are separated. The PEX PCR method substantially and significantly improved the evenness of recovery of sequences from a mock community of known composition, and allowed for amplification of templates with introduced mismatches near the 3’ end of the primer annealing sites. When the PEX PCR method was applied to genomic DNA extracted from complex environmental samples, a significant shift in the observed microbial community was detected. Furthermore, the PEX PCR method provides a mechanism to identify which primers in a primer pool are annealing to target gDNA. Primer utilization patterns revealed that at high annealing temperatures in the PEX PCR method, perfect match annealing predominates, while at lower annealing temperatures, primers with up to four mismatches with templates can contribute substantially to amplification. The PEX PCR method is simple to perform, is limited to PCR mixes and a single exonuclease step which can be performed without reaction cleanup, and is recommended for reactions in which degenerate primer pools are used or when mismatches between primers and template are possible.

5-color multiplexed microwave-accelerated metal-enhanced fluorescence: detection and analysis of multiple DNA sequences from within one sample well within a few seconds

We present a potentially highly sensitive and selective bio-assay for the potential detection of any five different DNA sequences from one sample in one well. The assay is based on a DNA "rapid catch and signal" (DNA-RCS) technology developed for the detection of different DNA sequences from a sample well area. Our signal amplification utilizes the metal-enhanced fluorescence (MEF) of dyes attached to the probe-DNAs, which hybridizes with the pre-formed mixture of anchor-DNA scaffolds on silver island films (SiFs). Low-power microwave irradiation accelerates both the formation of the anchor-DNA scaffold on the SiF-surface and anchor/probe DNA hybridization, i.e. "rapid catch" of target DNAs from a bulk solution, decreasing the assay run time from hours to only a few seconds. Localization of signaling dye-labels close to the SiFs make them extremely photostable, which allows for collecting/integrating the signal over a long time period. To demonstrate a 5 color DNA assay (5-plex) we have used a range of readily available Alexa™ dyes. Advantages and perspectives of the RCS-technologies ability to detect 5 different DNA sequences from within one plate-well are discussed.

Targeted amplicon sequencing (TAS): a scalable next-gen approach to multilocus, multitaxa phylogenetics

Next-gen sequencing technologies have revolutionized data collection in genetic studies and advanced genome biology to novel frontiers. However, to date, next-gen technologies have been used principally for whole genome sequencing and transcriptome sequencing. Yet many questions in population genetics and systematics rely on sequencing specific genes of known function or diversity levels. Here, we describe a targeted amplicon sequencing (TAS) approach capitalizing on next-gen capacity to sequence large numbers of targeted gene regions from a large number of samples. Our TAS approach is easily scalable, simple in execution, neither time-nor labor-intensive, relatively inexpensive, and can be applied to a broad diversity of organisms and/or genes. Our TAS approach includes a bioinformatic application, BarcodeCrucher, to take raw next-gen sequence reads and perform quality control checks and convert the data into FASTA format organized by gene and sample, ready for phylogenetic analyses. We demonstrate our approach by sequencing targeted genes of known phylogenetic utility to estimate a phylogeny for the Pancrustacea. We generated data from 44 taxa using 68 different 10-bp multiplexing identifiers. The overall quality of data produced was robust and was informative for phylogeny estimation. The potential for this method to produce copious amounts of data from a single 454 plate (e.g., 325 taxa for 24 loci) significantly reduces sequencing expenses incurred from traditional Sanger sequencing. We further discuss the advantages and disadvantages of this method, while offering suggestions to enhance the approach.

Strategy for extracting DNA from clay soil and detecting a specific target sequence via selective enrichment and real-time (quantitative) PCR amplification

We present a simple strategy for isolating and accurately enumerating target DNA from high-clay-content soils: desorption with buffers, an optional magnetic capture hybridization step, and quantitation via real-time PCR. With the developed technique, microg quantities of DNA were extracted from mg samples of pure kaolinite and a field clay soil.

Real-time PCR detection and quantification of nine potential sources of fecal contamination by analysis of mitochondrial cytochrome b targets

We designed and tested real-time PCR probe/primer sets to detect and quantify Cytochrome b sequences of mitochondrial DNA (mtDNA) from nine vertebrate species of pet (dog), farm (cow, chicken, sheep, horse, pig), wildlife (Canada goose, white-tailed deer), and human. Linear ranges of the assays were from 10(1) to 10(8) copies/microl. To formally test the performance of the assays, twenty blinded fecal suspension samples were analyzed by real-time PCR to identify the source of the feces. Sixteen of the twenty samples were correctly and unambiguously identified. Average sensitivity was calculated to be 0.850, while average specificity was found to be 0.994. One beef cow sample was not detected, but mtDNA from 11 other beef cattle of both sexes and varying physiological states was found in concentrations similar (3.45 x 10(7) copies/g) to thatfound in human feces (1.1 x 10(7) copies/g). Thus, environmental conditions and sample handling are probably important factors for successful detection of fecal mtDNA. When sewage samples were analyzed, only human mtDNA (7.2 x 10(4) copies/ 100 mL) was detected. With a detection threshold of 250 copies/reaction, an efficient concentration and purification method resulted in a final detection limit for human feces of 1.8 mg/ 100 mL water.