Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

Reconstruction of ribosomal RNA genes from metagenomic data

Direct sequencing of environmental DNA (metagenomics) has a great potential for describing the 16S rRNA gene diversity of microbial communities. However current approaches using this 16S rRNA gene information to describe community diversity suffer from low taxonomic resolution or chimera problems. Here we describe a new strategy that involves stringent assembly and data filtering to reconstruct full-length 16S rRNA genes from metagenomicpyrosequencing data. Simulations showed that reconstructed 16S rRNA genes provided a true picture of the community diversity, had minimal rates of chimera formation and gave taxonomic resolution down to genus level. The strategy was furthermore compared to PCR-based methods to determine the microbial diversity in two marine sponges. This showed that about 30% of the abundant phylotypes reconstructed from metagenomic data failed to be amplified by PCR. Our approach is readily applicable to existing metagenomic datasets and is expected to lead to the discovery of new microbial phylotypes.     

Analysis of 23S rRNA genes in metagenomes - a case study from the Global Ocean Sampling Expedition

As an evolutionary marker, 23S ribosomal RNA (rRNA) offers more diagnostic sequence stretches and greater sequence variation than 16S rRNA. However, 23S rRNA is still not as widely used. Based on 80 metagenome samples from the Global Ocean Sampling (GOS) Expedition, the usefulness and taxonomic resolution of 23S rRNA were compared to those of 16S rRNA. Since 23S rRNA is approximately twice as large as 16S rRNA, twice as many 23S rRNA gene fragments were retrieved from the GOS reads than 16S rRNA gene fragments, with 23S rRNA gene fragments being generally about 100 bp longer. Datasets for 16S and 23S rRNA sequences revealed similar relative abundances for major marine bacterial and archaeal taxa. However, 16S rRNA sequences had a better taxonomic resolution due to their significantly larger reference database.Reevaluation of the specificity of previously published PCR amplification primers and group specific fluorescence in situ hybridization probes on this metagenomic set of non-amplified 23S rRNA sequences revealed that out of 16 primers investigated, only two had more than 90% target group coverage. Evaluations of two probes, BET42a and GAM42a, were in accordance with previous evaluations, with a discrepancy in the target group coverage of the GAM42a probe when evaluated against the GOS metagenomic dataset.     

OTU Analysis Using Metagenomic Shotgun Sequencing Data

Because of technological limitations, the primer and amplification biases in targeted sequencing of 16S rRNA genes have veiled the true microbial diversity underlying environmental samples. However, the protocol of metagenomic shotgun sequencing provides 16S rRNA gene fragment data with natural immunity against the biases raised during priming and thus the potential of uncovering the true structure of microbial community by giving more accurate predictions of operational taxonomic units (OTUs). Nonetheless, the lack of statistically rigorous comparison between 16S rRNA gene fragments and other data types makes it difficult to interpret previously reported results using 16S rRNA gene fragments. Therefore, in the present work, we established a standard analysis pipeline that would help confirm if the differences in the data are true or are just due to potential technical bias. This pipeline is built by using simulated data to find optimal mapping and OTU prediction methods. The comparison between simulated datasets revealed a relationship between 16S rRNA gene fragments and full-length 16S rRNA sequences that a 16S rRNA gene fragment having a length >150 bp provides the same accuracy as a full-length 16S rRNA sequence using our proposed pipeline, which could serve as a good starting point for experimental design and making the comparison between 16S rRNA gene fragment-based and targeted 16S rRNA sequencing-based surveys possible.     

内参基因加标法定量土壤微生物目标基因绝对拷贝数

【目的】通过荧光定量PCR技术对土壤微生物目标基因进行绝对定量,其定量结果的准确性容易受到DNA提取得率以及腐殖酸抑制性的影响。【方法】采用内参基因加标法,利用构建的突变质粒DNA,对供试水稻土壤样品中的微生物16S r RNA目标基因的绝对拷贝数进行荧光定量PCR检测,用来表征该样品中细菌群落总体丰度。在定量前通过双向引物扩增方法验证突变质粒中的内参基因对供试土壤的特异性。【结果】不同水稻土壤样品的DNA提取量在样品间差异较大。通过内参基因加标法对DNA提取量进行校正,显著提高了16S r RNA基因绝对定量的精确度。不同水稻土壤样品间的变异系数为17.8,与未加标处理相比降低了66.7%。在此基础上,进一步通过内参基因加标法对土壤有机质和含水率均呈现典型空间特征差异的6处亚热带湿地土壤样品中的16S r RNA基因进行绝对定量。16S r RNA基因绝对拷贝数与土壤微生物生物量碳具有显著的线性相关性(R2=0.694,P0.001),表明内参校正后的16S r RNA基因绝对拷贝数可以准确反映单位质量土壤中微生物的丰度。【结论】内参基因加标法可以对DNA提取得率以及腐殖酸对PCR扩增的抑制性进行校正,从而提高绝对定量的准确性。基于内参基因加标法的目标基因绝对定量PCR检测,可作为土壤微生物生物量测量,以及微生物功能基因绝对丰度定量的一种核酸检测方法。     

16S rRNA基因在微生物生态学中的应用

16S rRNA(Small subunit ribosomal RNA)基因是对原核微生物进行系统进化分类研究时最常用的分子标志物(Biomarker),广泛应用于微生物生态学研究中。近些年来随着高通量测序技术及数据分析方法等的不断进步,大量基于16S rRNA基因的研究使得微生物生态学得到了快速发展,然而使用16S rRNA基因作为分子标志物时也存在诸多问题,比如水平基因转移、多拷贝的异质性、基因扩增效率的差异、数据分析方法的选择等,这些问题影响了微生物群落组成和多样性分析时的准确性。对当前使用16S rRNA基因分析微生物群落组成和多样性的进展情况做一总结,重点讨论当前存在的主要问题以及各种分析方法的发展,尤其是与高通量测序技术有关的实验和数据处理问题。     

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.     

Turkey fecal microbial community structure and functional gene diversity revealed by 16S rRNA gene and metagenomic sequences

The primary goal of this study was to better understand the microbial composition and functional genetic diversity associated with turkey fecal communities. To achieve this, 16S rRNA gene and metagenomic clone libraries were sequenced from turkey fecal samples. The analysis of 382 16S rRNA gene sequences showed that the most abundant bacteria were closely related to Lactobacillales (47%), Bacillales (31%), and Clostridiales (11%). Actinomycetales, Enterobacteriales, and Bacteroidales sequences were also identified, but represented a smaller part of the community. The analysis of 379 metagenomic sequences showed that most clones were similar to bacterial protein sequences (58%). Bacteriophage (10%) and avian viruses (3%) sequences were also represented. Of all metagenomic clones potentially encoding for bacterial proteins, most were similar to low G+C Gram-positive bacterial proteins, particularly from Lactobacillales (50%), Bacillales (11%), and Clostridiales (8%). Bioinformatic analyses suggested the presence of genes encoding for membrane proteins, lipoproteins, hydrolases, and functional genes associated with the metabolism of nitrogen and sulfur containing compounds. The results from this study further confirmed the predominance of Firmicutes in the avian gut and highlight the value of coupling 16S rRNA gene and metagenomic sequencing data analysis to study the microbial composition of avian fecal microbial communities.     

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.     

Metagenomic analysis of microbial communities across a transect from low to highly hydrocarbon-contaminated soils in King George Island,Maritime Antarctica

Soil samples from a transect from low to highly hydrocarbon-contaminated soils were collected around the Brazilian Antarctic Station Comandante Ferraz (EACF), located at King George Island, Antarctica. Quantitative PCR (qPCR) analysis of bacterial 16S rRNA genes, 16S rRNA gene (iTag), and shotgun metagenomic sequencing were used to characterize microbial community structure and the potential for petroleum degradation by indigenous microbes. Hydrocarbon contamination did not affect bacterial abundance in EACF soils (bacterial 16S rRNA gene qPCR). However, analysis of 16S rRNA gene sequences revealed a successive change in the microbial community along the pollution gradient. Microbial richness and diversity decreased with the increase of hydrocarbon concentration in EACF soils. The abundance of Cytophaga, Methyloversatilis, Polaromonas, and Williamsia was positively correlated (p-value = <.05) with the concentration of total petroleum hydrocarbons (TPH) and/or polycyclic aromatic hydrocarbons (PAH). Annotation of metagenomic data revealed that the most abundant hydrocarbon degradation pathway in EACF soils was related to alkyl derivative-PAH degradation (mainly methylnaphthalenes) via the CYP450 enzyme family. The abundance of genes related to nitrogen fixation increased in EACF soils as the concentration of hydrocarbons increased. The results obtained here are valuable for the future of bioremediation of petroleum hydrocarbon-contaminated soils in polar environments.     

Questionable 16S ribosomal RNA gene annotations are frequent in completed microbial genomes

According to recent reports, many ribosomal RNA gene annotations are still questionable, and the use of inappropriate tools for annotation has been blamed. However, we believe that the abundant 16S rRNA partial sequence in the databases, mainly created by culture-independent PCR methods, is another main cause of the ambiguous annotations of 16S rRNA. To examine the current status of 16S rRNA gene annotations in complete microbial genomes, we used as a criterion the conserved anti-SD sequence, located at the 3′ end of the 16S rRNA gene, which is commonly overlooked by culture-independent PCR methods. In our large survey, 859 16S rRNA gene sequences from 252 different species of the microbial complete genomes were inspected. 67 species (234 genes) were detected with ambiguous annotations. The common anti-SD sequence and other conserved 16S rRNA sequence features could be detected in the downstream-intergenic regions for almost every questionable sequence, indicating that many of the 16S rRNA genes were annotated incorrectly. Furthermore, we found that more than 91.5% of the 93,716 sequences of the available 16S rRNA in the main databases are partial sequences. We also performed BLAST analysis for every questionable rRNA sequence, and most of the best hits in the analysis were rRNA partial sequences. This result indicates that partial sequences are prevalent in the databases, and that these sequences have significantly affected the accuracy of microbial genomic annotation. We suggest that the annotation of 16S rRNA genes in newly complete microbial genomes must be done in more detail, and that revision of questionable rRNA annotations should commence as soon as possible.     

Quantifying Microbial Diversity: Morphotypes, 16S rRNA Genes, and Carotenoids of Oxygenic Phototrophs in Microbial Mats

We quantified the diversity of oxygenic phototrophic microorganisms present in eight hypersaline microbial mats on the basis of three cultivation-independent approaches. Morphological diversity was studied by microscopy. The diversity of carotenoids was examined by extraction from mat samples and high-pressure liquid chromatography analysis. The diversity of 16S rRNA genes from oxygenic phototrophic microorganisms was investigated by extraction of total DNA from mat samples, amplification of 16S rRNA gene segments from cyanobacteria and plastids of eukaryotic algae by phylum-specific PCR, and sequence-dependent separation of amplification products by denaturing-gradient gel electrophoresis. A numerical approach was introduced to correct for crowding the results of chromatographic and electrophoretic analyses. Diversity estimates typically varied up to twofold among mats. The congruence of richness estimates and Shannon-Weaver indices based on numbers and proportional abundances of unique morphotypes, 16S rRNA genes, and carotenoids unveiled the underlying diversity of oxygenic phototrophic microorganisms in the eight mat communities studied.     

Species-specific detection of hydrocarbon-utilizing bacteria

Rapid detection and quantitative assessment of specific microbial species in environmental samples is desirable for monitoring changes in ecosystems and for tracking natural or introduced microbial species during bioremediation of contaminated sites. In the interests of developing rapid tests for hydrocarbon-degrading bacteria, species-specific PCR primer sets have been developed for Pseudomonas aeruginosa, Stentrophomonas (Xanthomonas) maltophilia, and Serratia marsescens. Highly variable regions of the 16S rRNA gene were used to design these primer sets. The amplification products of these primer sets have been verified and validated with hemi-nested PCR and with ligase chain reaction (LCR) techniques, and have been applied to the analyses of environmental water samples. These species-specific primer sets were also chosen to amplify in conjunction with a universal set of PCR primers chosen from highly conserved neighboring sequences in the same gene. These multiplex or competitive PCR procedures enable testing with an internal marker and/or the quantitative estimation of the relative proportion of the microbial community that any one of these species occupies. In addition, this universal PCR primer set amplified the same size amplicon from a wide spectrum of procaryotic and eucaryotic organisms and may have potential in earth biota analyses.     

Molecular identification of filterable bacteria and archaea in the water of acidic lakes of northern Russia

Wetland ecosystems are the natural centers of freshwater formation in northern Russia lowland landscapes. The humic acidic waters formed in bogs feed the numerous lakes of the northern regions. One milliliter of the water in these lakes contains up to 10⁴ ultrasmall microbial cells that pass through “bacterial” filters with a pore size of 0.22 μm. The vast majority of these cells do not grow on nutrient media and cannot be identified by routine cultivation-based approaches. Their identification was performed by analysis of clone libraries obtained by PCR amplification of archaeal and bacterial 16S rRNA genes from the fraction of cells collected from water filtrates of acidic lakes. Most of the obtained bacterial 16S rRNA gene sequences represented the class Betaproteobacteria and exhibited the highest homology of (94–99%) with 16S rRNA genes of representatives of the genera Herbaspirillum, Herminiimonas, Curvibacter, and Burkholderia. The archaeal 16S rRNA gene clone library comprised genes of Euryarchaeota representatives. One-third of these genes exhibited 97–99% homology to the 16S rRNA genes of taxonomically described organisms of the orders Methanobacteriales and Methanosarcinales. The rest of the cloned archaeal 16S rRNA genes were only distantly related (71–74% homology) to those in all earlier characterized archaea.     

Diversity and distribution in hypersaline microbial mats of bacteria related to Chloroflexus spp

Filamentous bacteria containing bacteriochlorophylls c and a were enriched from hypersaline microbial mats. Based on phylogenetic analyses of 16S rRNA gene sequences, these organisms form a previously undescribed lineage distantly related to Chloroflexus spp. We developed and tested a set of PCR primers for the specific amplification of 16S rRNA genes from filamentous phototrophic bacteria within the kingdom of "green nonsulfur bacteria." PCR products recovered from microbial mats in a saltern in Guerrero Negro, Mexico, were subjected to cloning or denaturing gradient gel electrophoresis and then sequenced. We found evidence of a high diversity of bacteria related to Chloroflexus which exhibit different distributions along a gradient of salinity from 5.5 to 16%.     

Selective phylogenetic analysis targeted at 16S rRNA genes of thermophiles and hyperthermophiles in deep-subsurface geothermal environments

Deep-subsurface samples obtained by deep drilling are likely to be contaminated with mesophilic microorganisms in the drilling fluid, and this could affect determination of the community structure of the geothermal microflora using 16S rRNA gene clone library analysis. To eliminate possible contamination by PCR-amplified 16S rRNA genes from mesophiles, a combined thermal denaturation and enzyme digestion method, based on a strong correlation between the G+C content of the 16S rRNA gene and the optimum growth temperatures of most known prokaryotic cultures, was used prior to clone library construction. To validate this technique, hot spring fluid (76 degrees C) and river water (14 degrees C) were used to mimic a deep-subsurface sample contaminated with drilling fluid. After DNA extraction and PCR amplification of the 16S rRNA genes from individual samples separately, the amplified products from river water were observed to be denatured at 82 degrees C and completely digested by exonuclease I (Exo I), while the amplified products from hot spring fluid remained intact after denaturation at 84 degrees C and enzyme digestion with Exo I. DNAs extracted from the two samples were mixed and used as a template for amplification of the 16S rRNA genes. The amplified rRNA genes were denatured at 84 degrees C and digested with Exo I before clone library construction. The results indicated that the 16S rRNA gene sequences from the river water were almost completely eliminated, whereas those from the hot spring fluid remained.     

Bacterial diversity based on 16S rRNA and gyrB genes at Yinshan mine, China

The diversity of bacterial communities at three sites impacted by acid mine drainage (AMD) from the Yinshan Mine in China was studied using comparative sequence analysis of two molecular markers, the 16S rRNA and gyrB genes. The phylogenetic analyses retrieved sequences from six classes of bacteria, Nitrospira, Alphaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Acidobacteria, and Actinobacteria, as well as sequences related to the plastid of the cyanobacterium Cyanidium acidocaldarium and also some unknown bacteria. The results of phylogenetic analyses based on gyrB and 16S rRNA were compared. This confirmed that gyrB gene analysis may be a useful tool, in addition to the comparative sequence analysis of the 16S rRNA gene, for the analysis of microbial community compositions. Moreover, the Mantel test showed that the geochemical characteristics, especially the pH value and the concentration of iron, strongly influenced the composition of the microbial communities.     

Nonradioactive method to study genetic profiles of natural bacterial communities by PCR-single-strand-conformation polymorphism.

We describe a new method for studying the structure and diversity of bacterial communities in the natural ecosystem. Our approach is based on single-strand-conformation polymorphism (SSCP) analysis of PCR products of 16S rRNA genes from complex bacterial populations. A pair of eubacterial universal primers for amplification of the variable V3 region were designed from the 16S rRNA sequences of 1,262 bacterial strains. The PCR conditions were optimized by using genomic DNAs from five gram-positive and seven gram-negative strains. The SSCP analysis of the PCR products demonstrated that a bacterial strain generated its characteristic band pattern and that other strains generated other band patterns, so that the relative diversity in bacterial communities could be measured. In addition, this method was sensitive enough to detect a bacterial population that made up less than 1.5% of a bacterial community. The distinctive differences between bacterial populations were observed in an oligotrophic lake and a eutrophic pond in a field study. The method presented here, using combined PCR amplification and SSCP pattern analyses of 16S rRNA genes, provides a useful tool to study bacterial community structures in various ecosystems.