Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

A robust PCR primer design platform applied to the detection of Acidobacteria Group 1 in soil

Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33,684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO(2). Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.     

Evaluating Bias of Illumina-Based Bacterial 16S rRNA Gene Profiles

Massively parallel sequencing of 16S rRNA genes enables the comparison of terrestrial, aquatic, and host-associated microbial communities with sufficient sequencing depth for robust assessments of both alpha and beta diversity. Establishing standardized protocols for the analysis of microbial communities is dependent on increasing the reproducibility of PCR-based molecular surveys by minimizing sources of methodological bias. In this study, we tested the effects of template concentration, pooling of PCR amplicons, and sample preparation/interlane sequencing on the reproducibility associated with paired-end Illumina sequencing of bacterial 16S rRNA genes. Using DNA extracts from soil and fecal samples as templates, we sequenced pooled amplicons and individual reactions for both high (5- to 10-ng) and low (0.1-ng) template concentrations. In addition, all experimental manipulations were repeated on two separate days and sequenced on two different Illumina MiSeq lanes. Although within-sample sequence profiles were highly consistent, template concentration had a significant impact on sample profile variability for most samples. Pooling of multiple PCR amplicons, sample preparation, and interlane variability did not influence sample sequence data significantly. This systematic analysis underlines the importance of optimizing template concentration in order to minimize variability in microbial-community surveys and indicates that the practice of pooling multiple PCR amplicons prior to sequencing contributes proportionally less to reducing bias in 16S rRNA gene surveys with next-generation sequencing.     

New Insights into the Diversity of Marine Picoeukaryotes

Over the last decade, culture-independent surveys of marine picoeukaryotic diversity based on 18S ribosomal DNA clone libraries have unveiled numerous sequences of novel high-rank taxa. This newfound diversity has significantly altered our understanding of marine microbial food webs and the evolution of eukaryotes. However, the current picture of marine eukaryotic biodiversity may be significantly skewed by PCR amplification biases, occurrence of rDNA genes in multiple copies within a single cell, and the capacity of DNA to persist as extracellular material. In this study we performed an analysis of the metagenomic dataset from the Global Ocean Survey (GOS) expedition, seeking eukaryotic ribosomal signatures. This PCR-free approach revealed similar phylogenetic patterns to clone library surveys, suggesting that PCR steps do not impose major biases in the exploration of environmental DNA. The different cell size fractions within the GOS dataset, however, displayed a distinct picture. High protistan diversity in the <0.8 µm size fraction, in particular sequences from radiolarians and ciliates (and their absence in the 0.8–3 µm fraction), suggest that most of the DNA in this fraction comes from extracellular material from larger cells. In addition, we compared the phylogenetic patterns from rDNA and reverse transcribed rRNA 18S clone libraries from the same sample harvested in the Mediterranean Sea. The libraries revealed major differences, with taxa such as pelagophytes or picobiliphytes only detected in the 18S rRNA library. MAST (Marine Stramenopiles) appeared as potentially prominent grazers and we observed a significant decrease in the contribution of alveolate and radiolarian sequences, which overwhelmingly dominated rDNA libraries. The rRNA approach appears to be less affected by taxon-specific rDNA copy number and likely better depicts the biogeochemical significance of marine protists.     

Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients.

A group-specific primer, F243 (positions 226 to 243, Escherichia coli numbering), was developed by comparison of sequences of genes encoding 16S rRNA (16S rDNA) for the detection of actinomycetes in the environment with PCR and temperature or denaturing gradient gel electrophoresis (TGGE or DGGE, respectively). The specificity of the forward primer in combination with different reverse ones was tested with genomic DNA from a variety of bacterial strains. Most actinomycetes investigated could be separated by TGGE and DGGE, with both techniques giving similar results. Two strategies were employed to study natural microbial communities. First, we used the selective amplification of actinomycete sequences (E. coli positions 226 to 528) for direct analysis of the products in denaturing gradients. Second, a nested PCR providing actinomycete-specific fragments (E. coli positions 226 to 1401) was used which served as template for a PCR when conserved primers were used. The products (E. coli positions 968 to 1401) of this indirect approach were then separated by use of gradient gels. Both approaches allowed detection of actinomycete communities in soil. The second strategy allowed the estimation of the relative abundance of actinomycetes within the bacterial community. Mixtures of PCR-derived 16S rDNA fragments were used as model communities consisting of five actinomycetes and five other bacterial species. Actinomycete products were obtained over a 100-fold dilution range of the actinomycete DNA in the model community by specific PCR; detection of the diluted actinomycete DNA was not possible when conserved primers were used. The methods tested for detection were applied to monitor actinomycete community changes in potato rhizosphere and to investigate actinomycete diversity in different soils.     

Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.     

Succession of microbial communities during hot composting as detected by PCR-single-strand-conformation polymorphism-based genetic profiles of small-subunit rRNA genes

A cultivation-independent technique for genetic profiling of PCR-amplified small-subunit rRNA genes (SSU rDNA) was chosen to characterize the diversity and succession of microbial communities during composting of an organic agricultural substrate. PCR amplifications were performed with DNA directly extracted from compost samples and with primers targeting either (i) the V4-V5 region of eubacterial 16S rRNA genes, (ii) the V3 region in the 16S rRNA genes of actinomycetes, or (iii) the V8-V9 region of fungal 18S rRNA genes. Homologous PCR products were converted to single-stranded DNA molecules by exonuclease digestion and were subsequently electrophoretically separated by their single-strand-conformation polymorphism (SSCP). Genetic profiles obtained by this technique showed a succession and increasing diversity of microbial populations with all primers. A total of 19 single products were isolated from the profiles by PCR reamplification and cloning. DNA sequencing of these molecular isolates showed similarities in the range of 92.3 to 100% to known gram-positive bacteria with a low or high G+C DNA content and to the SSU rDNA of gamma-Proteobacteria. The amplified 18S rRNA gene sequences were related to the respective gene regions of Candida krusei and Candida tropicalis. Specific molecular isolates could be attributed to different composting stages. The diversity of cultivated bacteria isolated from samples taken at the end of the composting process was low. A total of 290 isolates were related to only 6 different species. Two or three of these species were also detectable in the SSCP community profiles. Our study indicates that community SSCP profiles can be highly useful for the monitoring of bacterial diversity and community successions in a biotechnologically relevant process.     

Bacterial diversity in an industrial wastewater bioreactor

Industrial wastewater bioreactors are potentially important sources of novel biocatalysts. However, the microbial populations in these bioreactors are not well characterized. The microbial community in an industrial wastewater bioreactor was surveyed by extracting DNA from a sample of activated sludge, followed by PCR amplification and sequencing of cloned 16S rRNA genes. A total of 407 cloned 16S rRNA gene sequences were compared with 88 bacterial isolates cultured from the same sample of sludge using a variety of standard media. Most of the bacteria detected by the PCR-based approach were -subdivision Proteobacteria, whereas most of the cultured bacteria were -subdivision Proteobacteria. Only a few types of bacteria were detected by both approaches. These observations indicate that multiple techniques are necessary to characterize the microbial diversity in any complex ecosystem.     

Genotypic characterization of soil bacteria in the Umm Al-Namil Island,Kuwait

Microflora is an integral part of soil ecosystem, in which bacteria are the largest group of soil microbes. This is a pioneer study for establishing baseline data on the diversity of soil bacteria among different regions in Kuwait. The aim is to understand biodiversity in different settings, how bacteria adapt to different niches in the environment as well as in different hosts. The identification of bacterial 16S rRNA molecules from environmental soil samples was investigated. Genomic Deoxyribonucleic acid DNA was extracted from 25 soil samples derived from five different test regions in the Umm Al-Namil Island, Kuwait. After amplification of bacterial 16S rRNA molecules by the Polymerase chain reaction PCR, the products were characterized and complex band patterns were obtained, indicating high bacterial diversity. A sample of the 16 s rRNA amplicons were sequenced in order to identify the species. The spatial distribution of bacterial taxa in the different soil samples was homogeneous, suggesting a stable and widespread community. Forty-nine isolates from Umm Al-Namil island were identified by comparative analysis of partial 16S rRNA gene sequences. Phylogenetic analysis was carried out in order to study the connection between the isolates to identify species. A large proportion of these isolates represent correspond to known or novel species within the Pseudomonus and Bacillus genera, which are common soil bacteria. Our results provided a reference for future studies to facilitate bacterial identification and ecological research in Kuwait.     

Influence of Commonly Used Primer Systems on Automated Ribosomal Intergenic Spacer Analysis of Bacterial Communities in Environmental Samples

Due to the high diversity of bacteria in many ecosystems, their slow generation times, specific but mostly unknown nutrient requirements and syntrophic interactions, isolation based approaches in microbial ecology mostly fail to describe microbial community structure. Thus, cultivation independent techniques, which rely on directly extracted nucleic acids from the environment, are a well-used alternative. For example, bacterial automated ribosomal intergenic spacer analysis (B-ARISA) is one of the widely used methods for fingerprinting bacterial communities after PCR-based amplification of selected regions of the operon coding for rRNA genes using community DNA. However, B-ARISA alone does not provide any taxonomic information and the results may be severely biased in relation to the primer set selection. Furthermore, amplified DNA stemming from mitochondrial or chloroplast templates might strongly bias the obtained fingerprints. In this study, we determined the applicability of three different B-ARISA primer sets to the study of bacterial communities. The results from in silico analysis harnessing publicly available sequence databases showed that all three primer sets tested are specific to bacteria but only two primers sets assure high bacterial taxa coverage (1406f/23Sr and ITSF/ITSReub). Considering the study of bacteria in a plant interface, the primer set ITSF/ITSReub was found to amplify (in silico) sequences of some important crop species such as Sorghum bicolor and Zea mays. Bacterial genera and plant species potentially amplified by different primer sets are given. These data were confirmed when DNA extracted from soil and plant samples were analyzed. The presented information could be useful when interpreting existing B-ARISA results and planning B-ARISA experiments, especially when plant DNA can be expected.     

Diversity of bacteria associated with grassland soil nematodes of different feeding groups

The bacterial diversity associated with soil nematodes and its relationship with their feeding habits are as yet poorly understood. In the present study the diversity and abundance of bacteria from nematodes and their surrounding soil were analysed and compared. The nematodes were collected from a grassland soil and sorted into bacterial, fungal, plant, predatory and omnivore feeding groups and assigned to taxonomic groups. Total DNA was extracted from the nematodes and partial bacterial 16S rRNA genes were PCR amplified, cloned and sequenced. The abundance and composition of bacterial taxa differed between and within feeding groups. The lowest bacterial diversity was found in the predatory nematodes Prionchulus sp., whereas the highest bacterial diversity was associated with the bacterial-feeding nematode Acrobeles sp. The soil had a more diverse bacterial community than the communities found in the nematode groups. The 16S rRNA gene sequences of bacteria associated with nematodes did not overlap with those detected in soil as determined using the cloning screening approach. However, bacterial sequences identified from nematodes could be detected in the soil with targeted PCR. Our data suggest that the nematodes do not feed on the most abundant bacteria present in soil. Furthermore, several nematodes contained suspected bacterial symbionts and parasites.     

Bacterial Diversity of the Broadbalk ‘Classical’ Winter Wheat Experiment in Relation to Long-Term Fertilizer Inputs

With more than 160 years of contrasting fertilizer regimes, the Broadbalk winter wheat experiment represents a unique experimental resource for studying the effects of long-term fertilizer application on microbial population diversity. Using DGGE and clone library analysis, we report here on eubacterial species diversity (16S rRNA gene) and diversity within two sets of gene products associated with microbial N acquisition: NifH (nitrogen fixation) and AmtB (ammonium transport). Comparisons were made within and between soils treated with mineral N fertilizer, farmyard manure or receiving no fertilizer. Analysis of 16S rRNA gene DGGE profiles showed no clear patterns to qualitatively distinguish bacterial community structure between the three different treatments (P > 0.05), with all samples containing a range of eubacterial taxa similar to those that are characteristic of soil bacteria reported elsewhere. Intra-plot heterogeneity was high and of a similar magnitude to that between treatments. This lack of qualitative between plot differences was echoed in the representative sequences of 16S rRNA, nifH, and amtB genes in the various samples. Taken together, both phylogenetic and functional gene analyses showed bacterial communities in the Broadbalk-trial soil were very stable and relatively non-responsive to long-term management of balanced fertilizer inputs.     

Diversity of polyketide synthase genes from bacteria associated with the marine sponge Pseudoceratina clavata: culture-dependent and culture-independent approaches

Diverse ketosynthase (KS) genes were retrieved from the microbial community associated with the Great Barrier Reef sponge Pseudoceratina clavata. Bacterial isolation and metagenomic approaches were employed. Phylogenetic analysis of 16S rRNA of culturable sponge-associated bacterial communities comprised eight groups over four phyla. Ten KS domains were amplified from four genera of isolates and phylogenetics demonstrated that these KS domains were located in three clusters (actinobacterial, cyanobacterial and trans-AT type). Metagenomic DNA of the sponge microbial community was extracted to explore community KS genes by two approaches: direct amplification of KS domains and construction of fosmid libraries for KS domain screening. Five KS domains were retrieved from polymerase chain reaction (PCR) amplification using sponge metagenome DNA as template and five fosmid clones containing KS domains found using multiplex PCR screening. Analysis of selected polyketide synthase (PKS) from one fosmid showed that the PKS consists of two modules. Open reading frames located up- and downstream of the PKS displayed similarity with membrane synthesis-related proteins such as cardiolipin synthase. Metagenome approaches did not detect KS domains found in sponge isolates. All KS domains from both metagenome approaches formed a single cluster with KS domains originating from metagenomes derived from other sponge species from other geographical regions.     

Diversity of bacteria and polyketide synthase associated with marine sponge Haliclona sp.

The microbial community associated with a marine sponge (Haliclona sp.) collected from Tateyama city, Japan was studied using 16S rRNA gene clone libraries. Two DNA templates were prepared using methods recommended for Gram-positive and Gram-negative bacteria in the Qiagen kit manual. From each DNA template, two 16S rRNA genes were PCR amplified, using the combination of universal bacterial primer 27f and primers 1385r and 1492r, respectively. A total of 347 clones were sequenced and compared with those available in DNA data banks. These sequences were members of ten bacterial phyla. Interestingly, more than 30 % of the clones represent novel sequences. A comparison of these sequences with sequences in a library prepared from DNA extracted from the surrounding water shows minimum DNA contamination. Taxonomically, the highest diversity was detected in the clone library prepared using a combination of primers 27f and 1492r and DNA isolated using the Gram-positive bacteria protocol. The potential of Haliclona sp.-associated bacteria to produce secondary metabolites was studied by cloning and sequencing the polyketide synthase (PKS, type 1) gene using the same DNA samples. Analysis of partial sequences derived from the sponge metagenome revealed 27 unique ketosynthase domains of PKS type I. This study suggests strongly that this Haliclona sp. plays host to diverse novel bacteria with a potential to produce novel polyketides.     

Estimating Bacterial Diversity for Ecological Studies: Methods,Metrics, and Assumptions

Methods to estimate microbial diversity have developed rapidly in an effort to understand the distribution and diversity of microorganisms in natural environments. For bacterial communities, the 16S rRNA gene is the phylogenetic marker gene of choice, but most studies select only a specific region of the 16S rRNA to estimate bacterial diversity. Whereas biases derived from from DNA extraction, primer choice and PCR amplification are well documented, we here address how the choice of variable region can influence a wide range of standard ecological metrics, such as species richness, phylogenetic diversity, β-diversity and rank-abundance distributions. We have used Illumina paired-end sequencing to estimate the bacterial diversity of 20 natural lakes across Switzerland derived from three trimmed variable 16S rRNA regions (V3, V4, V5). Species richness, phylogenetic diversity, community composition, β-diversity, and rank-abundance distributions differed significantly between 16S rRNA regions. Overall, patterns of diversity quantified by the V3 and V5 regions were more similar to one another than those assessed by the V4 region. Similar results were obtained when analyzing the datasets with different sequence similarity thresholds used during sequences clustering and when the same analysis was used on a reference dataset of sequences from the Greengenes database. In addition we also measured species richness from the same lake samples using ARISA Fingerprinting, but did not find a strong relationship between species richness estimated by Illumina and ARISA. We conclude that the selection of 16S rRNA region significantly influences the estimation of bacterial diversity and species distributions and that caution is warranted when comparing data from different variable regions as well as when using different sequencing techniques.     

Use of rpoB and 16S rRNA genes to analyse bacterial diversity of a tropical soil using PCR and DGGE

AIM: To evaluate the rpoB gene as a biomarker for PCR-DGGE microbial analyses using soil DNA from the Cerrado, Brazil. METHODS: DNA extraction from soil was followed by Polymerase Chain Reaction (PCR) amplification of rpoB and 16S rRNA genes. PCR products were compared by Denaturing Gradient Gel Electrophoresis (DGGE) to compare gene/community profiles. RESULTS: The rpoB DGGE profiles comprised fewer bands than the 16S rDNA profiles and were easier to delineate and therefore to analyse. Comparison of the community profiles revealed that the methods were complementary. CONCLUSIONS, SIGNIFICANCE AND IMPACT OF THE STUDY: The gene for the beta subunit of the RNA polymerase, rpoB, is a single copy gene unlike 16S rDNA. Multiple copies of 16S rRNA genes in bacterial genomes complicate diversity assessments made from DGGE profiles. Using the rpoB gene offers a better alternative to the commonly used 16S rRNA gene for microbial community analyses based on DGGE.     

应用rpoB和16S rDNA基因的变性梯度凝胶电泳技术对山羊瘤胃细菌多样性的研究

采用免培养的rpoB和16S rDNA基因的变性梯度凝胶电泳技术(DGGE)对3种山羊(波尔山羊,内蒙古绒山羊,四川南江黄羊)瘤胃细菌优势菌群结构进行了比较分析。研究结果显示rpoBDGGE图谱中条带数目少于16S rDNA图谱,并且条带分离效果明显,更有利于分析瘤胃细菌群落组成。从两种DGGE图谱中均可以发现3种山羊瘤胃细菌具有一定的相似性,种内个体间相似性明显高于种间相似性,这说明寄主品种是影响瘤胃细菌种群构成的一个重要因素。同时进行了部分优势细菌16S rDNA基因V6-V8区序列的系统发育分析。基因序列分析表明,DGGE图谱中优势条带的16S rDNA基因序列中有4条克隆的序列与基因库最相似菌的相似性大于97%,余下的克隆序列相似性在89%～96%之间,其中13条序列的与之相似性最高的序列均来自于未被鉴定的瘤胃细菌。     

Influence of an oyster reef on development of the microbial heterotrophic community of an estuarine biofilm

We characterized microbial biofilm communities developed over two very closely located but distinct benthic habitats in the Pensacola Bay estuary using two complementary cultivation-independent molecular techniques. Biofilms were grown for 7 days on glass slides held in racks 10 to 15 cm over an oyster reef and an adjacent muddy sand bottom. Total biomass and optical densities of dried biofilms showed dramatic differences for oyster reef versus non-oyster reef biofilms. This study assessed whether the observed spatial variation was reflected in the heterotrophic prokaryotic species composition. Genomic biofilm DNA from both locations was isolated and served as a template to amplify 16S rRNA genes with universal eubacterial primers. Fluorescently labeled PCR products were analyzed by terminal restriction fragment length polymorphism, creating a genetic fingerprint of the composition of the microbial communities. Unlabeled PCR products were cloned in order to construct a clone library of 16S rRNA genes. Amplified ribosomal DNA restriction analysis was used to screen and define ribotypes. Partial sequences from unique ribotypes were compared with existing database entries to identify species and to construct phylogenetic trees representative of community structures. A pronounced difference in species richness and evenness was observed at the two sites. The biofilm community structure from the oyster reef setting had greater evenness and species richness than the one from the muddy sand bottom. The vast majority of the bacteria in the oyster reef biofilm were related to members of the gamma- and delta-subdivisions of Proteobacteria, the Cytophaga-Flavobacterium -Bacteroides cluster, and the phyla Planctomyces and Holophaga-Acidobacterium. The same groups were also present in the biofilm harvested at the muddy sand bottom, with the difference that nearly half of the community consisted of representatives of the Planctomyces phylum. Total species richness was estimated to be 417 for the oyster reef and 60 for the muddy sand bottom, with 10.5% of the total unique species identified being shared between habitats. The results suggest dramatic differences in habitat-specific microbial diversity that have implications for overall microbial diversity within estuaries.