Novel Major Bacterial Candidate Division within a Municipal Anaerobic Sludge Digester

In a previous study, we analyzed the molecular diversity of Planctomycetales by PCR amplification and sequencing of 16S rRNA clone libraries generated from a municipal wastewater plant, using planctomycete-specific and universal primer sets (R. Chouari, D. Le Paslier, P. Daegelen, P. Ginestet, J. Weissenbach, and A. Sghir, Appl. Environ. Microbiol. 69:7354-7363, 2003). Only a small fraction (4%) of the 16S rRNA gene sequences of the digester clone library corresponded to the Planctomycetales division. Importantly, 85.9% of the digester clone sequences are grouped into two different clusters named WWE1 (81.4% of the sequences) and WWE2 (4.5%) and are distantly affiliated with unidentified bacterial sequences retrieved from a methanogenic reactor community and from a termite gut, respectively. In phylogenetic analysis using 16S rRNA gene sequence representatives of the main phylogenetic bacterial divisions, the two clusters are monophyletic, branch apart from each other, and are distantly related to Planctomycetales and other bacterial divisions. A novel candidate division is proposed for WWE1, while the WWE2 cluster strongly affiliates with the recently proposed Lentisphearae phylum. We designed and validated a 16S rRNA probe targeting WWE1 16S rRNA sequences by both fluorescent in situ hybridization (FISH) and dot blot hybridization (DBH). Results of FISH analysis show that WWE1 representative microorganisms are rods or filamentous shaped, while DBH shows that WWE1 accounts for 12% of the total bacterial rRNA within the anaerobic digester. The remaining 16S rRNA gene sequences are affiliated with Verrucomicrobia or recently described candidate divisions with no known pure culture representatives, such as OD1, BRC1, or NBL-UPA2, making up less than 3.5% of the clone library, respectively. This inventory expands the known diversity of the latter bacterial division-level lineages.     

Novel major bacterial candidate division within a municipal anaerobic sludge digester

In a previous study, we analyzed the molecular diversity of Planctomycetales by PCR amplification and sequencing of 16S rRNA clone libraries generated from a municipal wastewater plant, using planctomycete-specific and universal primer sets (R. Chouari, D. Le Paslier, P. Daegelen, P. Ginestet, J. Weissenbach, and A. Sghir, Appl. Environ. Microbiol. 69:7354-7363, 2003). Only a small fraction (4%) of the 16S rRNA gene sequences of the digester clone library corresponded to the Planctomycetales division. Importantly, 85.9% of the digester clone sequences are grouped into two different clusters named WWE1 (81.4% of the sequences) and WWE2 (4.5%) and are distantly affiliated with unidentified bacterial sequences retrieved from a methanogenic reactor community and from a termite gut, respectively. In phylogenetic analysis using 16S rRNA gene sequence representatives of the main phylogenetic bacterial divisions, the two clusters are monophyletic, branch apart from each other, and are distantly related to Planctomycetales and other bacterial divisions. A novel candidate division is proposed for WWE1, while the WWE2 cluster strongly affiliates with the recently proposed Lentisphearae phylum. We designed and validated a 16S rRNA probe targeting WWE1 16S rRNA sequences by both fluorescent in situ hybridization (FISH) and dot blot hybridization (DBH). Results of FISH analysis show that WWE1 representative microorganisms are rods or filamentous shaped, while DBH shows that WWE1 accounts for 12% of the total bacterial rRNA within the anaerobic digester. The remaining 16S rRNA gene sequences are affiliated with Verrucomicrobia or recently described candidate divisions with no known pure culture representatives, such as OD1, BRC1, or NBL-UPA2, making up less than 3.5% of the clone library, respectively. This inventory expands the known diversity of the latter bacterial division-level lineages.     

Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons

The diversity of the equine fecal bacterial community was evaluated using pyrosequencing of 16S rRNA gene amplicons. Fecal samples were obtained from horses fed cool-season grass hay. Fecal bacteria were characterized by amplifying the V4 region of bacterial 16S rRNA gene. Of 5898 mean unique sequences, a mean of 1510 operational taxonomic units were identified in the four fecal samples. Equine fecal bacterial richness was higher than that reported in humans, but lower than that reported in either cattle feces or soil. Bacterial classified sequences were assigned to 16 phyla, of which 10 were present in all samples. The largest number of reads belonged to Firmicutes (43.7% of total bacterial sequences), Verrucomicrobia (4.1%), Proteobacteria (3.8%), and Bacteroidetes (3.7%). The less abundant Actinobacteria, Cyanobacteria, and TM7 phyla presented here have not been previously described in the gut contents or feces of horses. Unclassified sequences represented 38.1% of total bacterial sequences; therefore, the equine fecal microbiome diversity is likely greater than that described. This is the first study to characterize the fecal bacterial community in horses by the use of 16S rRNA gene amplicon pyrosequencing, expanding our knowledge of the fecal microbiota of forage-fed horses.     

Use of real-time PCR for group-specific quantification of aceticlastic methanogens in anaerobic processes: population dynamics and community structures

The TaqMan quantitative PCR (QPCR) method was used to detect and quantify the 16S rRNA genes of aceticlastic methanogens at different taxonomic levels. Three different sets of primers coupled with a TaqMan probe for QPCR assays to detect the 16S rRNA genes of the order Methanosarcinales, as well as the families Methanosarcinaceae and Methanosaetaceae, were separately used. Using these primer and probe sets, the 16S rRNA genes of aceticlastic methanogens in samples from various anaerobic processes (i.e., nine pure cultures, batch experiment, and three different continuous processes including a full-scale digester), were monitored and quantified by QPCR assays. A batch experiment cultivating a mixture of aceticlastic methanogens, was conducted to monitor their population dynamics. Using this group-specific quantification method, the dynamics of a competition between two aceticlastic populations, as modulated by the acetate concentration, could well be described. The target 16S rRNA genes in environmental samples, collected from three different anaerobic processes treating sludge, cheese whey, and synthetic wastewaters, were additionally quantified. The quantified 16S rRNA gene concentrations for all samples successfully represented the community structures of the target methanogens, which were correlated accurately with the operational parameters of the anaerobic processes. It was also successful to demonstrate probe nesting of aceticlastic methanogens at the levels of order and family.     

Selection of primers for optimal taxonomic classification of environmental 16S rRNA gene sequences

Microbial community profiling using 16S rRNA gene sequences requires accurate taxonomy assignments. ‘Universal'' primers target conserved sequences and amplify sequences from many taxa, but they provide variable coverage of different environments, and regions of the rRNA gene differ in taxonomic informativeness—especially when high-throughput short-read sequencing technologies (for example, 454 and Illumina) are used. We introduce a new evaluation procedure that provides an improved measure of expected taxonomic precision when classifying environmental sequence reads from a given primer. Applying this measure to thousands of combinations of primers and read lengths, simulating single-ended and paired-end sequencing, reveals that these choices greatly affect taxonomic informativeness. The most informative sequence region may differ by environment, partly due to variable coverage of different environments in reference databases. Using our Rtax method of classifying paired-end reads, we found that paired-end sequencing provides substantial benefit in some environments including human gut, but not in others. Optimal primer choice for short reads totaling 96 nt provides 82–100% of the confident genus classifications available from longer reads.     

Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs

Massively parallel high throughput sequencing technologies allow us to interrogate the microbial composition of biological samples at unprecedented resolution. The typical approach is to perform high-throughout sequencing of 16S rRNA genes, which are then taxonomically classified based on similarity to known sequences in existing databases. Current technologies cause a predicament though, because although they enable deep coverage of samples, they are limited in the length of sequence they can produce. As a result, high-throughout studies of microbial communities often do not sequence the entire 16S rRNA gene. The challenge is to obtain reliable representation of bacterial communities through taxonomic classification of short 16S rRNA gene sequences. In this study we explored properties of different study designs and developed specific recommendations for effective use of short-read sequencing technologies for the purpose of interrogating bacterial communities, with a focus on classification using naïve Bayesian classifiers. To assess precision and coverage of each design, we used a collection of ∼8,500 manually curated 16S rRNA gene sequences from cultured bacteria and a set of over one million bacterial 16S rRNA gene sequences retrieved from environmental samples, respectively. We also tested different configurations of taxonomic classification approaches using short read sequencing data, and provide recommendations for optimal choice of the relevant parameters. We conclude that with a judicious selection of the sequenced region and the corresponding choice of a suitable training set for taxonomic classification, it is possible to explore bacterial communities at great depth using current technologies, with only a minimal loss of taxonomic resolution.     

Bacterial community in the crop of the hoatzin, a neotropical folivorous flying bird

The hoatzin is unique among known avian species because of the fermentative function of its enlarged crop. A small-bodied flying foregut fermenter is a paradox, and this bird provides an interesting model to examine how diet selection and the gut microbiota contribute to maximizing digestive efficiency. Therefore, we characterized the bacterial population in the crop of six adult hoatzins captured from the wild. A total of 1,235 16S rRNA gene sequences were grouped into 580 phylotypes (67% of the pooled species richness sampled, based on Good's coverage estimator, with C(ACE) and Chao1 estimates of 1,709 and 1,795 species-level [99% identity] operational taxonomic units, respectively). Members of 9 of the approximately 75 known phyla in Bacteria were identified in this gut habitat; the Firmicutes were dominant (67% of sequences, belonging to the classes Clostridia, Mollicutes, and Bacilli), followed by the Bacteroidetes (30%, mostly in the order Bacteroidales), Proteobacteria (1.8%), and Lentisphaerae, Verrucomicrobia, TM7, Spirochaetes, Actinobacteria, and Aminanaerobia (all <0.1%). The novelty in this ecosystem is great; 94% of the phylotypes were unclassified at the "species" level and thus likely include novel cellulolytic lineages.     

Comparative metagenomic analysis of bacterial populations in three full-scale mesophilic anaerobic manure digesters

While the use of anaerobic digestion to generate methane as a source of bioenergy is increasing worldwide, our knowledge of the microbial communities that perform biomethanation is very limited. Using next-generation sequencing, bacterial population profiles were determined in three full-scale mesophilic anaerobic digesters operated on dairy farms in the state of Vermont (USA). To our knowledge, this is the first report of a metagenomic analysis on the bacterial population of anaerobic digesters using dairy manure as their main substrate. A total of 20,366 non-chimeric sequence reads, covering the V1-V2 hypervariable regions of the bacterial 16S rRNA gene, were assigned to 2,176 operational taxonomic units (OTUs) at a genetic distance cutoff value of 5 %. Based on their limited sequence identity to validly characterized species, the majority of OTUs identified in our study likely represented novel bacterial species. Using a naïve Bayesian classifier, 1,624 anaerobic digester OTUs could be assigned to 16 bacterial phyla, while 552 OTUs could not be classified and may belong to novel bacterial taxonomic groups that have yet to be described. Firmicutes, Bacteroidetes, and Chloroflexi were the most highly represented bacteria overall, with Bacteroidetes and Chloroflexi showing the least and the most variation in abundance between digesters, respectively. All digesters shared 132 OTUs, which as a “core” group represented 65.4 to 70.6 % of sequences in individual digesters. Our results show that bacterial populations from microbial communities of anaerobic manure digesters can display high levels of diversity despite sharing a common core substrate.     

Novel high-rank phylogenetic lineages within a sulfur spring (Zodletone Spring, Oklahoma), revealed using a combined pyrosequencing-sanger approach

The utilization of high-throughput sequencing technologies in 16S rRNA gene-based diversity surveys has indicated that within most ecosystems, a significant fraction of the community could not be assigned to known microbial phyla. Accurate determination of the phylogenetic affiliation of such sequences is difficult due to the short-read-length output of currently available high-throughput technologies. This fraction could harbor multiple novel phylogenetic lineages that have so far escaped detection. Here we describe our efforts in accurate assessment of the novelty and phylogenetic affiliation of selected unclassified lineages within a pyrosequencing data set generated from source sediments of Zodletone Spring, a sulfide- and sulfur-rich spring in southwestern Oklahoma. Lineage-specific forward primers were designed for 78 putatively novel lineages identified within the pyrosequencing data set, and representative nearly full-length small-subunit (SSU) rRNA gene sequences were obtained by pairing those primers with reverse universal bacterial primers. Of the 78 lineages tested, amplifiable products were obtained for 52, 32 of which had at least one nearly full-length sequence that was representative of the lineage targeted. Analysis of phylogenetic affiliation of the obtained Sanger sequences identified 5 novel candidate phyla and 10 novel candidate classes (within Fibrobacteres, Planctomycetes, and candidate phyla BRC1, GN12, TM6, TM7, LD1, WS2, and GN06) in the data set, in addition to multiple novel orders and families. The discovery of multiple novel phyla within a pilot study of a single ecosystem clearly shows the potential of the approach in identifying novel diversities within the rare biosphere.     

metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data

The ribosomal rRNA genes are widely used as genetic markers for taxonomic identification of microbes. Particularly the small subunit (SSU; 16S/18S) rRNA gene is frequently used for species‐ or genus‐level identification, but also the large subunit (LSU; 23S/28S) rRNA gene is employed in taxonomic assignment. The metaxa software tool is a popular utility for extracting partial rRNA sequences from large sequencing data sets and assigning them to an archaeal, bacterial, nuclear eukaryote, mitochondrial or chloroplast origin. This study describes a comprehensive update to metaxa – metaxa 2 – that extends the capabilities of the tool, introducing support for the LSU rRNA gene, a greatly improved classifier allowing classification down to genus or species level, as well as enhanced support for short‐read (100 bp) and paired‐end sequences, among other changes. The performance of metaxa 2 was compared to other commonly used taxonomic classifiers, showing that metaxa 2 often outperforms previous methods in terms of making correct predictions while maintaining a low misclassification rate. metaxa 2 is freely available from http://microbiology.se/software/metaxa2/ .     

Molecular Evidence for the Broad Distribution of Anaerobic Ammonium-Oxidizing Bacteria in Freshwater and Marine Sediments

Previously available primer sets for detecting anaerobic ammonium-oxidizing (anammox) bacteria are inefficient, resulting in a very limited database of such sequences, which limits knowledge of their ecology. To overcome this limitation, we designed a new primer set that was 100% specific in the recovery of ~700-bp 16S rRNA gene sequences with >96% homology to the “Candidatus Scalindua” group of anammox bacteria, and we detected this group at all sites studied, including a variety of freshwater and marine sediments and permafrost soil. A second primer set was designed that exhibited greater efficiency than previous primers in recovering full-length (1,380-bp) sequences related to “Ca. Scalindua,” “Candidatus Brocadia,” and “Candidatus Kuenenia.” This study provides evidence for the widespread distribution of anammox bacteria in that it detected closely related anammox 16S rRNA gene sequences in 11 geographically and biogeochemically diverse freshwater and marine sediments.     

From Genus to Phylum: Large-Subunit and Internal Transcribed Spacer rRNA Operon Regions Show Similar Classification Accuracies Influenced by Database Composition

We compared the classification accuracy of two sections of the fungal internal transcribed spacer (ITS) region, individually and combined, and the 5′ section (about 600 bp) of the large-subunit rRNA (LSU), using a naive Bayesian classifier and BLASTN. A hand-curated ITS-LSU training set of 1,091 sequences and a larger training set of 8,967 ITS region sequences were used. Of the factors evaluated, database composition and quality had the largest effect on classification accuracy, followed by fragment size and use of a bootstrap cutoff to improve classification confidence. The naive Bayesian classifier and BLASTN gave similar results at higher taxonomic levels, but the classifier was faster and more accurate at the genus level when a bootstrap cutoff was used. All of the ITS and LSU sections performed well (>97.7% accuracy) at higher taxonomic ranks from kingdom to family, and differences between them were small at the genus level (within 0.66 to 1.23%). When full-length sequence sections were used, the LSU outperformed the ITS1 and ITS2 fragments at the genus level, but the ITS1 and ITS2 showed higher accuracy when smaller fragment sizes of the same length and a 50% bootstrap cutoff were used. In a comparison using the larger ITS training set, ITS1 and ITS2 had very similar accuracy classification for fragments between 100 and 200 bp. Collectively, the results show that any of the ITS or LSU sections we tested provided comparable classification accuracy to the genus level and underscore the need for larger and more diverse classification training sets.     

Analysis of gut microbiota in three species belonging to different genera (Hemitragus,Pseudois, and Ovis) from the subfamily Caprinae in the absence of environmental variance

This study aimed to identify the effects of host species on the gut microbial flora in three species (Hemitragus jemlahicus, Pseudois nayaur, and Ovis orientalis) from the subfamily Caprinae, by excluding the impact of environment factors. We investigated the differences in intestinal flora of three species belonging to Caprinae, which were raised in identical conditions. Fecal samples were collected from tahr, mouflon, and bharal, and the V3–V4 region of the 16S ribosomal RNA gene was analyzed by high‐throughput sequencing. The analysis of 16S rRNA gene sequences reveals that fecal samples were mainly composed of four phyla: Firmicutes, Bacteroidetes, Spirochaetes, and Proteobacteria. The most abundant phyla included Firmicutes and Bacteroidetes accounting for >90% of the bacteria, and a higher Firmicutes/Bacteroidetes ratio was observed in tahrs. Moreover, significant differences existed at multiple levels of classifications in the relative abundance of intestinal flora, differing greatly between species. Phylogenetic analyses based on 16S rRNA gene indicated that mouflon is closely related to bharal, and it is inconsistent with previous reports in the species evolutionary relationships. In this study, we demonstrated that the gut microbiota in tahr had a stronger ability to absorb and store energy from the diet compared with mouflon and bharal, and the characteristics of host–microbiome interactions were not significant.     

Detection of Stachybotrys chartarum using rRNA, tri5, and β-tubulin primers and determining their relative copy number by real-time PCR

Highly conserved regions are attractive targets for detection and quantitation by PCR, but designing species-specific primer sets can be difficult. Ultimately, almost all primer sets are designed based upon literature searches in public domain databases, such as the National Center for Biotechnology Information (NCBI). Prudence suggests that the researcher needs to evaluate as many sequences as available for designing species-specific PCR primers. In this report, we aligned 11, 9, and 16 DNA sequences entered for Stachybotrys spp. rRNA, tri5, and β-tubulin regions, respectively. Although we were able to align and determine consensus primer sets for the 9 tri5 and the 16 β-tubulin sequences, there was no consensus sequence that could be derived from alignment of the 11 rRNA sequences. However, by judicious clustering of the sequences that aligned well, we were able to design three sets of primers for the rRNA region of S. chartarum. The two primer sets for tri5 and β-tubulin produced satisfactory PCR results for all four strains of S. chartarum used in this study whereas only one rRNA primer set of three produced similar satisfactory results. Ultimately, we were able to show that rRNA copy number is approximately 2-log greater than for tri5 and β-tubulin in the four strains of S. chartarum tested.     

Effect of the Earthworms Lumbricus terrestris and Aporrectodea caliginosa on Bacterial Diversity in Soil

Earthworms ingest large amounts of soil and have the potential to radically alter the biomass, activity, and structure of the soil microbial community. In this study, the diversity of eight bacterial groups from fresh soil, gut, and casts of the earthworms Lumbricus terrestris and Aporrectodea caliginosa were studied by single-strand conformation polymorphism (SSCP) analysis using both newly designed 16S rRNA gene-specific primer sets targeting Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Bacteroidetes, Verrucomicrobia, Planctomycetes, and Firmicutes and a conventional universal primer set for SSCP, with RNA and DNA as templates. In parallel, the study of the relative abundance of these taxonomic groups in the same samples was performed using fluorescence in situ hybridization. Bacteroidetes, Alphaproteobacteria, and Betaproteobacteria were predominant in communities from the soil and worm cast samples. Representatives of classes Flavobacteria and Sphingobacteria (Bacteroidetes) and Pseudomonas spp. (low-abundant Gammaproteobacteria) were detected in soil and worm cast samples with conventional and taxon-targeting SSCP and through the sequence analysis of 16S rRNA clone libraries. Physiologically active unclassified Sphingomonadaceae (Alphaproteobacteria) and Alcaligenes spp. (Betaproteobacteria) also maintained their diversities during transit through the earthworm intestine and were found on taxon-targeting SSCP profiles from the soil and worm cast samples. In conclusion, our results suggest that some specific bacterial taxonomic groups maintain their diversity and even increase their relative numbers during transit through the gastrointestinal tract of earthworms.     

Members of the Candidate Division OP10 are spread in a variety of environments

Our planet holds a huge bacterial diversity. Most of these bacteria have only been detected by their 16S rRNA gene sequences remaining to be cultured. Many are classified within Candidate Divisions. One them is the Candidate Division OP10. Analysis of environmental 16S rRNA gene sequences available in public repositories revealed the existence of numerous sequences clustering within the Candidate Division OP10 but currently unclassified or assigned to other bacterial phyla. Newly proposed 16S rRNA sequences multiply several fold the reported sequences for the Candidate Division OP10. This study showed that the Candidate Division OP10 is a diverse and broadly distributed bacterial phylum and represents a stable microbial component in different natural environments.     

Discovery and characterization of a new bacterial candidate division by an anaerobic sludge digester metagenomic approach

We have constructed a large fosmid library from a mesophilic anaerobic digester and explored its 16S rDNA diversity using a high-density filter DNA–DNA hybridization procedure. We identified a group of 16S rDNA sequences forming a new bacterial lineage named WWE3 (Waste Water of Evry 3). Only one sequence from the public databases shares a sequence identity above 80% with the WWE3 group which hence cannot be affiliated to any known or candidate prokaryotic division. Despite representing a non-negligible fraction (5% of the 16S rDNA sequences) of the bacterial population of this digester, the WWE3 bacteria could not have been retrieved using the conventional 16S rDNA amplification procedure due to their unusual 16S rDNA gene sequence. WWE3 bacteria were detected by polymerase chain reaction (PCR) in various environments (anaerobic digesters, swine lagoon slurries and freshwater biofilms) using newly designed specific PCR primer sets. Fluorescence in situ hybridization (FISH) analysis of sludge samples showed that WWE3 microorganisms are oval-shaped and located deep inside sludge flocs. Detailed phylogenetic analysis showed that WWE3 bacteria form a distinct monophyletic group deeply branching apart from all known bacterial divisions. A new bacterial candidate division status is proposed for this group.     

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.