Reevaluation and reduction of a PCR bias caused by reannealing of templates

Multi-factorial examination of amplicon sequencing workflows from sample preparation to bioinformatic analysis

BACKGROUND

The development of sequencing technologies to evaluate bacterial microbiota composition has allowed new insights into the importance of microbial ecology. However, the variety of methodologies used among amplicon sequencing workflows leads to uncertainty about best practices as well as reproducibility and replicability among microbiome studies. Using a bacterial mock community composed of 37 soil isolates, we performed a comprehensive methodological evaluation of workflows, each with a different combination of methodological factors spanning sample preparation to bioinformatic analysis to define sources of artifacts that affect coverage, accuracy, and biases in the resulting compositional profiles.

RESULTS

Of the workflows examined, those using the V4-V4 primer set enabled the highest level of concordance between the original mock community and resulting microbiome sequence composition. Use of a high-fidelity polymerase, or a lower-fidelity polymerase with an increased PCR elongation time, limited chimera formation. Bioinformatic pipelines presented a trade-off between the fraction of distinct community members identified (coverage) and fraction of correct sequences (accuracy). DADA2 and QIIME2 assembled V4-V4 reads amplified by Taq polymerase resulted in the highest accuracy (100%) but had a coverage of only 52%. Using mothur to assemble and denoise V4-V4 reads resulted in a coverage of 75%, albeit with marginally lower accuracy (99.5%).

CONCLUSIONS

Optimization of microbiome workflows is critical for accuracy and to support reproducibility and replicability among microbiome studies. These considerations will help reveal the guiding principles of microbial ecology and impact the translation of microbiome research to human and environmental health.

Coexistence patterns of soil methanogens are closely tied to methane generation and community assembly in rice paddies

BACKGROUND

Soil methanogens participate in complex interactions, which determine the community structures and functions. Studies continue to seek the coexistence patterns of soil methanogens, influencing factors and the contribution to methane (CH4) production, which are regulated primarily by species interactions, and the functional significance of these interactions. Here, methane emissions were measured in rice paddies across the Asian continent, and the complex interactions involved in coexistence patterns of methanogenic archaeal communities were represented as pairwise links in co-occurrence networks.

RESULTS

The network topological properties, which were positively correlated with mean annual temperature, were the most important predictor of CH4 emissions among all the biotic and abiotic factors. The methanogenic groups involved in commonly co-occurring links among the 39 local networks contributed most to CH4 emission (53.3%), much higher than the contribution of methanogenic groups with endemic links (36.8%). The potential keystone taxa, belonging to Methanobacterium, Methanocella, Methanothrix, and Methanosarcina, possessed high linkages with the methane generation functional genes mcrA, fwdB, mtbA, and mtbC. Moreover, the commonly coexisting taxa showed a very different assembly pattern, with ~ 30% determinism and ~ 70% stochasticity. In contrast, a higher proportion of stochasticity (93~99%) characterized the assembly of endemically coexisting taxa.

CONCLUSIONS

These results suggest that the coexistence patterns of microbes are closely tied to their functional significance, and the potential importance of common coexistence further imply that complex networks of interactions may contribute more than species diversity to soil functions. Video abstract.

Patterns of PCR Amplification Artifacts of the Fungal Barcode Marker in a Hybrid Mushroom

The polymerase chain reaction (PCR) is widely used in modern biology and medicine. However, PCR artifacts can complicate the interpretation of PCR-based results. The internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster is the consensus fungal barcode marker and suspected PCR artifacts have been reported in many studies, especially for the analyses of environmental fungal samples. At present, the patterns of PCR artifacts in the whole fungal ITS region (ITS1+5.8S+ITS2) are not known. In this study, we analyzed the error rates of PCR at three template complexity levels using the divergent copies of ITS from the mushroom Agaricus subrufescens. Our results showed that PCR using the Phusion^® High-Fidelity DNA Polymerase has a per nucleotide error rate of about 4 × 10^–6 per replication. Among the detected mutations, transitions were much more frequent than transversions, insertions, and deletions. When divergent alleles were mixed as templates in the same reaction, a significant proportion (∼30%) of recombinant molecules were detected. The in vitro mixed-template results were comparable to those obtained from using the genomic DNA of the original mushroom specimen as template. Our results indicate that caution should be in place when interpreting ITS sequences from individual fungal specimens, especially those containing divergent ITS copies. Similar results could also happen to PCR-based analyses of other multicopy DNA fragments as well as single-copy DNA sequences with divergent alleles in diploid organisms.

Impact of salinity on the population dynamics of microorganisms in a membrane bioreactor treating produced water

Biomass characteristics may change subject to a salinity increase when treating high strength wastewater. In this study, the impact of salinity in a membrane bioreactor (MBR) was investigated for the treatment of produced water (PW). MBR was operated as a pre-treatment prior to nanofiltration (NF) and reverse osmosis (RO). Mixed PW, that was originated from different oil, gas, and oil-gas wells, were subjected to pre-treatment for 146 days including three different operational phases: Low-salinity (~10 mS/cm), gradual increased salinity (10-40 mS/cm) and high salinity (~40 mS/cm). The results of this study showed that microorganisms could adapt using real wastewater and treat PW up to a certain level of the salinity. At high salinity levels of PW, the floc structure started to be disrupted and membrane fouling was accelerated. PCR-DGGE analysis confirmed the changes in microbial communities' composition in relation with high salinity. The results of the final treatment experiments presented that NF and RO treatment produced high-quality effluents that could be suitable for reuse.

Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis

The purpose of this review is to present the most common and emerging DNA-based methods used to generate data for biodiversity and biomonitoring studies. As environmental assessment and monitoring programmes may require biodiversity information at multiple levels, we pay particular attention to the DNA metabarcoding method and discuss a number of bioinformatic tools and considerations for producing DNA-based indicators using operational taxonomic units (OTUs), taxa at a variety of ranks and community composition. By developing the capacity to harness the advantages provided by the newest technologies, investigators can "scale up" by increasing the number of samples and replicates processed, the frequency of sampling over time and space, and even the depth of sampling such as by sequencing more reads per sample or more markers per sample. The ability to scale up is made possible by the reduced hands-on time and cost per sample provided by the newest kits, platforms and software tools. Results gleaned from broad-scale monitoring will provide opportunities to address key scientific questions linked to biodiversity and its dynamics across time and space as well as being more relevant for policymakers, enabling science-based decision-making, and provide a greater socio-economic impact. As genomic approaches are continually evolving, we provide this guide to methods used in biodiversity genomics.

High-frequency, low-coverage "false positives" mutations may be true in GS Junior sequencing studies

The GS Junior sequencer provides simplified procedures for library preparation and data processing. Errors in pyrosequencing generate some biases during library construction and emulsion PCR amplification. False-positive mutations are identified by related characteristics described in the manufacturer’s manual, and some detected mutations may have ‘borderline’ characteristics when they are detected in few reads or at low frequency. Among these mutations, however, some may be true positives. This study aimed to improve the accuracy of identifying true positives among mutations with borderline false-positive characteristics detected with GS Junior sequencing. Mutations with the borderline features were tested for validity with Sanger sequencing. We examined 10 mutations detected in coverages <20-fold at frequencies >30% (group A) and 16 mutations detected in coverages >20-fold at frequencies < 30% (group B). In group A, two mutations were not confirmed, and two mutations with 100% frequency were confirmed as heterozygous alleles. No mutation in group B was confirmed. The two groups had significantly different false-positive prevalences (p = 0.001). These results suggest that mutations detected at frequencies less than 30% can be confidently identified as false-positives but that mutations detected at frequencies over 30%, despite coverages less than 20-fold, should be verified with Sanger sequencing.

The Use of Degenerate Primers in qPCR Analysis of Functional Genes Can Cause Dramatic Quantification Bias as Revealed by Investigation of nifH Primer Performance

The measurement of functional gene abundance in diverse microbial communities often employs quantitative PCR (qPCR) with highly degenerate oligonucleotide primers. While degenerate PCR primers have been demonstrated to cause template-specific bias in PCR applications, the effect of such bias on qPCR has been less well explored. We used a set of diverse, full-length nifH gene standards to test the performance of several universal nifH primer sets in qPCR. We found significant template-specific bias in all but the PolF/PolR primer set. Template-specific bias caused more than 1000-fold mis-estimation of nifH gene copy number for three of the primer sets and one primer set resulted in more than 10,000-fold mis-estimation. Furthermore, such template-specific bias will cause qPCR estimates to vary in response to beta-diversity, thereby causing mis-estimation of changes in gene copy number. A reduction in bias was achieved by increasing the primer concentration. We conclude that degenerate primers should be evaluated across a range of templates, annealing temperatures, and primer concentrations to evaluate the potential for template-specific bias prior to their use in qPCR.

Molecular Analysis of the Cyanobacterial Community in Gastric Contents of Egrets with Symptoms of Steatitis

Many deaths of wild birds that have drunk water contaminated with hepatotoxic microcystin-producing cyanobacteria have been reported. A mass death of egrets and herons with steatitis were found at the agricultural reservoir occurring cyanobacterial waterblooms. This study aimed to verify a hypothesis that the egrets and herons which died in the reservoir drink microcystin-producing cyanobacteria and microcystin involves in the cause of death as well as the symptoms of steatitis. The cyanobacterial community in gastric contents of egrets and herons that died from steatitis was assessed using cyanobacterial 16S rRNA-based terminal-restriction fragment length polymorphism (T-RFLP) profiling and a cyanobacterial 16S rRNA-based clone library analysis. In addition, PCR amplification of the mcyB-C region and the mcyG gene, involved in microcystin biosynthesis, was examined. The cyanobacterial community in the gastric contents of two birds showed a simplistic composition. A comparison of cyanobacterial T-RFLP profiling and cloned sequences suggested that the genus Microcystis predominated in both samples of egrets died. Although we confirmed that two egrets which died in the reservoir have taken in cyanobacterial waterblooms containing the genus Microcystis, no mcy gene was detected in both samples according to the mcy gene-based PCR analysis. This study is the first to show the profiling and traceability of a cyanobacterial community in the gastric contents of wild birds by molecular analysis. Additionally, we consider causing symptoms of steatitis in the dead egrets.

Assessment of a 16S rRNA amplicon Illumina sequencing procedure for studying the microbiome of a symbiont-rich aphid genus

The bacterial communities inhabiting arthropods are generally dominated by a few endosymbionts that play an important role in the ecology of their hosts. Rather than comparing bacterial species richness across samples, ecological studies on arthropod endosymbionts often seek to identify the main bacterial strains associated with each specimen studied. The filtering out of contaminants from the results and the accurate taxonomic assignment of sequences are therefore crucial in arthropod microbiome studies. We aimed here to validate an Illumina 16S rRNA gene sequencing protocol and analytical pipeline for investigating endosymbiotic bacteria associated with aphids. Using replicate DNA samples from 12 species (Aphididae: Lachninae, Cinara) and several controls, we removed individual sequences not meeting a minimum threshold number of reads in each sample and carried out taxonomic assignment for the remaining sequences. With this approach, we show that (i) contaminants accounted for a negligible proportion of the bacteria identified in our samples; (ii) the taxonomic composition of our samples and the relative abundance of reads assigned to a taxon were very similar across PCR and DNA replicates for each aphid sample; in particular, bacterial DNA concentration had no impact on the results. Furthermore, by analysing the distribution of unique sequences across samples rather than aggregating them into operational taxonomic units (OTUs), we gained insight into the specificity of endosymbionts for their hosts. Our results confirm that Serratia symbiotica is often present in Cinara species, in addition to the primary symbiont, Buchnera aphidicola. Furthermore, our findings reveal new symbiotic associations with Erwinia- and Sodalis-related bacteria. We conclude with suggestions for generating and analysing 16S rRNA gene sequences for arthropod-endosymbiont studies.

Computational analysis of stochastic heterogeneity in PCR amplification efficiency revealed by single molecule barcoding

The polymerase chain reaction (PCR) is one of the most widely used techniques in molecular biology. In combination with High Throughput Sequencing (HTS), PCR is widely used to quantify transcript abundance for RNA-seq, and in the context of analysis of T and B cell receptor repertoires. In this study, we combine DNA barcoding with HTS to quantify PCR output from individual target molecules. We develop computational tools that simulate both the PCR branching process itself, and the subsequent subsampling which typically occurs during HTS sequencing. We explore the influence of different types of heterogeneity on sequencing output, and compare them to experimental results where the efficiency of amplification is measured by barcodes uniquely identifying each molecule of starting template. Our results demonstrate that the PCR process introduces substantial amplification heterogeneity, independent of primer sequence and bulk experimental conditions. This heterogeneity can be attributed both to inherited differences between different template DNA molecules, and the inherent stochasticity of the PCR process. The results demonstrate that PCR heterogeneity arises even when reaction and substrate conditions are kept as constant as possible, and therefore single molecule barcoding is essential in order to derive reproducible quantitative results from any protocol combining PCR with HTS.

Primer and platform effects on 16S rRNA tag sequencing

Sequencing of 16S rRNA gene tags is a popular method for profiling and comparing microbial communities. The protocols and methods used, however, vary considerably with regard to amplification primers, sequencing primers, sequencing technologies; as well as quality filtering and clustering. How results are affected by these choices, and whether data produced with different protocols can be meaningfully compared, is often unknown. Here we compare results obtained using three different amplification primer sets (targeting V4, V6–V8, and V7–V8) and two sequencing technologies (454 pyrosequencing and Illumina MiSeq) using DNA from a mock community containing a known number of species as well as complex environmental samples whose PCR-independent profiles were estimated using shotgun sequencing. We find that paired-end MiSeq reads produce higher quality data and enabled the use of more aggressive quality control parameters over 454, resulting in a higher retention rate of high quality reads for downstream data analysis. While primer choice considerably influences quantitative abundance estimations, sequencing platform has relatively minor effects when matched primers are used. Beta diversity metrics are surprisingly robust to both primer and sequencing platform biases.

Improved pipeline for reducing erroneous identification by 16S rRNA sequences using the Illumina MiSeq platform

The cost of DNA sequencing has decreased due to advancements in Next Generation Sequencing. The number of sequences obtained from the Illumina platform is large, use of this platform can reduce costs more than the 454 pyrosequencer. However, the Illumina platform has other challenges, including bioinformatics analysis of large numbers of sequences and the need to reduce erroneous nucleotides generated at the 3'-ends of the sequences. These erroneous sequences can lead to errors in analysis of microbial communities. Therefore, correction of these erroneous sequences is necessary for accurate taxonomic identification. Several studies that have used the Illumina platform to perform metagenomic analyses proposed curating pipelines to increase accuracy. In this study, we evaluated the likelihood of obtaining an erroneous microbial composition using the MiSeq 250 bp paired sequence platform and improved the pipeline to reduce erroneous identifications. We compared different sequencing conditions by varying the percentage of control phiX added, the concentration of the sequencing library, and the 16S rRNA gene target region using a mock community sample composed of known sequences. Our recommended method corrected erroneous nucleotides and improved identification accuracy. Overall, 99.5% of the total reads shared 95% similarity with the corresponding template sequences and 93.6% of the total reads shared over 97% similarity. This indicated that the MiSeq platform can be used to analyze microbial communities at the genus level with high accuracy. The improved analysis method recommended in this study can be applied to amplicon studies in various environments using high-throughput reads generated on the MiSeq platform.

Bias in Whole Genome Amplification: Causes and Considerations

Whole genome amplification (WGA) is a widely used molecular technique that is becoming increasingly necessary in genetic research on a range of sample types including individual cells, fossilized remains and entire ecosystems. Multiple methods of WGA have been developed, each with specific strengths and weaknesses, but with a common defect in that each method distorts the initial template DNA during the course of amplification. The type, extent, and circumstance of the bias vary with the WGA method and particulars of the template DNA. In this review, we endeavor to discuss the types of bias introduced, the susceptibility of common WGA techniques to these bias types, and the interdependence between bias and characteristics of the template DNA. Finally, we attempt to illustrate some of the criteria specific to the analytical platform and research application that should be considered to enable combination of the appropriate WGA method, template DNA, sequencing platform, and intended use for optimal results.

Description of an unusual Neisseria meningitidis isolate containing and expressing Neisseria gonorrhoeae-Specific 16S rRNA gene sequences

An apparently rare Neisseria meningitidis isolate containing one copy of a Neisseria gonorrhoeae 16S rRNA gene is described herein. This isolate was identified as N. meningitidis by biochemical identification methods but generated a positive signal with Gen-Probe Aptima assays for the detection of Neisseria gonorrhoeae. Direct 16S rRNA gene sequencing of the purified isolate revealed mixed bases in signature regions that allow for discrimination between N. meningitidis and N. gonorrhoeae. The mixed bases were resolved by sequencing individually PCR-amplified single copies of the genomic 16S rRNA gene. A total of 121 discrete sequences were obtained; 92 (76%) were N. meningitidis sequences, and 29 (24%) were N. gonorrhoeae sequences. Based on the ratio of species-specific sequences, the N. meningitidis strain seems to have replaced one of its four intrinsic 16S rRNA genes with the gonococcal gene. Fluorescence in situ hybridization (FISH) probes specific for meningococcal and gonococcal rRNA were used to demonstrate the expression of the rRNA genes. Interestingly, the clinical isolate described here expresses both N. meningitidis and N. gonorrhoeae 16S rRNA genes, as shown by positive FISH signals with both probes. This explains why the probes for N. gonorrhoeae in the Gen-Probe Aptima assays cross-react with this N. meningitidis isolate. The N. meningitidis isolate described must have obtained N. gonorrhoeae-specific DNA through interspecies recombination.

The existence and diversity of myxobacteria in lake mud - a previously unexplored myxobacteria habitat

Myxobacteria are widely distributed in soil and oceanic sediment with a phylogeographic separation at high levels of classification. However, it is unclear whether freshwater environments, from which there has been no isolation report of myxobacteria since 1981, are habitats for myxobacteria. In this study, we investigated the presence of myxobacteria in lake mud using a two-step strategy. First, we constructed two universal bacterial libraries from the V3-V4 (V34) and V6-V8 (V678) hypervariable regions of 16S rRNA gene sequences. High-throughput 454 pyrosequencing revealed that myxobacteria were one of the major bacterial groups in the lake mud. They accounted for 5.77% of the total sequences and 7.52% of the total operational taxonomic units (OTUs) at a phylogenetic distance of 0.03. The community composition and taxonomic structure of the mud myxobacterial community were further analysed using myxobacteria-enriched libraries targeting the V34 and V678 regions, which were amplified with Cystobacterineae- and Sorangineae-specific primer pairs respectively. Phylogenetic analysis showed that the limnetic myxobacteria exhibited closer relationships to their soil than their marine relatives, but there were also exclusive taxa of limnetic myxobacteria detected. These results, together with a survey on available GenBank data, indicate that lake mud is a primary habitat for myxobacteria.

Archaeal Diversity of Upland Rice Field Soils Assessed by the Terminal Restriction Fragment Length Polymorphism Method Combined with Real Time Quantitative-PCR and a Clone Library Analysis

The PCR amplification-based analysis of microbial diversity is subject to potential problems. In this study, to minimize the bias toward a 1:1 ratio in multitemplate PCR, a real-time PCR assay was carried out using a quenching fluorescence dye primer and amplification efficiency was monitored. Then terminal-restriction fragment length polymorphism (T-RFLP) profiling was performed using the PCR product with minimized PCR bias. This method was applied to an analysis of the diversity of the archaeal community in an upland rice field under different tillage systems and winter cover cropping. Terminal restriction fragments (T-RFs) of PCR-amplified archaeal 16S rRNA genes were assigned to the gene sequences recovered from the same soil by using an archaeal 16S rRNA gene clone library. Our results indicated that soil archaeal members were not influenced but the relative abundance of archaeal species particularly those belonging to Crenarchaeota which changed between the tillage and non-tillage treatments.

Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

Analysis of microbial communities by high-throughput pyrosequencing of SSU rRNA gene PCR amplicons has transformed microbial ecology research and led to the observation that many communities contain a diverse assortment of rare taxa-a phenomenon termed the Rare Biosphere. Multiple studies have investigated the effect of pyrosequencing read quality on operational taxonomic unit (OTU) richness for contrived communities, yet there is limited information on the fidelity of community structure estimates obtained through this approach. Given that PCR biases are widely recognized, and further unknown biases may arise from the sequencing process itself, a priori assumptions about the neutrality of the data generation process are at best unvalidated. Furthermore, post-sequencing quality control algorithms have not been explicitly evaluated for the accuracy of recovered representative sequences and its impact on downstream analyses, reducing useful discussion on pyrosequencing reads to their diversity and abundances. Here we report on community structures and sequences recovered for in vitro-simulated communities consisting of twenty 16S rRNA gene clones tiered at known proportions. PCR amplicon libraries of the V3-V4 and V6 hypervariable regions from the in vitro-simulated communities were sequenced using the Roche 454 GS FLX Titanium platform. Commonly used quality control protocols resulted in the formation of OTUs with >1% abundance composed entirely of erroneous sequences, while over-aggressive clustering approaches obfuscated real, expected OTUs. The pyrosequencing process itself did not appear to impose significant biases on overall community structure estimates, although the detection limit for rare taxa may be affected by PCR amplicon size and quality control approach employed. Meanwhile, PCR biases associated with the initial amplicon generation may impose greater distortions in the observed community structure.

Massively parallel haplotyping on microscopic beads for the high-throughput phase analysis of single molecules

In spite of the many advances in haplotyping methods, it is still very difficult to characterize rare haplotypes in tissues and different environmental samples or to accurately assess the haplotype diversity in large mixtures. This would require a haplotyping method capable of analyzing the phase of single molecules with an unprecedented throughput. Here we describe such a haplotyping method capable of analyzing in parallel hundreds of thousands single molecules in one experiment. In this method, multiple PCR reactions amplify different polymorphic regions of a single DNA molecule on a magnetic bead compartmentalized in an emulsion drop. The allelic states of the amplified polymorphisms are identified with fluorescently labeled probes that are then decoded from images taken of the arrayed beads by a microscope. This method can evaluate the phase of up to 3 polymorphisms separated by up to 5 kilobases in hundreds of thousands single molecules. We tested the sensitivity of the method by measuring the number of mutant haplotypes synthesized by four different commercially available enzymes: Phusion, Platinum Taq, Titanium Taq, and Phire. The digital nature of the method makes it highly sensitive to detecting haplotype ratios of less than 1:10,000. We also accurately quantified chimera formation during the exponential phase of PCR by different DNA polymerases.

Simple absolute quantification method correcting for quantitative PCR efficiency variations for microbial community samples

Real-time quantitative PCR (qPCR) is a widely used technique in microbial community analysis, allowing the quantification of the number of target genes in a community sample. Currently, the standard-curve (SC) method of absolute quantification is widely employed for these kinds of analysis. However, the SC method assumes that the amplification efficiency (E) is the same for both the standard and the sample target template. We analyzed 19 bacterial strains and nine environmental samples in qPCR assays, targeting the nifH and 16S rRNA genes. The E values of the qPCRs differed significantly, depending on the template. This has major implications for the quantification. If the sample and standard differ in their E values, quantification errors of up to orders of magnitude are possible. To address this problem, we propose and test the one-point calibration (OPC) method for absolute quantification. The OPC method corrects for differences in E and was derived from the ΔΔC(T) method with correction for E, which is commonly used for relative quantification in gene expression studies. The SC and OPC methods were compared by quantifying artificial template mixtures from Geobacter sulfurreducens (DSM 12127) and Nostoc commune (Culture Collection of Algae and Protozoa [CCAP] 1453/33), which differ in their E values. While the SC method deviated from the expected nifH gene copy number by 3- to 5-fold, the OPC method quantified the template mixtures with high accuracy. Moreover, analyzing environmental samples, we show that even small differences in E between the standard and the sample can cause significant differences between the copy numbers calculated by the SC and the OPC methods.

Soil microbial community successional patterns during forest ecosystem restoration

Soil microbial community characterization is increasingly being used to determine the responses of soils to stress and disturbances and to assess ecosystem sustainability. However, there is little experimental evidence to indicate that predictable patterns in microbial community structure or composition occur during secondary succession or ecosystem restoration. This study utilized a chronosequence of developing jarrah (Eucalyptus marginata) forest ecosystems, rehabilitated after bauxite mining (up to 18 years old), to examine changes in soil bacterial and fungal community structures (by automated ribosomal intergenic spacer analysis [ARISA]) and changes in specific soil bacterial phyla by 16S rRNA gene microarray analysis. This study demonstrated that mining in these ecosystems significantly altered soil bacterial and fungal community structures. The hypothesis that the soil microbial community structures would become more similar to those of the surrounding nonmined forest with rehabilitation age was broadly supported by shifts in the bacterial but not the fungal community. Microarray analysis enabled the identification of clear successional trends in the bacterial community at the phylum level and supported the finding of an increase in similarity to nonmined forest soil with rehabilitation age. Changes in soil microbial community structure were significantly related to the size of the microbial biomass as well as numerous edaphic variables (including pH and C, N, and P nutrient concentrations). These findings suggest that soil bacterial community dynamics follow a pattern in developing ecosystems that may be predictable and can be conceptualized as providing an integrated assessment of numerous edaphic variables.

Mathematical modeling of 16S ribosomal DNA amplification reveals optimal conditions for the interrogation of complex microbial communities with phylogenetic microarrays

MOTIVATION

Many current studies of complex microbial communities rely on the isolation of community genomic DNA, amplification of 16S ribosomal RNA genes (rDNA) and subsequent examination of community structure through interrogation of the amplified 16S rDNA pool by high-throughput sequencing, phylogenetic microarrays or quantitative PCR.

RESULTS

Here we describe the development of a mathematical model aimed to simulate multitemplate amplification of 16S ribosomal DNA sample and subsequent detection of these amplified 16S rDNA species by phylogenetic microarray. Using parameters estimated from the experimental results obtained in the analysis of intestinal microbial communities with Microbiota Array, we show that both species detection and the accuracy of species abundance estimates depended heavily on the number of PCR cycles used to amplify 16S rDNA. Both parameters initially improved with each additional PCR cycle and reached optimum between 15 and 20 cycles of amplification. The use of more than 20 cycles of PCR amplification and/or more than 50 ng of starting genomic DNA template was, however, detrimental to both the fraction of detected community members and the accuracy of abundance estimates. Overall, the outcomes of the model simulations matched well available experimental data. Our simulations also showed that species detection and the accuracy of abundance measurements correlated positively with the higher sample-wide PCR amplification rate, lower template-to-template PCR bias and lower number of species in the interrogated community. The developed model can be easily modified to simulate other multitemplate DNA mixtures as well as other microarray designs and PCR amplification protocols.

Optimizing the analysis of human intestinal microbiota with phylogenetic microarray

Phylogenetic microarrays present an attractive strategy to high-throughput interrogation of complex microbial communities. In this work, we present several approaches to optimize the analysis of intestinal microbiota with the recently developed Microbiota Array. First, we determined how 16S rDNA-specific PCR amplification influenced bacterial detection and the consistency of measured abundance values. Bacterial detection improved with an increase in the number of PCR amplification cycles, but 25 cycles were sufficient to achieve the maximum possible detection. A PCR-caused deviation in the measured abundance values was also observed. We also developed two mathematical algorithms that aimed to account for a predicted cross-hybridization of 16S rDNA fragments among different species, and to adjust the measured hybridization signal based on the number of 16S rRNA gene copies per species genome. The 16S rRNA gene copy adjustment indicated that the presence of members of the class Clostridia might be overestimated in some 16S rDNA-based studies. Finally, we show that the examination of total community RNA with phylogenetic microarray can provide estimates of the relative metabolic activity of individual community members. Complementary profiling of genomic DNA and total RNA isolated from the same sample presents an opportunity to assess population structure and activity in the same microbial community.

Nested PCR bias: a case study of Pseudomonas spp. in soil microcosms

Nested PCR methods combined with denaturing gradient gel electrophoresis (DGGE) are widely used for the detection of low copy number species or for the analysis of group-specific community profiles. With an appropriate number of PCR cycles during the first round of amplification, initial differences in the copy numbers of different DNA fragments that are targeted can be maintained during the second round without significant bias. However, if an excessive number of cycles in used in the first round, relative differences in the copy numbers of the targeted sequences can be obscured. Here we demonstrate the effect of "nested PCR bias" in a case study with PCR-DGGE of 16S rRNA gene sequences targeting Pseudomonas spp. following exposure of soil to naphthalene vapors. Our results demonstrate artifacts caused by nested PCR bias can be substantially minimized by calibrating the number of first round PCR cycles, thereby preserving the ability to obtain semiquantitative data for evaluating changes in gene copy numbers over time.

Influence of long-term land application of Class B biosolids on soil bacterial diversity

AIM

To evaluate the effect of long-term annual land applications of Class B biosolids on soil bacterial diversity at University of Arizona Marana Agricultural Field Center, Tucson, Arizona.

METHODS AND RESULTS

Following the final of 20 consecutive years of application of Class B biosolids in March 2005, followed by cotton growth from April to November 2005 surface soil samples (0-30 cm) were collected from control (unamended) and biosolid-amended plots. Total bacterial community DNA was extracted, amplified using 16S rRNA primers, cloned, and sequenced. All 16S rRNA sequences were identified by 16S rRNA sequence analysis and comparison to known sequences in GenBank (NCBI BlastN and Ribosomal Database Project II, RDP). Results showed that the number of known genera (identifiable > 96%) increased in the high rate biosolid plots compared to control plots. Biosolids-amended soils had a broad phylogenetic diversity comprising more than four major phyla: Proteobacteria (32%), Acidobacteria (21%), Actinobacteria (16%), Firmicutes (7%), and Bacteroidetes (6%) which were typical to bacterial diversity found in the unamended arid southwestern soils.

CONCLUSION

Bacterial diversity was either enhanced or was not negatively impacted following 20 years of land application of Class B biosolids.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study illustrates that long-term land application of biosolids to arid southwestern desert soils has no deleterious effect on soil microbial diversity.

Reducing the impact of PCR-mediated recombination in molecular evolution and environmental studies using a new-generation high-fidelity DNA polymerase

PCR-mediated recombination can greatly impact estimates of diversity, both in environmental studies and in analyses of gene family evolution. Here we measure chimera (PCR-mediated recombinant) formation by analyzing a mixture of eight partial actin sequences isolated from the amoeba Arcella hemisphaerica amplified under a variety of conditions that mimic standard laboratory situations. We further compare a new-generation proofreading processivity-enhanced polymerase to both a standard proofreading enzyme and previously published results. Proofreading polymerases are preferred over other polymerases in instances where evolutionary inferences must be made. Our analyses reveal that reducing the initial template concentration is as critical as reducing the number of cycles for decreasing chimera formation and improving accuracy. Furthermore, assessing the efficiency of recovery of original haplotypes demonstrates that multiple PCR reactions are required to capture the actual genetic diversity of a sample. Finally, the experiments confirm that processivity-enhanced polymerases enable a substantial decrease in PCR-mediated recombination through reducing starting template concentration, without compromising the robustness of PCR reactions.

Rapid estimation of microbial populations in fish samples by using terminal restriction fragment length polymorphism analysis of 16S rDNA

A rapid system using terminal restriction fragment length polymorphism (T-RFLP) analysis targeting 16S rDNA is described for microbial population analysis in edible fish samples. The defined terminal restriction fragment database was constructed by collecting 102 strains of bacteria representing 53 genera that are associated with fish. Digestion of these 102 strains with two restriction enzymes, HhaI and MspI, formed 54 pattern groups with discrimination to the genus level. This T-RFLP system produced results comparable to those from a culture-based method in six natural fish samples with a qualitative correspondence of 71.4 to 92.3%. Using the T-RFLP system allowed an estimation of the microbial population within 7 h. Rapid assay of the microbial population is advantageous for food manufacturers and testing laboratories; moreover, the strategy presented here allows adaptation to specific testing applications.

Polymerase chain reaction primers miss half of rRNA microbial diversity

The rRNA approach is the principal tool to study microbial diversity, but it has important biases. These include polymerase chain reaction (PCR) primers bias, and relative inefficiency of DNA extraction techniques. Such sources of potential undersampling of microbial diversity are well known, but the scale of the undersampling has not been quantified. Using a marine tidal flat bacterial community as a model, we show that even with unlimited sampling and sequencing effort, a single combination of PCR primers/DNA extraction technique enables theoretical recovery of only half of the richness recoverable with three such combinations. This shows that different combinations of PCR primers/DNA extraction techniques recover in principle different species, as well as higher taxa. The majority of earlier estimates of microbial richness seem to be underestimates. The combined use of multiple PCR primer sets, multiple DNA extraction techniques, and deep community sequencing will minimize the biases and recover substantially more species than prior studies, but we caution that even this--yet to be used--approach may still leave an unknown number of species and higher taxa undetected.

Chemically modified primers for improved multiplex polymerase chain reaction

Multiplex polymerase chain reaction (PCR), the amplification of multiple targets in a single reaction, presents a new set of challenges that further complicate more traditional PCR setups. These complications include a greater probability for nonspecific amplicon formation and for imbalanced amplification of different targets, each of which can compromise quantification and detection of multiple targets. Despite these difficulties, multiplex PCR is frequently used in applications such as pathogen detection, RNA quantification, mutation analysis, and (recently) next generation DNA sequencing. Here we investigated the utility of primers with one or two thermolabile 4-oxo-1-pentyl phosphotriester modifications in improving multiplex PCR performance. Initial endpoint and real-time analyses revealed a decrease in off-target amplification and a subsequent increase in amplicon yield. Furthermore, the use of modified primers in multiplex setups revealed a greater limit of detection and more uniform amplification of each target as compared with unmodified primers. Overall, the thermolabile modified primers present a novel and exciting avenue for improving multiplex PCR performance.

High-throughput quantitative analysis of the human intestinal microbiota with a phylogenetic microarray

Gut microbiota carry out key functions in health and participate in the pathogenesis of a growing number of diseases. The aim of this study was to develop a custom microarray that is able to identify hundreds of intestinal bacterial species. We used the Entrez nucleotide database to compile a data set of bacterial 16S rRNA gene sequences isolated from human intestinal and fecal samples. Identified sequences were clustered into separate phylospecies groups. Representative sequences from each phylospecies were used to develop a microbiota microarray based on the Affymetrix GeneChip platform. The designed microbiota array contains probes to 775 different bacterial phylospecies. In our validation experiments, the array correctly identified genomic DNA from all 15 bacterial species used. Microbiota array has a detection sensitivity of at least 1 pg of genomic DNA and can detect bacteria present at a 0.00025% level of overall sample. Using the developed microarray, fecal samples from two healthy children and two healthy adults were analyzed for bacterial presence. Between 227 and 232 species were detected in fecal samples from children, whereas 191 to 208 species were found in adult stools. The majority of identified phylospecies belonged to the classes Clostridia and Bacteroidetes. The microarray revealed putative differences between the gut microbiota of healthy children and adults: fecal samples from adults had more Clostridia and less Bacteroidetes and Proteobacteria than those from children. A number of other putative differences were found at the genus level.

Quantitative analysis of methanogenic community dynamics in three anaerobic batch digesters treating different wastewaters

Quantitative changes in methanogenic community structures, associated with performance data, were investigated in three anaerobic batch digesters treating synthetic glucose medium, whey permeate, and liquefied sewage sludge. All digesters were initially seeded with anaerobic sludge obtained from a local municipal wastewater treatment plant. Dynamics of methanogenic populations were monitored, at order and family levels, using real-time PCR based on the 16S rRNA gene. The molecular monitoring revealed that, in each digester, the quantitative structure of methanogenic community varied continuously over treatment time and the variation corresponded well to the changes in chemical profiles. Biphasic production of methane, associated with successive increases in aceticlastic (mainly Methanosarcinaceae) and hydrogenotrophic (mainly Methanomicrobiales) methanogenic groups, was observed in each digester. This corresponded to the diauxic utilization of acetate and longer-chain volatile fatty acids (C(3)-C(6)), mainly propionate. Additionally, the non-metric multidimensional scaling (NMDS) analysis of the quantification results demonstrated that the community shift patterns in three digesters were totally different from each other. Considering that the operating conditions in all trials were identical except substrates, the differences in quantitative shift profiles were suggested to be due to the different substrate compositions. This implied that the composition of wastewater could affect the evolution of quantitative methanogenic community structure in an anaerobic process. Overall, our results suggested that more attention to quantitative as well as qualitative approaches on microbial communities is needed for fundamental understanding of anaerobic processes, particularly under dynamic or transitional conditions.

Environmental rRNA inventories miss over half of protistan diversity

BACKGROUND

The main tool to discover novel microbial eukaryotes is the rRNA approach. This approach has important biases, including PCR discrimination against certain rRNA gene species, which makes molecular inventories skewed relative to the source communities. The degree of this bias has not been quantified, and it remains unclear whether species missed from clone libraries could be recovered by increasing sequencing efforts, or whether they cannot be detected in principle. Here we attempt to discriminate between these possibilities by statistically analysing four protistan inventories obtained using different general eukaryotic PCR primers.

RESULTS

We show that each PCR primer set-specific clone library is not a sample from the community diversity but rather from a fraction of this diversity. Therefore, even sequencing such clone libraries to saturation would only recover that fraction, which, according to the parametric models, varies between 17 +/- 4% to 49 +/- 10%, depending on the set of primers. The pooled data is thus qualitatively richer than individual libraries, even if normalized to the same sequencing effort.

CONCLUSION

The use of a single pair of primers leads to significant underestimation of the true community richness at all levels of taxonomic hierarchy. The majority of available protistan rRNA gene surveys likely sampled less than half of the target diversity, and might have completely missed the rest. The use of multiple PCR primers reduces this bias but does not necessarily eliminate it.

Rapid phylogenetic dissection of prokaryotic community structure in tidal flat using pyrosequencing

Dissection of prokaryotic community structure is prerequisite to understand their ecological roles. Various methods are available for such a purpose which amplification and sequencing of 16S rRNA genes gained its popularity. However, conventional methods based on Sanger sequencing technique require cloning process prior to sequencing, and are expensive and labor-intensive. We investigated prokaryotic community structure in tidal flat sediments, Korea, using pyrosequencing and a subsequent automated bioinformatic pipeline for the rapid and accurate taxonomic assignment of each amplicon. The combination of pyrosequencing and bioinformatic analysis showed that bacterial and archaeal communities were more diverse than previously reported in clone library studies. Pyrosequencing analysis revealed 21 bacterial divisions and 37 candidate divisions. Proteobacteria was the most abundant division in the bacterial community, of which Gamma-and Delta-Proteobacteria were the most abundant. Similarly, 4 archaeal divisions were found in tidal flat sediments. Euryarchaeota was the most abundant division in the archaeal sequences, which were further divided into 8 classes and 11 unclassified euryarchaeota groups. The system developed here provides a simple, in-depth and automated way of dissecting a prokaryotic community structure without extensive pretreatment such as cloning.

Advances in the use of terminal restriction fragment length polymorphism (T-RFLP) analysis of 16S rRNA genes to characterize microbial communities

Terminal restriction fragment length polymorphism (T-RFLP) analysis is a popular high-throughput fingerprinting technique used to monitor changes in the structure and composition of microbial communities. This approach is widely used because it offers a compromise between the information gained and labor intensity. In this review, we discuss the progress made in T-RFLP analysis of 16S rRNA genes and functional genes over the last 10 years and evaluate the performance of this technique when used in conjunction with different statistical methods. Web-based tools designed to perform virtual polymerase chain reaction and restriction enzyme digests greatly facilitate the choice of primers and restriction enzymes for T-RFLP analysis. Significant improvements have also been made in the statistical analysis of T-RFLP profiles such as the introduction of objective procedures to distinguish between signal and noise, the alignment of T-RFLP peaks between profiles, and the use of multivariate statistical methods to detect changes in the structure and composition of microbial communities due to spatial and temporal variation or treatment effects. The progress made in T-RFLP analysis of 16S rRNA and genes allows researchers to make methodological and statistical choices appropriate for the hypotheses of their studies.

Bacterial taxa abundance pattern in an industrial wastewater treatment system determined by the full rRNA cycle approach

The description of the diversity and structure of microbial communities through quantification of the constituent populations is one of the major objectives in environmental microbiology. The implications of models for community assembly are practical as well as theoretical, because the extent of biodiversity is thought to influence the function of ecosystems. Current attempts to predict species diversity in different environments derive the numbers of individuals for each operational taxonomic unit (OTU) from the frequency of clones in 16S rDNA gene libraries, which are subjected to a number of inherent biases and artefacts. We show that diversity of the bacterial community present in a complex microbial ensemble can be estimated by fitting the data of the full-cycle rRNA approach to a model of species abundance distribution. Sequences from a 16S rDNA gene library from activated sludge were reliably assigned to OTUs at a genetic distance of 0.04. A group of 17 newly designed rRNA-targeted oligonucleotide probes were used to quantify by fluorescence in situ hybridization, OTUs represented with more than three clones in the 16S rDNA clone library. Cell abundance distribution was best described by a geometric series, after the goodness of fit was evaluated by the Kolmogorov-Smirnov test. Although a complete mechanistic understanding of all the ecological processes involved is still not feasible, describing the distribution pattern of a complex bacterial assemblage model can shed light on the way bacterial communities operate.

Reliability for detecting composition and changes of microbial communities by T-RFLP genetic profiling

Terminal restriction fragment length polymorphism (T-RFLP) analysis is commonly used for profiling microbial communities in various environments. However, it may suffer from biases during the analytic process. This study addressed the potential of T-RFLP profiles (1) to reflect real community structures and diversities, as well as (2) to reliably detect changing components of microbial community structures. For this purpose, defined artificial communities of 30 SSU rRNA gene clones, derived from nine bacterial phyla, were used. PCR amplification efficiency was one primary bias with a maximum variability factor of 3.5 among clones. PCR downstream analyses such as enzymatic restriction and capillary electrophoresis introduced a maximum bias factor of 4 to terminal restriction fragment (T-RF) signal intensities, resulting in a total maximum bias factor of 14 in the final T-RFLP profiles. In addition, the quotient between amplification efficiency and T-RF size allowed predicting T-RF abundances in the profiles with high accuracy. Although these biases impaired detection of real community structures, the relative changes in structures and diversities were reliably reflected in the T-RFLP profiles. These data support the suitability of T-RFLP profiling for monitoring effects on microbial communities.

Influence of high temperature and humidity on rumen bacterial diversity in Holstein heifers

The effect of heat and humidity stresses on the rumen bacterial molecular diversity of heifers was studied. No statistically significant changes in the rumen microbiota composition were found in the first experiment (average body mass 250kg) while in the second and third experiments (additional variables included the relative humidity and body weight), the microbiota composition was significantly different at elevated environmental temperatures and humidity. These shifts were accompanied by the decrease in concentration of short-chain fatty acids in the rumen.

Distinct protistan assemblages characterize the euphotic zone and deep sea (2500 m) of the western North Atlantic (Sargasso Sea and Gulf Stream)

Protistan diversity was characterized at three locations in the western North Atlantic (Sargasso Sea and Gulf Stream) by sequencing 18S rRNA genes in samples from euphotic (< or = 125 m) and bathypelagic depths (2500 m). A total of 923 partial-length protistan sequences were analysed, revealing 324 distinct operational taxonomic units (OTUs) determined by an automated OTU-calling program set to 95% sequence similarity. Most OTUs were comprised of only one or two sequences suggesting a large but rare pool of protistan diversity. Many OTUs from both depth strata were associated with recently described novel alveolate and stramenopile lineages while many OTUs from the bathypelagic were affiliated with Acantharea, Polycystinea and Euglenozoa and were not observed in euphotic zone libraries. Protistan assemblages from the euphotic zone and the deep sea were largely composed of distinct OTUs; only 28 of the 324 protistan OTUs were detected in both shallow and deep sea clone libraries. The diversity of protistan assemblages in the deep sea was distinctly lower than the diversity of euphotic zone assemblages. Protistan assemblages from the Gulf Stream were the most diverse for either depth strata. Overall, protistan assemblages from different stations but comparable depths were more similar than the assemblages from different depths at the same station. These data suggest that particular groups of protistan OTUs formed distinct 'shallow' and 'deep-sea' assemblages across widely spaced oceanic locales.

Introduction of a boost of Legionella pneumophila into a stagnant-water model by heat treatment

An environmentally representative stagnant-water model was developed to monitor the growth dynamics of Legionella pneumophila. This model was evaluated for three distinct water treatments: untreated tap water, heat-treated tap water, and heat-treated tap water supplemented with Pseudomonas putida, a known biofilm-forming bacterium. Bringing heat-treated tap water after subsequent cooling into contact with a densely formed untreated biofilm was found to promote the number of L. pneumophila by 4 log units within the biofilm, while the use of untreated water only sustained the L. pneumophila levels. Subsequent colonization of the water phase by L. pneumophila was noticed in the heat-treated stagnant-water models, with concentrations as high as 1 x 10(10) mip gene copies L(-1) stagnant water. Denaturing gradient gel electrophoresis in combination with clustering analysis of the prokaryotic community in the water phase and in the biofilm phase suggests that the different water treatments induced different communities. Moreover, boosts of L. pneumophila arising from heat treatment of water were accompanied by shifts to a more diverse eukaryotic community. Stimulated growth of L. pneumophila after heating of the water may explain the rapid recolonization of L. pneumophila in water systems. These results highlight the need for additional or alternative measures to heat treatment of water in order to prevent or abate potential outbreaks of L. pneumophila.

Application of a high-density oligonucleotide microarray approach to study bacterial population dynamics during uranium reduction and reoxidation

Reduction of soluble uranium U(VI) to less-soluble uranium U(IV) is a promising approach to minimize migration from contaminated aquifers. It is generally assumed that, under constant reducing conditions, U(IV) is stable and immobile; however, in a previous study, we documented reoxidation of U(IV) under continuous reducing conditions (Wan et al., Environ. Sci. Technol. 2005, 39:6162-6169). To determine if changes in microbial community composition were a factor in U(IV) reoxidation, we employed a high-density phylogenetic DNA microarray (16S microarray) containing 500,000 probes to monitor changes in bacterial populations during this remediation process. Comparison of the 16S microarray with clone libraries demonstrated successful detection and classification of most clone groups. Analysis of the most dynamic groups of 16S rRNA gene amplicons detected by the 16S microarray identified five clusters of bacterial subfamilies responding in a similar manner. This approach demonstrated that amplicons of known metal-reducing bacteria such as Geothrix fermentans (confirmed by quantitative PCR) and those within the Geobacteraceae were abundant during U(VI) reduction and did not decline during the U(IV) reoxidation phase. Significantly, it appears that the observed reoxidation of uranium under reducing conditions occurred despite elevated microbial activity and the consistent presence of metal-reducing bacteria. High-density phylogenetic microarrays constitute a powerful tool, enabling the detection and monitoring of a substantial portion of the microbial population in a routine, accurate, and reproducible manner.

Assessment of whole genome amplification-induced bias through high-throughput, massively parallel whole genome sequencing

BACKGROUND

Whole genome amplification is an increasingly common technique through which minute amounts of DNA can be multiplied to generate quantities suitable for genetic testing and analysis. Questions of amplification-induced error and template bias generated by these methods have previously been addressed through either small scale (SNPs) or large scale (CGH array, FISH) methodologies. Here we utilized whole genome sequencing to assess amplification-induced bias in both coding and non-coding regions of two bacterial genomes. Halobacterium species NRC-1 DNA and Campylobacter jejuni were amplified by several common, commercially available protocols: multiple displacement amplification, primer extension pre-amplification and degenerate oligonucleotide primed PCR. The amplification-induced bias of each method was assessed by sequencing both genomes in their entirety using the 454 Sequencing System technology and comparing the results with those obtained from unamplified controls.

RESULTS

All amplification methodologies induced statistically significant bias relative to the unamplified control. For the Halobacterium species NRC-1 genome, assessed at 100 base resolution, the D-statistics from GenomiPhi-amplified material were 119 times greater than those from unamplified material, 164.0 times greater for Repli-G, 165.0 times greater for PEP-PCR and 252.0 times greater than the unamplified controls for DOP-PCR. For Campylobacter jejuni, also analyzed at 100 base resolution, the D-statistics from GenomiPhi-amplified material were 15 times greater than those from unamplified material, 19.8 times greater for Repli-G, 61.8 times greater for PEP-PCR and 220.5 times greater than the unamplified controls for DOP-PCR.

CONCLUSION

Of the amplification methodologies examined in this paper, the multiple displacement amplification products generated the least bias, and produced significantly higher yields of amplified DNA.

Development and characterization of stable sediment-free anaerobic bacterial enrichment cultures that dechlorinate aroclor 1260

We have developed sediment-free anaerobic enrichment cultures that dechlorinate a broad spectrum of highly chlorinated polychlorinated biphenyls (PCBs). The cultures were developed from Aroclor 1260-contaminated sediment from the Housatonic River in Lenox, MA. Sediment slurries were primed with 2,6-dibromobiphenyl to stimulate Process N dechlorination (primarily meta dechlorination), and sediment was gradually removed by successive transfers (10%) to minimal medium. The cultures grow on pyruvate, butyrate, or acetate plus H(2). Gas chromatography-electron capture detector analysis demonstrated that the cultures extensively dechlorinate 50 to 500 mug/ml of Aroclor 1260 at 22 to 24 degrees C by Dechlorination Process N. Triplicate cultures of the eighth transfer without sediment dechlorinated 76% of the hexa- through nonachlorobiphenyls in Aroclor 1260 (250 mug/ml) to tri- through pentachlorobiphenyls in 110 days. At least 64 PCB congeners, all of which are chlorinated on both rings and 47 of which have six or more chlorines, were substrates for this dechlorination. To characterize the bacterial diversity in the enrichments, we used eubacterial primers to amplify and clone 16S rRNA genes from DNA extracted from cultures grown on acetate plus H(2). Restriction fragment length polymorphism analysis of 107 clones demonstrated the presence of Thauera-like Betaproteobacteria, Geobacter-like Deltaproteobacteria, Pseudomonas species, various Clostridiales, Bacteroidetes, Dehalococcoides of the Chloroflexi group, and unclassified Eubacteria. Our development of highly enriched, robust, stable, sediment-free cultures that extensively dechlorinate a highly chlorinated commercial PCB mixture is a major and unprecedented breakthrough in the field. It will enable intensive study of the organisms and genes responsible for a major PCB dechlorination process that occurs in the environment and could also lead to effective remediation applications.

PCR-induced sequence artifacts and bias: insights from comparison of two 16S rRNA clone libraries constructed from the same sample

The contribution of PCR artifacts to 16S rRNA gene sequence diversity from a complex bacterioplankton sample was estimated. Taq DNA polymerase errors were found to be the dominant sequence artifact but could be constrained by clustering the sequences into 99% sequence similarity groups. Other artifacts (chimeras and heteroduplex molecules) were significantly reduced by employing modified amplification protocols. Surprisingly, no skew in sequence types was detected in the two libraries constructed from PCR products amplified for different numbers of cycles. Recommendations for modification of amplification protocols and for reporting diversity estimates at 99% sequence similarity as a standard are given.

Differentiation of important and closely related cereal plant species (Poaceae) in food by hybridization to an oligonucleotide array

We report the development of an oligonucleotide microarray for the simultaneous detection of six important cereal food plant species from the Poaceae based on the chloroplast trnL intron sequence. We used universal primers to amplify the trnL intron from wheat, rye, barley, oat, rice, and maize, followed by a cyclic labeling of oligonucleotides probes and subsequent hybridization to an oligonucleotide microarray. In single taxon analyses, positive signals were produced with a high signal-to-noise ratio. The assay also enabled the analysis of mixed samples. The results obtained for real food samples were in agreement with the ingredient labels, but positive results for grains not declared on the ingredients list were observed in three out of 10 samples, which indicates that the final products and/or the declared ingredients were probably botanically impure or contaminated. The combination of the sensitivity of a universal polymerase chain reaction with the specificity of the labeling reaction allows this protocol to be applied in routine analyses of food samples, as demonstrated by successful analysis of processed composite food products.