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.     

Comparison of bacterial diversity along the human intestinal tract by direct cloning and sequencing of 16S rRNA genes

Bacterial diversity of the mucosal biopsies from human jejunum, distal ileum, ascending colon and rectum were compared by analysis of PCR-amplified 16S rDNA clone libraries. A total of 347 clones from the mucosal biopsies were partially sequenced and assigned to six phylogenetic phyla of the domain Bacteria: Firmicutes, Bacteroidetes, Proteobacteria, Fusobacteria, Verrucomicrobia, and Actinobacteria. The jejunum sample had least microbial diversity compared to the other samples and a trend towards highest diversity in ascending colon was observed. The clone libraries of distal ileum, ascending colon and rectum were not significantly different from each other (P>0.0043), but they differed significantly from the jejunum library (P=0.001). The population of sequences retrieved from jejunal biopsies was dominated by sequences closely related to Streptococcus (67%), while the population of sequences derived from distal ileum, ascending colon and rectum were dominated by sequences affiliated with Bacteroidetes (27-49%), and Clostridium clusters XIVa (20-34%) and IV (7-13%). The results indicate that the microbial community in jejunum is different from those in distal ileum, ascending colon and rectum, and that the major phylogenetic groups are similar from distal ileum to rectum.     

Subgrouping of bacterial populations by cellular fatty acid composition

Abstract The cellular fatty acid composition of six bacterial species isolated from the seeds and leaves of sugar beet ( Beta vulgaris ) and from soil were analysed. The quantitative data from the fatty acid methyl ester (FAME) profiles were highly reproducible. Numerical analysis of Xanthomonas maltophilia . FAME profiles sub-grouped strains according to when they were isolated in the growing season. The analytical method used was sensitive enough to differentiate strains of Klebsiella terrigena isolated from either soil or leaves. The results from this study confirm reports that analyses of bacterial FAME composition were rapid to perform, specific and allowed differentiation of strains within the same species.     

Sequence-specific bacterial growth inhibition by peptide nucleic acid targeted to the mRNA binding site of 16S rRNA

To compare oxidative dissolution rates of chalcopyrite by different consortia of moderately thermophilic acidophiles, various defined mixed cultures of three bacteria Acidithiobacillus caldus s2, Leptospirillum ferriphilum YSK, and Sulfobacillus sp. LN and one archaeon Ferroplasma thermophilum L1 were studied in batch shake flask cultures incubated at 45 °C. Chalcopyrite dissolution was determined by measuring variations of soluble copper, ferric iron, and pH. Microbial population dynamics involved in bioleaching process were monitored using real-time quantitative polymerase chain reaction (PCR) technology. The complex consortia containing both chemoautotrophic (L. ferriphilum and At. caldus) and chemomixotrophic (Sulfobacillus LN and F. thermophilum) moderate thermophiles were found to be the most efficient in all of those tested. Mutualistic interactions between physiologically distinct moderately thermophilic acidophiles, involving transformations of iron and sulfur and transfer of organic compound, were considered to play a critical role in promoting chalcopyrite dissolution. The real-time PCR assay was reliable to analyze population dynamics of moderate thermophiles in bioleaching systems, and the analysis results were consistent with physiological characteristics of these strains.     

Investigation of chicken intestinal bacterial communities by 16S rRNA targeted fluorescence in situ hybridization

The aim of the investigation was to quantify selected dominant bacterial groups in the chicken intestinal tract. Conventional production was used as model and the effect of the supplement with Salinomycin was evaluated. Hybridization conditions were optimized for published probes with respect to a panel of reference bacteria. In chicken intestinal samples bacteria were quantified by fluorescence in situ hybridization with 16S rRNA oligonucleotides directed towards bacteria related to Lactobacillus, Bacillus, Enterococcus-Streptococcus-Lactococcus, Enterobacteriaceae, Bacteroides, Clostridium and the domain Bacteria in lumen of ileum and cecum as well as on the intestinal wall including mucus of four individuals. Salinomycin in feed reduced counts of the Lactobacillus-, Enterobacteriaceae- and Clostridium-like bacteria in lumen of ileum compared to the conventional control. Increased or decreased bacterial counts were registered by Salinomycin in the ceca compared to the control. Relatively higher counts of Bacteroides- and Clostridium-like bacteria were found on the intestinal wall including mucus compared to lumen. The increase in numbers of some bacterial groups as well as the expected reduction by Salinomycin and the observed difference in the relative distribution of bacteria between lumen and intestinal wall are new observations that will need further investigation.     

Analysis of diversity and activity of sulfate-reducing bacterial communities in sulfidogenic bioreactors using 16S rRNA and dsrB genes as molecular markers

Here we describe the diversity and activity of sulfate-reducing bacteria (SRB) in sulfidogenic bioreactors by using the simultaneous analysis of PCR products obtained from DNA and RNA of the 16S rRNA and dissimilatory sulfite reductase (dsrAB) genes. We subsequently analyzed the amplified gene fragments by using denaturing gradient gel electrophoresis (DGGE). We observed fewer bands in the RNA-based DGGE profiles than in the DNA-based profiles, indicating marked differences in the populations present and in those that were metabolically active at the time of sampling. Comparative sequence analyses of the bands obtained from rRNA and dsrB DGGE profiles were congruent, revealing the same SRB populations. Bioreactors that received either ethanol or isopropanol as an energy source showed the presence of SRB affiliated with Desulfobulbus rhabdoformis and/or Desulfovibrio sulfodismutans, as well as SRB related to the acetate-oxidizing Desulfobacca acetoxidans. The reactor that received wastewater containing a diverse mixture of organic compounds showed the presence of nutritionally versatile SRB affiliated with Desulfosarcina variabilis and another acetate-oxidizing SRB, affiliated with Desulfoarculus baarsii. In addition to DGGE analysis, we performed whole-cell hybridization with fluorescently labeled oligonucleotide probes to estimate the relative abundances of the dominant sulfate-reducing bacterial populations. Desulfobacca acetoxidans-like populations were most dominant (50 to 60%) relative to the total SRB communities, followed by Desulfovibrio-like populations (30 to 40%), and Desulfobulbus-like populations (15 to 20%). This study is the first to identify metabolically active SRB in sulfidogenic bioreactors by using the functional gene dsrAB as a molecular marker. The same approach can also be used to infer the ecological role of coexisting SRB in other habitats.     

Analysis of cellular fatty acids in Orientia tsutsugamushi as taxonomic markers

Six Orientia strains including 3 prototype strains such as Gilliam, Karp, and Kato, and 3 strains (Boryong, Pajoo, and Yongworl) isolated in Korea, were studied for the profiles of their cellular fatty acids. All tested strains contained octadecenoic acid C (18: 1) omega 9 c(57.3 +/- 3.5%), octadecanoic acid C (18: 0) (15.3 +/- 1.5%), and hexadecanoic acid C (16: 0) (12.7 +/- 1.7%) as major components; however, interestingly, eicosenoic acid C (20: 1) omega 9 c(2.6 +/- 0.6%) was found in all strains except the Yongworl strain. Furthermore none of the strains contained 3-hydroxy fatty acids. The ratios of total saturated fatty acid (SFA) to total unsaturated fatty acid (UFA) were within the range of 0.34 to 0.54. These results showed that the cellular fatty acid profile should provide more reliable information for the identification of these bacteria.     

Polymerase chain reaction detection of bacterial 16S rRNA gene in human blood

Bacterial 16S ribosomal RNA genes (rDNA) were detected in blood samples from two healthy individuals by PCR under conditions involving 30 cycles that did not produce any visible products from negative control saline. Even from control samples, PCR involving 35-40 cycles yielded visible bands. Major clones detected in the blood samples, but not in control, were the Aquabacterium subgroup, Stenotrophomonas subgroup, Budvicia subgroup, Serratia subgroup, Bacillus subgroup and Flavobacteria subgroup. No clone was located within the bacteroides-clostridium-lactobacillus cluster, which is indigenous to gastrointestinal flora.     

Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys

Taxonomic classification of the thousands–millions of 16S rRNA gene sequences generated in microbiome studies is often achieved using a naïve Bayesian classifier (for example, the Ribosomal Database Project II (RDP) classifier), due to favorable trade-offs among automation, speed and accuracy. The resulting classification depends on the reference sequences and taxonomic hierarchy used to train the model; although the influence of primer sets and classification algorithms have been explored in detail, the influence of training set has not been characterized. We compared classification results obtained using three different publicly available databases as training sets, applied to five different bacterial 16S rRNA gene pyrosequencing data sets generated (from human body, mouse gut, python gut, soil and anaerobic digester samples). We observed numerous advantages to using the largest, most diverse training set available, that we constructed from the Greengenes (GG) bacterial/archaeal 16S rRNA gene sequence database and the latest GG taxonomy. Phylogenetic clusters of previously unclassified experimental sequences were identified with notable improvements (for example, 50% reduction in reads unclassified at the phylum level in mouse gut, soil and anaerobic digester samples), especially for phylotypes belonging to specific phyla (Tenericutes, Chloroflexi, Synergistetes and Candidate phyla TM6, TM7). Trimming the reference sequences to the primer region resulted in systematic improvements in classification depth, and greatest gains at higher confidence thresholds. Phylotypes unclassified at the genus level represented a greater proportion of the total community variation than classified operational taxonomic units in mouse gut and anaerobic digester samples, underscoring the need for greater diversity in existing reference databases.     

Alterations in growth and fatty acid profiles under stress conditions of Hansenula polymorpha defective in polyunsaturated fatty acid synthesis

Using chemical mutagenesis, mutants of Hansenula polymorpha that were defective in fatty acid synthesis were selected based on their growth requirements on saturated fatty acid mixtures. One mutant (S7) was incapable of synthesizing polyunsaturated fatty acids (PUFA), linoleic and α-linolenic acids. A genetic analysis demonstrated that the S7 strain had a double lesion affecting fatty acid synthesis and Δ¹²-desaturation. A segregant with a defect in PUFA synthesis (H69-2C) displayed normal growth characteristics in the temperature range of 20–42 °C through a modulation of the cellular fatty acid composition. Compared with the parental strain, this yeast mutant had increased sensitivity at low and high temperatures (15 and 48 °C, respectively) with an increased tolerance to oxidative stress. The responses to ethanol stress were similar for the parental and PUFA-defective strains. Myristic acid was also determined to play an essential role in the cell growth of H. polymorpha. These findings suggest that both the type of cellular fatty acids and the composition of fatty acids might be involved in the stress responsive mechanisms in this industrially important yeast.     

Evaluation of the 16S and 12S rRNA genes as universal markers for the identification of commercial fish species in South Africa

The development of DNA-based methods for the identification of fish species is important for fisheries research and control, as well as for the detection of unintentional or fraudulent species substitutions in the marketplace. The aim of this study was to generate a comprehensive reference database of DNA sequences from the mitochondrial 16S and 12S ribosomal RNA (rRNA) genes for 53 commercial fish species in South Africa and to evaluate the applicability of these genetic markers for the identification of fish at the species level. The DNA extracted from all target species was readily amplified using universal primers targeting both rRNA gene regions. Sequences from the 16S and 12S rRNA genes were submitted to GenBank for the first time for 34% and 53% of the fish species, respectively. Cumulative analysis of the 16S rRNA gene sequences revealed mean conspecific, congeneric and confamilial Kimura two parameter (K2P) distances of 0.03%, 0.70% and 5.10% and the corresponding values at the 12S level were 0.03%, 1.00% and 5.57%. K2P neighbour-joining trees based on both sequence datasets generally clustered species in accordance with their taxonomic classifications. The nucleotide variation in both the 16S and 12S sequences was suitable for identifying the large majority of the examined fish specimens to at least the level of genus, but was found to be less useful for the explicit differentiation of certain congeneric fish species. It is recommended that one or more faster-evolving DNA regions be analysed to confirm the identities of closely-related fish species in South Africa.     

Evaluation of broiler litter with reference to the microbial composition as assessed by using 16S rRNA and functional gene markers

Very little is known about the microbial composition of animal bedding wastes, including poultry litter, and what is known has been deduced from standard culture methods, by which some fastidious organisms that exist in the environment may not be detected. We evaluated the bacterial composition of poultry litter by using a combination of culture and molecular detection. Total aerobic bacteria in poultry litter were detected by culture at 10(9) CFU/g of material. Enteric bacteria such as Enterococcus spp. and coliforms composed 0.1 and 0.01%, respectively, of the total aerobic cultivatable bacteria in poultry litter; no Salmonella strains were detected by culture. In order to characterize the most abundant bacterial groups, we sequenced 16S ribosomal DNA (rDNA) genes amplified by PCR with microbial community DNA isolated from poultry litter as the template. From the 16S rDNA library, 31 genera were identified. Twelve families or groups were identified with lactobacilli and Salinococcus spp. forming the most abundant groups. In fact, 82% of the total sequences were identified as gram-positive bacteria with 62% of total belonging to low G+C gram-positive groups. In addition to detection of 16S rDNA sequences associated with the expected fecal bacteria present in manure, we detected many bacterial sequences for organisms, such as Globicatella sulfidofaciens, Corynebacterium ammoniagenes, Corynebacterium urealyticum, Clostridium aminovalericum, Arthrobacter sp., and Denitrobacter permanens, that may be involved in the degradation of wood and cycling of nitrogen and sulfur. Several sequences were identified in the library for bacteria associated with disease in humans and poultry such as clostridia, staphylococci, and Bordetella spp. However, specific PCR targeting other human and veterinary pathogens did not detect the presence of Salmonella, pathogenic Escherichia coli, Campylobacter spp., Yersinia spp., Listeria spp., or toxigenic staphylococci. PCR and DNA hybridization revealed the presence of class 1 integrons with gene cassettes that specify resistance to aminoglycosides and chloramphenicol. Only from understanding the microbial community of animal wastes such as poultry litter can we manage animal disease and limit the impact of animal waste on the environment and human and animal health.     

Distinctness of spore and vegetative cellular fatty acid profiles of some aerobic endospore-forming bacilli

A gas chromatographic analysis method was employed to determine the cellular fatty acid (CFA) profiles of spores and vegetative cells of some aerobic endospore-forming bacilli. The harvests of experimental strains were processed to obtain pure spores and acquire whole cell fatty acid methyl esters for the subsequent gas chromatographic analysis, and the corresponding vegetative cells were set as control. Evaluation of reproducibility of spore CFA components revealed that, provided under standardized experimental procedure, spore CFA composition was stable enough for research purposes. Fatty acids recovered in spores in greater quantities were saturated branched-chain acids containing 15 and 17 carbon atoms, similar to the vegetative cells. Commonly, the proportions of saturated branched-chain acids in spores were greater than in vegetative cells. The dendrograms obtained by cluster analysis provided some meaningful taxonomic information of the experimental strains. The fatty acids analysis of spores seems to be a promising supplementary tool for the chemotaxonomic research of aerobic endospore-forming bacilli.     

Characterization of Bacteroides gingivalis by direct fluorescent antibody staining and cellular fatty acid profiles

Bacteroides gingivalis is a newly proposed species which includes strains isolated from the mouth. Thirteen strains of B. gingivalis isolated from three geographic locations in the United States and France were examined with direct fluorescent antibody staining and analysis of total cellular fatty acids and compared with 16 strains of B. asaccharolyticus of nonoral origin by the same methods. Bacteroides gingivalis strains reacted with the B. gingivalis conjugate (fluorescein isothiocyanate labeled antibody reagent) only, while the B, asaccharolyticus strains reacted with the B. asaccharolyticus conjugate only. The B. gingivalis strains showed negative fluorescence with fluorescein isothiocyanate conjugates for other black-pigmented Bacteroides species. The specificity of the B. gingivalis conjugate was demonstrated by its failure to stain 88 strains of aerobic and anaerobic bacteria other than B. gingivalis. The fatty acid profiles of B. gingivalis and B. asaccharolyticus were readily distinguishable. The B. gingivalis profile was also distinguishable from those of other pigmenting Bacteroides species on the basis of concentration ratios among the characteristic components. These results support the species separation of B. gingivalis and B. asaccharolyticus. Further, they indicate the usefulness of cellular fatty acid profiles as an adjunct to the use of specific fluorescent antibody conjugates for identification of Bacteroides species.     

Identification of bacterial populations in dairy wastewaters by use of 16S rRNA gene sequences and other genetic markers

Hydraulic flush waste removal systems coupled to solid/liquid separators and circulated treatment lagoons are commonly utilized to manage the large amounts of animal waste produced on high-intensity dairy farms. Although these systems are common, little is known about the microbial populations that inhabit them or how they change as they traverse the system. Using culture-based and non-culture-based methods, we characterized the microbial community structure of manure, water from the separator pit, and water from the circulated treatment lagoon from a large dairy in the San Joaquin Valley of California. Our results show that both total bacterial numbers and bacterial diversity are highest in manure, followed by the separator pit water and the lagoon water. The most prevalent phylum in all locations was the Firmicutes (low-G+C, gram-positive bacteria). The most commonly occurring operational taxonomic unit (OTU) had a 16S rRNA gene (rDNA) sequence 96 to 99% similar to that of Clostridium lituseburense and represented approximately 6% of the manure derived sequences, 14% of the separator pit-derived sequences and 20% of the lagoon-derived sequences. Also highly prevalent was an OTU with a 16S rDNA sequence 97 to 100% similar to that of Eubacterium tenue, comprising approximately 3% of the manure-derived sequences, 6% of the separator pit-derived sequences and 9% of the lagoon-derived sequences. Taken together, these sequences represent approximately one-third of the total organisms in the lagoon waters, suggesting that they are well adapted to this environment.     

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.     

Key bacterial taxa and metabolic pathways affecting gut short-chain fatty acid profiles in early life

Infant gut microbiota development affects the host physiology throughout life, and short-chain fatty acids (SCFAs) are promising key metabolites mediating microbiota-host relationships. Here, we investigated dense longitudinally collected faecal samples from 12 subjects during the first 2 years (n = 1048) to identify early life gut SCFA patterns and their relationships with the microbiota. Our results revealed three distinct phases of progression in the SCFA profiles: early phase characterised by low acetate and high succinate, middle-phase characterised by high lactate and formate and late-phase characterised by high propionate and butyrate. Assessment of the SCFA–microbiota relationships revealed that faecal butyrate is associated with increased Clostridiales and breastfeeding cessation, and that diverse and personalised assemblage of Clostridiales species possessing the acetyl-CoA pathway play major roles in gut butyrate production. We also found an association between gut formate and some infant-type bifidobacterial species, and that human milk oligosaccharides (HMO)-derived fucose is the substrate for formate production during breastfeeding. We identified genes upregulated in fucose and fucosylated HMO utilisation in infant-type bifidobacteria. Notably, bifidobacteria showed interspecific and intraspecific variation in the gene repertoires, and cross-feeding of fucose contributed to gut formate production. This study provides an insight into early life SCFA–microbiota relationships, which is an important step for developing strategies for modulating lifelong health.Subject terms: Microbial ecology, Microbiome, Bacterial genetics