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.     

Comparison of 16S rRNA and protein-coding genes as molecular markers for assessing microbial diversity (Bacteria and Archaea) in ecosystems

PCR amplification of the rRNA gene is the most popular method for assessing microbial diversity. However, this molecular marker is often present in multiple copies in cells presenting, in addition, an intragenomic heterogeneity. In this context, housekeeping genes may be used as taxonomic markers for ecological studies. However, the efficiency of these protein-coding genes compared to 16S rRNA genes has not been tested on environmental data. For this purpose, five protein marker genes for which primer sets are available, were selected (rplB, pyrG, fusA, leuS and rpoB) and compared with 16S rRNA gene results from PCR amplification or metagenomic data from aquatic ecosystems. Analysis of the major groups found in these ecosystems, such as Actinobacteria, Bacteroides, Proteobacteria and Cyanobacteria, showed good agreement between the protein markers and the results given by 16S rRNA genes from metagenomic reads. However, with the markers it was possible to detect minor groups among the microbial assemblages, providing more details compared to 16S rRNA results from PCR amplification. In addition, the use of a set of protein markers made it possible to deduce a mean copy number of rRNA operons. This average estimate is essentially lower than the one estimated in sequenced genomes.     

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

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

Intragenomic Heterogeneity of 16S rRNA Genes Causes Overestimation of Prokaryotic Diversity

Ever since Carl Woese introduced the use of 16S rRNA genes for determining the phylogenetic relationships of prokaryotes, this method has been regarded as the “gold standard” in both microbial phylogeny and ecology studies. However, intragenomic heterogeneity within 16S rRNA genes has been reported in many investigations and is believed to bias the estimation of prokaryotic diversity. In the current study, 2,013 completely sequenced genomes of bacteria and archaea were analyzed and intragenomic heterogeneity was found in 952 genomes (585 species), with 87.5% of the divergence detected being below the 1% level. In particular, some extremophiles (thermophiles and halophiles) were found to harbor highly divergent 16S rRNA genes. Overestimation caused by 16S rRNA gene intragenomic heterogeneity was evaluated at different levels using the full-length and partial 16S rRNA genes usually chosen as targets for pyrosequencing. The result indicates that, at the unique level, full-length 16S rRNA genes can produce an overestimation of as much as 123.7%, while at the 3% level, an overestimation of 12.9% for the V6 region may be introduced. Further analysis showed that intragenomic heterogeneity tends to concentrate in specific positions, with the V1 and V6 regions suffering the most intragenomic heterogeneity and the V4 and V5 regions suffering the least intragenomic heterogeneity in bacteria. This is the most up-to-date overview of the diversity of 16S rRNA genes within prokaryotic genomes. It not only provides general guidance on how much overestimation can be introduced when applying 16S rRNA gene-based methods, due to its intragenomic heterogeneity, but also recommends that, for bacteria, this overestimation be minimized using primers targeting the V4 and V5 regions.     

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.     

16S rRNA GENE HETEROGENEITY IN THE FILAMENTOUS MARINE CYANOBACTERIAL GENUS LYNGBYA1

The SSU (16S) rRNA gene was used to investigate the phylogeny of the cyanobacterial genus Lyngbya as well as examined for its capacity to discriminate between different marine species of Lyngbya. We show that Lyngbya forms a polyphyletic genus composed of a marine lineage and a halophilic/brackish/freshwater lineage. In addition, we found morphological and genetic evidence that Lyngbya spp. often grow in association with other microorganisms, in particular smaller filamentous cyanobacteria such as Oscillatoria, and propose that these associated microorganisms have led to extensive phylogenetic confusion in identification of Lyngbya spp. At the species level, the phylogenetic diversity obtained from the comparison of 16S rRNA genes exceeded morphological diversity in Lyngbya. However, the expectation that this improved phylogeny would be useful to species and subspecies identification was eliminated by the fact that phylogenetic species did not correlate in any respect with the species obtained from current taxonomic systems. In addition, phylogenetic identification was adversely affected by the presence of multiple gene copies within individual Lyngbya colonies. Analysis of clonal Lyngbya cultures and multiple displacement amplified (MDA) single‐cell genomes revealed that Lyngbya genomes contain two 16S rRNA gene copies, and that these typically are of variable sequence. Furthermore, intragenomic and interspecies 16S rRNA gene heterogeneity was approximately of the same magnitude. Hence, the intragenomic heterogeneity of the 16S rRNA gene overestimates the microdiversity of different strains and does not accurately reflect speciation within cyanobacteria, including the genus Lyngbya.     

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-based microbial community analysis avoids limitations inherent in 16S rRNA gene intraspecies heterogeneity

Contemporary microbial community analysis frequently involves PCR-amplified sequences of the 16S rRNA gene (rDNA). However, this technology carries the inherent problem of heterogeneity between copies of the 16S rDNA in many species. As an alternative to 16S rDNA sequences in community analysis, we employed the gene for the RNA polymerase beta subunit (rpoB), which appears to exist in one copy only in bacteria. In the present study, the frequency of 16S rDNA heterogeneity in bacteria isolated from the marine environment was assessed using bacterial isolates from the red alga Delisea pulchra and from the surface of a marine rock. Ten strains commonly used in our laboratory were also assessed for the degree of heterogeneity between the copies of 16S rDNA and were used to illustrate the effect of this heterogeneity on microbial community pattern analysis. The rock isolates and the laboratory strains were also used to confirm nonheterogeneity of rpoB, as well as to investigate the versatility of the primers. In addition, a comparison between 16S rDNA and rpoB PCR-DGGE (denaturing gradient gel electrophoresis)-based community analyses was performed using a DNA mixture of nine isolates from D. pulchra. Eight out of 14 isolates from D. pulchra, all rock isolates, and 6 of 10 laboratory strains displayed multiple bands for 16S rDNA when analyzed by DGGE. There was no indication of heterogeneity for either the rock isolates or the laboratory strains when rpoB was used for PCR-DGGE analysis. Microbial community pattern analysis using 16S rDNA PCR-DGGE showed an overestimation of the number of laboratory strains in the sample, while some strains were not represented. Therefore, the 16S rDNA PCR-DGGE-based community analysis was proven to be severely limited by 16S rDNA heterogeneity. The mixture of isolates from D. pulchra proved to be more accurately described using rpoB, compared to the 16S rDNA-based PCR-DGGE.     

Multiple copies of 16s rRNA gene affect the restriction patterns and DGGE profile as revealed by analysis of genome database

The use of 16S rRNA gene has been a "golden" method to determine the diversity of microbial communities in environmental samples, phylogenetic relationships of prokaryotes and taxonomic position of newly isolated organisms. However due to the presence of multiple heterogeneous 16S rRNA gene copies in many strains, the interpretation of microbial ecology via 16S rRNA sequences is complicated. Purpose of present paper is to demonstrate the extent to which the multiple heterogeneous 16S rRNA gene copies affect RFLP patterns and DGG E profiles by using the genome database. In present genome database, there are 782 bacterial strains in total whose genomes have been completely sequenced and annotated. Among the total strains, 639 strains (82%) possess multiple 16S rRNA gene copies, 415 strains (53%) whose multiple copies are heterogeneous in sequences as revealed by alignment, 236 strains (30%) whose multiple copies show different restrict patterns by CSP61 + Hinfl, MspI + Rsal or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G + C% and phylogentic distances based on the sequences of V3 region, which are linked to DGGE patters. Polymorphisms of a few strains were shown as examples. Using artificial communities, it is demonstrated that the presence of multiple heterogeneous 16S rRNA gene copies potentially leads to over-estimation of the diversity of a community. It is suggested that care must be taken when interpreting 16S rRNA-based RFLP and DGGE data and profiling an environmental community.     

A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity

Background

The 16S rRNA gene is the gold standard in molecular surveys of bacterial and archaeal diversity, but it has the disadvantages that it is often multiple-copy, has little resolution below the species level and cannot be readily interpreted in an evolutionary framework. We compared the 16S rRNA marker with the single-copy, protein-coding rpoB marker by amplifying and sequencing both from a single soil sample. Because the higher genetic resolution of the rpoB gene prohibits its use as a universal marker, we employed consensus-degenerate primers targeting the Proteobacteria.

Methodology/Principal Findings

Pyrosequencing can be problematic because of the poor resolution of homopolymer runs. As these erroneous runs disrupt the reading frame of protein-coding sequences, removal of sequences containing nonsense mutations was found to be a valuable filter in addition to flowgram-based denoising. Although both markers gave similar estimates of total diversity, the rpoB marker revealed more species, requiring an order of magnitude fewer reads to obtain 90% of the true diversity. The application of population genetic methods was demonstrated on a particularly abundant sequence cluster.

Conclusions/Significance

The rpoB marker can be a complement to the 16S rRNA marker for high throughput microbial diversity studies focusing on specific taxonomic groups. Additional error filtering is possible and tests for recombination or selection can be employed.     

Multiple copies of 16S rRNA gene affect the restriction patterns and DGGE profile revealed by analysis of genome database

The use of 16S rRNA gene has been a “golden” method to determine the diversity of microbial communities in environmental samples, phylogenetic relationships of prokaryotes and taxonomic position of newly isolated organisms. However due to the presence of multiple heterogeneous 16S rRNA gene copies in many strains, the interpretation of microbial ecology via 16S rRNA sequences is complicated. Purpose of present paper is to demonstrate the extent to which the multiple heterogeneous 16S rRNA gene copies affect RFLP patterns and DGGE profiles by using the genome database. In present genome database, there are 782 bacterial strains in total whose genomes have been completely sequenced and annotated. Among the total strains, 639 strains (82%) possess multiple 16S rRNA gene copies, 415 strains (53%) whose multiple copies are heterogeneous in sequences as revealed by alignment, 236 strains (30%) whose multiple copies show different restrict patterns by CSP6I+HinfI, MspI+RsaI or HhaI as analyzed in silico. Polymorphisms of the multiple copies in certain strains were further characterized by G+C% and phy-logentic distances based on the sequences of V3 region, which are linked to DGGE patters. Polymorphisms of a few strains were shown as examples. Using artificial communities, it is demonstrated that the presence of multiple heterogeneous 16S rRNA gene copies potentially leads to over-estimation of the diversity of a community. It is suggested that care must be taken when interpreting 16S rRNA-based RFLP and DGGE data and profiling an environmental community.     

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.     

台湾海峡鱼类绦虫的分子系统学研究

利用DNA测序技术对台湾海峡部分鱼类绦虫的16S rRNA和18S rRNA基因片段序列进行了分析。使用PAUP4·0b10软件构建的进化树显示,目前关于绦虫二叶目、锥吻目、假叶目、盘头目和四叶目的划分是比较合理的,绦虫进化基本遵循了头节形态从简单到复杂的进化规律。报道了国内首次发现的双叶目绦虫,进化树结果初步支持了巨槽属和棘头属的划分。此外,结果也支持了前孔属绦虫的分类地位。但是,对耳槽属绦虫与阶室属绦虫的形态学划分与分子系统学相矛盾,利用16S rRNA基因对盘头目各种的进化树分析与形态学差异很大,这些问题都需要更多研究来进行深入分析。     

Ribosomal Multi-Operon Diversity: An Original Perspective on the Genus Aeromonas

16S rRNA gene (rrs) is considered of low taxonomic interest in the genus Aeromonas. Here, 195 Aeromonas strains belonging to populations structured by multilocus phylogeny were studied using an original approach that considered Ribosomal Multi-Operon Diversity. This approach associated pulsed-field gel electrophoresis (PFGE) to assess rrn operon number and distribution across the chromosome and PCR-temporal temperature gel electrophoresis (TTGE) to assess rrs V3 region heterogeneity. Aeromonads harbored 8 to 11 rrn operons, 10 operons being observed in more than 92% of the strains. Intraspecific variability was low or nul except for A. salmonicida and A. aquariorum suggesting that large chromosomic rearrangements might occur in these two species while being extremely rarely encountered in the evolution of other taxa. rrn operon number at 8 as well as PFGE patterns were shown valuable for taxonomic purpose allowing resolution of species complexes. PCR-TTGE revealed a high rate of strains (41.5%) displaying intragenomic rrs heterogeneity. Strains isolated from human samples more frequently displayed intragenomic heterogeneity than strains recovered from non-human and environmental specimens. Intraspecific variability ranged from 0 to 76.5% of the strains. The observation of species-specific TTGE bands, the recovery of identical V3 regions in different species and the variability of intragenomic heterogeneity (1–13 divergent nucleotides) supported the occurrence of mutations and horizontal transfer in aeromonad rrs evolution. Altogether, the presence of a high number of rrn operon, the high proportion of strains harboring divergent rrs V3 region and the previously demonstrated high level of genetic diversity argued in favor of highly adaptative capabilities of aeromonads. Outstanding features observed for A. caviae supported the ongoing process of adaptation to a specialized niche represented by the gut, previously hypothesized. 16S rRNA gene is an informative marker in the genus Aeromonas for both evolutionary and polyphasic taxonomic studies provided that multi-operon fingerprinting approaches are used.     

Taxonomic Resolutions Based on 18S rRNA Genes: A Case Study of Subclass Copepoda

Biodiversity studies are commonly conducted using 18S rRNA genes. In this study, we compared the inter-species divergence of variable regions (V1–9) within the copepod 18S rRNA gene, and tested their taxonomic resolutions at different taxonomic levels. Our results indicate that the 18S rRNA gene is a good molecular marker for the study of copepod biodiversity, and our conclusions are as follows: 1) 18S rRNA genes are highly conserved intra-species (intra-species similarities are close to 100%); and could aid in species-level analyses, but with some limitations; 2) nearly-whole-length sequences and some partial regions (around V2, V4, and V9) of the 18S rRNA gene can be used to discriminate between samples at both the family and order levels (with a success rate of about 80%); 3) compared with other regions, V9 has a higher resolution at the genus level (with an identification success rate of about 80%); and 4) V7 is most divergent in length, and would be a good candidate marker for the phylogenetic study of Acartia species. This study also evaluated the correlation between similarity thresholds and the accuracy of using nuclear 18S rRNA genes for the classification of organisms in the subclass Copepoda. We suggest that sample identification accuracy should be considered when a molecular sequence divergence threshold is used for taxonomic identification, and that the lowest similarity threshold should be determined based on a pre-designated level of acceptable accuracy.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene ( rpoB  ) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

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.     

Phylogenetic relationships among Staphylococcus species and refinement of cluster groups based on multilocus data

ABSTRACT: BACKGROUND: Estimates of relationships among Staphylococcus species have been hampered by poor and inconsistent resolution of phylogenies based largely on single gene analyses incorporating only a limited taxon sample. As such, the evolutionary relationships and hierarchical classification schemes among species have not been confidently established. Here, we address these points through analyses of DNA sequence data from multiple loci (16S rRNA gene, dnaJ, rpoB, and tuf gene fragments) using multiple Bayesian and maximum likelihood phylogenetic approaches that incorporate nearly all recognized Staphylococcus taxa. RESULTS: We estimated the phylogeny of fifty-seven Staphylococcus taxa using partitioned-model Bayesian and maximum likelihood analysis, as well as Bayesian gene-tree species-tree methods. Regardless of methodology, we found broad agreement among methods that the current cluster groups require revision, although there was some disagreement among methods in resolution of higher order relationships. Based on our phylogenetic estimates, we propose a refined classification for Staphylococcus with species being classified into 15 cluster groups (based on molecular data) that adhere to six species groups (based on phenotypic properties) CONCLUSIONS: Our findings are in general agreement with gene tree-based reports of the staphylococcal phylogeny, although we identify multiple previously unreported relationships among species. Our results support the general importance of such multilocus assessments as a standard in microbial studies to more robustly infer relationships among recognized and newly discovered.     

Use of rpoB gene analysis for identification of nitrogen-fixing Paenibacillus species as an alternative to the 16S rRNA gene

AIM: To avoid the limitations of 16S rRNA-based phylogenetic analysis for Paenibacillus species, the usefulness of the RNA polymerase beta-subunit encoding gene (rpoB) was investigated as an alternative to the 16S rRNA gene for taxonomic studies. METHODS AND RESULTS: Partial rpoB sequences were generated for the type strains of eight nitrogen-fixing Paenibacillus species. The presence of only one copy of rpoB in the genome of P. graminis strain RSA19(T) was demonstrated by denaturing gradient gel electrophoresis and hybridization assays. A comparative analysis of the sequences of the 16S rRNA and rpoB genes was performed and the eight species showed between 91.6-99.1% (16S rRNA) and 77.9-97.3% (rpoB) similarity, allowing a more accurate discrimination between the different species using the rpoB gene. Finally, 24 isolates from the rhizosphere of different cultivars of maize previously identified as Paenibacillus spp. were assigned correctly to one of the nitrogen-fixing species. CONCLUSIONS, SIGNIFICANCE AND IMPACT OF THE STUDY: The data obtained in this study indicate that rpoB is a powerful identification tool, which can be used for the correct discrimination of the nitrogen-fixing species of agricultural and industrial importance within the genus Paenibacillus.     

Confident phylogenetic identification of uncultured prokaryotes through long read amplicon sequencing of the 16S-ITS-23S rRNA operon

Amplicon sequencing of the 16S rRNA gene is the predominant method to quantify microbial compositions and to discover novel lineages. However, traditional short amplicons often do not contain enough information to confidently resolve their phylogeny. Here we present a cost-effective protocol that amplifies a large part of the rRNA operon and sequences the amplicons with PacBio technology. We tested our method on a mock community and developed a read-curation pipeline that reduces the overall read error rate to 0.18%. Applying our method on four environmental samples, we captured near full-length rRNA operon amplicons from a large diversity of prokaryotes. The method operated at moderately high-throughput (22286–37,850 raw ccs reads) and generated a large amount of putative novel archaeal 23S rRNA gene sequences compared to the archaeal SILVA database. These long amplicons allowed for higher resolution during taxonomic classification by means of long (∼1000 bp) 16S rRNA gene fragments and for substantially more confident phylogenies by means of combined near full-length 16S and 23S rRNA gene sequences, compared to shorter traditional amplicons (250 bp of the 16S rRNA gene). We recommend our method to those who wish to cost-effectively and confidently estimate the phylogenetic diversity of prokaryotes in environmental samples at high throughput.