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

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

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.     

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.     

Quantification of mcrA by fluorescent PCR in methanogenic and methanotrophic microbial communities

A quantitative fluorogenic PCR method for detecting methanogenic and methanotrophic orders was established using a refined primer set for the methyl coenzyme M reductase subunit A gene (mcrA). The method developed was applied to several microbial communities in which diversity and abundance of methanogens or anaerobic methanotrophs (ANMEs) was identified by 16S rRNA gene clone analysis, and strong correlations between the copy numbers of mcrA with those of archaeal 16S rRNA genes in the communities were observed. The assay can be applied to detecting and assessing the abundance of methanogens and/or ANMEs in anoxic environments that could not be detected by 16S rRNA gene sequence analyses.     

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.     

引物应用策略：基于2对古菌16S rRNA基因通用引物对环境样品古菌多样性分析

【目的】找到适宜的16S rRNA基因通用引物应用策略,应对复杂环境微生物多样性调查,尤其目前高速发展的高通量测序技术带来的巨大挑战。【方法】用Oligocheck软件分别将两对应试的古菌16S rRNA基因通用引物与RDP(Ribosomal database project)数据库中古菌16S rRNA基因序列进行匹配比对。用两对应试引物分别构建海洋沉积物样品的古菌16S rRNA基因文库。【结果】软件匹配结果显示引物f109/r958与目的基因的匹配程度高于引物f21/r958。该结果与古菌16S rRNA基因文库RFLP分析、古菌多样性指数分析结果相吻合。数据还表明,2对引物的综合文库能更好满足该沉积物样品的古菌多样性分析。【结论】选用与数据库中目的基因匹配性高的通用引物和多个引物的联合使用,可以有效提高环境样品微生物多样性调查的分辨率。     

An interlaboratory comparison of 16S rRNA gene-based terminal restriction fragment length polymorphism and sequencing methods for assessing microbial diversity of seafloor basalts

We present an interlaboratory comparison between full-length 16S rRNA gene sequence analysis and terminal restriction fragment length polymorphism (TRFLP) for microbial communities hosted on seafloor basaltic lavas, with the goal of evaluating how similarly these two different DNA-based methods used in two independent labs would estimate the microbial diversity of the same basalt samples. Two samples were selected for these analyses based on differences detected in the overall levels of microbial diversity between them. Richness estimators indicate that TRFLP analysis significantly underestimates the richness of the relatively high-diversity seafloor basalt microbial community: at least 50% of species from the high-diversity site are missed by TRFLP. However, both methods reveal similar dominant species from the samples, and they predict similar levels of relative diversity between the two samples. Importantly, these results suggest that DNA-extraction or PCR-related bias between the two laboratories is minimal. We conclude that TRFLP may be useful for relative comparisons of diversity between basalt samples, for identifying dominant species, and for estimating the richness and evenness of low-diversity, skewed populations of seafloor basalt microbial communities, but that TRFLP may miss a majority of species in relatively highly diverse samples.     

Genetic diversity of bacterial strains isolated from soils, contaminated with polycyclic aromatic hydrocarbons, by 16S rRNA gene sequencing and amplified fragment length polymorphism fingerprinting

In order to study microbial diversity in a polycyclic aromatic hydrocarbon-impacted soil, 14 bacterial strains were analyzed by 16S rRNA gene sequencing and amplified fragment length polymorphism (AFLP) analysis. Bacterial strains isolated from two different hydrocarbon-polluted sites were identified to the species level by 16S rRNA full-gene sequencing using MicroSeq 16S rRNA gene sequencing. Their genome was subsequently analyzed by high-resolution genotyping with AFLP analysis, in order to monitor species variability and to differentiate closely related strains. Cluster analysis based on AFLP fingerprinting showed intra-specific polymorphism, even among strains with 100% 16S rRNA gene sequence identity. The results show that AFLP is a powerful, highly reproducible and discriminatory tool for revealing genetic relationships in bacterial populations. The ability to differentiate and track related closely microbes is fundamental for studying structure and dynamics of microbial communities in contaminated ecosystems.     

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.     

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.     

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.     

Individual hydrothermal vents at Axial Seamount harbor distinct subseafloor microbial communities

The microbial community structure of five geographically distinct hydrothermal vents located within the Axial Seamount caldera, Juan de Fuca Ridge, was examined over 6 years following the 1998 diking eruptive event. Terminal restriction fragment length polymorphism (TRFLP) and 16S rRNA gene sequence analyses were used to determine the bacterial and archaeal diversity, and the statistical software primer v6 was used to compare vent microbiology, temperature and fluid chemistry. Statistical analysis of vent fluid temperature and composition shows that there are significant differences between vents in any year, but that the fluid composition changes over time such that no vent maintains a chemical composition completely distinct from the others. In contrast, the subseafloor microbial communities associated with individual vents changed from year to year, but each location maintained a distinct community structure (based on TRFLP and 16S rRNA gene sequence analyses) that was significantly different from all other vents included in this study. Epsilonproteobacterial microdiversity is shown to be important in distinguishing vent communities, while archaeal microdiversity is less variable between sites. We propose that persistent venting at diffuse flow vents over time creates the potential to isolate and stabilize diverse microbial community structures between vents.     

基于16S rRNA基因测序分析微生物群落多样性

微生物群落多样性的研究对于挖掘微生物资源,探索微生物群落功能,阐明微生物群落与生境间的关系具有重要意义。随着宏基因组概念的提出以及测序技术的快速发展,16S rRNA基因测序在微生物群落多样性的研究中已被广泛应用。文中系统地介绍了16S rRNA基因测序分析流程中的四个重要环节,包括测序平台与扩增区的选择、测序数据预处理以及多样性分析方法,就其面临的问题与挑战进行了探讨并对未来的研究方向进行了展望,以期为微生物群落多样性相关研究提供参考。     

Influence of depth and time on diversity of free-living microbial community in the variably saturated zone of a granitic aquifer

We investigated the temporal and vertical changes in the microbial communities related to hydrological variations an aquifer (Brittany, France). Five water samplings were carried out, spanning three hydrological cycles in the variably and the permanently saturated zones. Seasonal variations in the major anion concentrations (NO3 -, SO4 2- and Cl(-) ) indicated that different physical processes occurred during the recharge process in the two zones. The variably saturated zone is mainly dominated by diffusion and advection processes from the soil, whereas the permanently saturated zone is controlled by moderate advective transfer from the variably saturated zone. Bacterial diversity was investigated by flow cytometry, 16S rRNA and narG genes analyses. Part of this diversity was new in that 6 of the 27 16S rRNA gene sequence phylotypes were unknown even at the class or phylum level. The narG gene analysis did not reveal any clear variation in time or depth within the nitrate reducers' community. In contrast, 16S rRNA gene analyses showed modifications of community composition that could be related to the hydrologic and chemical contrast between the two zones. It was concluded that the physical processes of water transfer could influence bacterial diversity at the soil-aquifer interface.     

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.     

八门湾红树林土壤芽胞杆菌分离与多样性分析

【目的】了解八门湾红树林海漆林区土壤中可培养芽胞杆菌资源的多样性。【方法】采用水浴处理与直接涂布相结合的方法选择性分离土壤中的芽胞杆菌;利用16S rDNA PCR-RFLP与16S rDNA序列分析技术研究可培养芽胞杆菌资源的遗传多样性和系统发育关系。【结果】16S rDNA PCR-RFLP酶切图谱UPGMA聚类分析表明,在100%的相似性水平上,分离的155株芽胞杆菌分属21个遗传类群,显示了较为丰富的遗传多样性;由21种遗传类型代表菌株的16S rDNA序列分析结果得知,这些芽胞杆菌主要分布在Bacillaceae和Paenibacillaceae科下的Bacillus、Halobacillus、Virgibacillus和Paenibacillus 4个属,其中Bacillus为优势属;有8株芽胞杆菌的16S rDNA序列与数据库中相应模式菌株的最大相似性在95.1%-99.0%之间。【结论】八门湾红树林土壤可培养芽胞杆菌有着较为丰富的遗传多样性,并存在新的芽胞杆菌物种资源。     

分子生物学方法在微生物多样性研究中的应用

微生物多样性是生物多样性的重要组成部分。由于微生物和大生物（动、植物）相比,存在着多种显著差异,因此其多样性,保护及利用也有所不同,尤其是研究方法亟待完善,提高。近年来,分子生物学方法广泛用于微生物多样性的研究并取得了一系列研究成果。本文从四个方面加以介绍：１）微生物总ＤＮＡ制备及其遗传多样性检测方法;２）１６ＳｒＲＮＡ基因序列研究;３）核酸杂交分析技术;４）ＤＮＡ动力学的研究。今后的发展趋势是加     

Combined use of cultivation-dependent and cultivation-independent methods indicates that members of most haloarchaeal groups in an Australian crystallizer pond are cultivable

Haloarchaea are the dominant microbial flora in hypersaline waters with near-saturating salt levels. The haloarchaeal diversity of an Australian saltern crystallizer pond was examined by use of a library of PCR-amplified 16S rRNA genes and by cultivation. High viable counts (10(6) CFU/ml) were obtained on solid media. Long incubation times (> or =8 weeks) appeared to be more important than the medium composition for maximizing viable counts and diversity. Of 66 isolates examined, all belonged to the family Halobacteriaceae, including members related to species of the genera Haloferax, Halorubrum, and Natronomonas. In addition, isolates belonging to a novel group (the ADL group), previously detected only as 16S rRNA genes in an Antarctic hypersaline lake (Deep Lake), were cultivated for the first time. The 16S rRNA gene library identified the following five main groups: Halorubrum groups 1 and 2 (49%), the SHOW (square haloarchaea of Walsby) group (33%), the ADL group (16%), and the Natronomonas group (2%). There were two significant differences between the organisms detected in cultivation and 16S rRNA sequence results. Firstly, Haloferax spp. were frequently isolated on plates (15% of all isolates) but were not detected in the 16S rRNA sequences. Control experiments indicated that a bias against Haloferax sequences in the generation of the 16S rRNA gene library was unlikely, suggesting that Haloferax spp. readily form colonies, even though they were not a dominant group. Secondly, while the 16S rRNA gene library identified the SHOW group as a major component of the microbial community, no isolates of this group were obtained. This inability to culture members of the SHOW group remains an outstanding problem in studying the ecology of hypersaline environments.     

Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such "heterogeneity hot spots" occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.