Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Comparative Analysis of Microbial Profiles in Cow Rumen Fed with Different Dietary Fiber by Tagged 16S rRNA Gene Pyrosequencing

The ruminal microbiome of cattle plays an important role not only in animal health and productivity but also in food safety and environment. Microbial profiles of rumen fluid obtained from dairy cows fed on three different fiber/starch diet compositions were characterized. Tagged 16S rRNA gene pyrosequencing and statistical analysis revealed that the dominant ruminal bacteria shared by all three sample groups belonged to phyla Bacteroidetes, Firmicutes, and Proteobacteria. However, the relative abundance of these bacterial groups was markedly affected by diet composition. In animals fed with a high fiber diet, the fibrolytic and cellulolytic bacteria Lachnospiraceae, Ruminococcaceae, and Fibrobacteraceae were found in highest abundance compared with animals fed other diets with lower fiber content. The polysaccharide-degrading Prevotellaceae and Flavobacteriaceae bacteria were most abundant in the rumen of cows fed on diet with the highest starch content. These data highlight the ruminal microbiome’s ability to adapt to feed composition and also provide a basis for the development of feed formulation systems designed to improve livestock productivity.     

Phylogenetic diversity and dietary association of rumen Treponema revealed using group-specific 16S rRNA gene-based analysis

Treponema spp. are a commonly detected bacterial group in the rumen that are involved in the degradation of soluble fibers. In this study, a ruminal Treponema group-specific PCR primer targeting the 16S rRNA gene was designed and used to assess the phylogenetic diversity and diet association of this group in sheep rumen. Total DNA was extracted from rumen digesta of three sheep fed a diet based on alfalfa/orchardgrass hay or concentrate. The real-time PCR quantification indicated that the relative abundance of the Treponema group in the total rumen bacteria was as high as 1.05%, while the known species Treponema bryantii accounted for only 0.02%. Fingerprints of the Treponema community determined by 16S rDNA-targeted denaturing gradient gel electrophoresis (DGGE) analysis tended to differ among the diets. Principal component analysis of the DGGE profiles distinguished those Treponema associated with either the hay or the concentrate diets. Analysis of a Treponema 16S rRNA gene clone library showed phylogenetically distinct operational taxonomic units for a specific dietary condition, and significant (P=0.001) differences in community composition were observed among clone libraries constructed from each dietary regimen. The majority of clones (75.4%) had <97% sequence similarity with known Treponema. These results suggest the predominance of uncultured Treponema that appear to have distinct members related to the digestion of either hay or concentrate diet.     

Adenosine Triphosphate and Deoxyribonucleic Acid in the Alimentary Tract of Cattle Fed Different Nitrogen Sources

Total microbial cell counts and the content of deoxyribonucleic acid (DNA) and adenosine triphosphate (ATP) were determined in ruminal and abomasal digesta and in faeces from heifers fed a basal diet of barley and straw supplemented with urea, casein, soybean-protein concentrate or feather meal. The concentration of ATP in the rumen content varied independently of the DNA content and total cell count, but dependently upon the nitrogen sources, the highest ATP/DNA ratio being obtained with casein. The ATP/DNA ratio in faeces was only one-tenth of that found in the digesta of the rumen and abomasum, indicating either an extremely low level of activity of the microbial cells or more probably a very large number of dead organisms. It is suggested that DNA and ATP might be useful indices of the microbial status in terms of biomass and metabolic activity. The total cell count should still be included in routine studies to determine the proportion of protozoa to bacteria.     

Detection of Clostridium proteoclasticum and Closely Related Strains in the Rumen by Competitive PCR

A competitive PCR technique was used to enumerate the proteolytic bacterium Clostridium proteoclasticum from the rumen. A PCR primer, which circumscribes this organism and several closely related strains, was designed for a variable region within their 16S rRNA genes and was used in conjunction with a universal forward primer. This primer pair was tested for specificity against 85 ruminal bacterial strains. An internal control DNA was constructed for use in competitive PCRs and was shown to amplify under the same reaction conditions and with the same amplification efficiency as the target DNA. DNA from a known number of C. proteoclasticum cells was coamplified with the internal control to construct a standard curve. Rumen samples were collected from eight dairy cows fed four diets in rotation: high nitrogen, high nitrogen supplemented with carbohydrate, low nitrogen, and low nitrogen supplemented with carbohydrate. DNA extracted from these and spiked with internal control DNA was amplified with the C. proteoclasticum primer pair. The relative intensities of the PCR products were used to quantitate the numbers of C. proteoclasticum cell equivalents from the rumen samples. The numbers ranged from 2.01 × 10⁶ ml⁻¹ to 3.12 × 10⁷ ml⁻¹. There was no significant effect on the numbers of C. proteoclasticum detected in rumen samples among cows fed the four diets. The utility of the competitive PCR approach for quantifying ruminal bacterial populations in vivo and the occurrence of C. proteoclasticum in forage-fed dairy cows are discussed.     

Natural variability and diurnal fluctuations within the bacteriophage population of the rumen.

To investigate the impact of nutritional and environmental factors on bacteriophage activity in the rumen, it is first valuable to determine the extent of natural variations and fluctuations in phage populations from different animal species, and from animals located together and separately, and variation in animals over time. Differences in phage populations between sheep on different diets, between sheep and goats, and within the rumen over time were investigated by using pulsed-field gel electrophoresis and comparing total phage DNA in ruminal fluid. It was found that no two individuals had similar DNA banding patterns, even when similarly fed and penned together, indicating there is considerable individual diversity in phage populations between animals. Despite these individual differences, the quantities, but not the banding patterns, of phage DNA were similar for animals within groups but varied between groups, suggesting that nutritional factors may influence overall phage activity in the rumen. In sheep fed once daily, a distinct diurnal variation in the phage population was observed. Two hours postfeeding, total phage DNA dropped to its lowest level. The phage population then increased, reaching a maximal level 8 to 10 h postfeeding before declining over the next 4 h to reach a stable concentration for the rest of the cycle. The general trend in phage DNA concentration appeared similar to previously recorded diurnal fluctuations in ruminal bacterial populations in cattle fed once daily.     

Phylogenetic diversity of gene sequences isolated from the rumen as analysed using a self‐organizing map (SOM)

Aims: To determine the origins of DNA sequences isolated from the rumen microbial ecosystem using a self‐organizing map (SOM). Methods and Results: DNA sequences other than 16S small subunit ribosomal RNA (SSU rRNA) gene sequences that were detected from the rumen were analysed by the SOM method reported by Abe et al. [2000, Self‐Organizing Map (SOM) unveils and visualizes hidden sequence characteristics of a wide range of eukaryote genomes. Gene 365, 27–34]. Because query sequences positioned by SOM were scattered on the master drawing of SOM, it was suggested that many DNA sequences isolated from the rumen were collected from a broad range of micro‐organisms. Although the results obtained by SOM were similar to those obtained by the neighbour‐joining (NJ) method, SOM was able to presume the phylotypes of the query sequences without information about the 16S SSU rRNA gene sequences and homology searches, and to reveal existence of novel micro‐organisms deduced to be cellulolytic bacteria, archaea and methanotrophic bacterium. Conclusions: As the SOM method defined phylotypes of unreported rumen micro‐organisms, it is presumed that these phylotypes would be involved in rumen fermentation in cooperation with known rumen micro‐organisms. Moreover, it is demonstrated that SOM is a useful tool for affiliating DNA sequences, which have no matches in databases. Significance and Impact of Study: Through SOM analysis, a better means of identifying rumen micro‐organisms and estimating their roles in rumen function was provided.     

Flaxseed supplementation decreases methanogenic gene abundance in the rumen of dairy cows

The objective of this study was to investigate the effects of a flaxseed-supplemented diet on archaeal abundance and gene expression of methanogens in the rumen of dairy cows. In all, 11 non-lactating dairy cows were randomly divided into two groups: group A (five cows) and B (six cows). The two diets fed were: (1) the control diet, a conventional dry cow ration; and (2) the flaxseed-supplemented diet, the conventional dry cow ration adjusted with 12.16% ground flaxseed incorporated into the total mixed ration. A cross-over experiment was performed with the two groups of cows fed the two different diets for five 21-day periods, which included the first adaptation period followed by two treatment and two wash out periods. At the end of each feeding period, rumen fluid samples were collected via rumenocentesis and DNA was extracted. Quantitative PCR was utilized to analyze the gene abundance of 16S ribosomal RNA (16S rRNA) targeting the ruminal archaea population and the mcrA gene coding for methyl coenzyme-M reductase subunit A, a terminal enzyme in the methanogenesis pathway. Results demonstrated a 49% reduction of 16S rRNA and 50% reduction of mcrA gene abundances in the rumen of dairy cows fed the flaxseed-supplemented diet in comparison with those fed the control diet. This shows flaxseed supplementation effectively decreases the methanogenic population in the rumen. Future studies will focus on the mechanisms for such reduction in the rumen of dairy cattle, as well as the relationship between methanogenic gene expression and methane production.     

Rumen Microbiome from Steers Differing in Feed Efficiency

The cattle rumen has a diverse microbial ecosystem that is essential for the host to digest plant material. Extremes in body weight (BW) gain in mice and humans have been associated with different intestinal microbial populations. The objective of this study was to characterize the microbiome of the cattle rumen among steers differing in feed efficiency. Two contemporary groups of steers (n=148 and n=197) were fed a ration (dry matter basis) of 57.35% dry-rolled corn, 30% wet distillers grain with solubles, 8% alfalfa hay, 4.25% supplement, and 0.4% urea for 63 days. Individual feed intake (FI) and BW gain were determined. Within contemporary group, the four steers within each Cartesian quadrant were sampled (n=16/group) from the bivariate distribution of average daily BW gain and average daily FI. Bacterial 16S rRNA gene amplicons were sequenced from the harvested bovine rumen fluid samples using next-generation sequencing technology. No significant changes in diversity or richness were indicated, and UniFrac principal coordinate analysis did not show any separation of microbial communities within the rumen. However, the abundances of relative microbial populations and operational taxonomic units did reveal significant differences with reference to feed efficiency groups. Bacteroidetes and Firmicutes were the dominant phyla in all ruminal groups, with significant population shifts in relevant ruminal taxa, including phyla Firmicutes and Lentisphaerae, as well as genera Succiniclasticum, Lactobacillus, Ruminococcus, and Prevotella. This study suggests the involvement of the rumen microbiome as a component influencing the efficiency of weight gain at the 16S level, which can be utilized to better understand variations in microbial ecology as well as host factors that will improve feed efficiency.     

PCR-DGGE analysis reveals a distinct diversity in the bacterial population attached to the rumen epithelium

Bacteria attached to the rumen epithelium (or epimural community) are not well characterised and their role in rumen functioning is not totally understood. There is just one published report of a clone library from one cow that suggests that this epimural community differs from the bacteria associated with the rumen digestive contents. However, this time-consuming approach is not adapted for examining microbial population changes in groups of animals. In in vivo studies, when samples from several animals have to be analysed simultaneously, a simpler technique has to be used. In this study, a genetic fingerprinting technique, polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE), was used to characterise the structure of the bacterial population attached to the rumen epithelium. This community was compared with that present in the solid and liquid phases of rumen content under two contrasting diets. Rumen samples were obtained from four forage-fed and four high-concentrate-fed (80 : 20, wheat grain : hay) 5-month-old lambs. After slaughter, samples from five epithelial sites and the solid and liquid digesta phases were taken for DNA extraction and analysis. Bacterial communities were profiled by PCR-DGGE using bacterial-specific 16S rDNA primers. Analysis of the fingerprint revealed that the epithelial community differed from those of rumen content in both diets. As expected, the nature of the feed influenced the bacterial communities from the solid and liquid rumen phases but no diet effect was observed in the rumen epithelial profiles suggesting a strong host effect on this bacterial population. Additionally, no differences were observed among the five epithelial sampling sites taken from each animal. The profile of the bacterial population attached to the rumen epithelium presented a high inter-animal variation, whether this difference has an influence in the function of this community remains to be determined.     

Advances in microbial ecosystem concepts and their consequences for ruminant agriculture

Microbial transformations in the rumen ecosystem have a major impact on our ability to meet the challenge of reducing the environmental footprint of ruminant livestock agriculture, as well as enhancing product quality. Current understanding of the rumen microbial ecosystem is limited, and affects our ability to manipulate rumen output. The view of ruminal fermentation as the sum of activities of the dominant rumen microbiota is no longer adequate, with a more holistic approach required. This paper reviews rumen functionality in the context of the microbiota of the rumen ecosystem, addressing ruminal fermentation as the product of an ecosystem while highlighting the consequences of this for ruminant agriculture. Microbial diversity in the rumen ecosystem enhances the resistance of the network of metabolic pathways present, as well as increasing the potential number of new pathways available. The resulting stability of rumen function is further promoted by the existence of rumen microbiota within biofilms. These protected, structured communities offer potential advantages, but very little is currently known about how ruminal microorganisms interact on feed-surfaces and how these communities develop. The temporal and spatial development of biofilms is strongly linked to the availability of dietary nutrients, the dynamics of which must also be given consideration, particularly in fresh-forage-based production systems. Nutrient dynamics, however, impact not only on pathway inputs but also the turnover and output of the whole ecosystem. Knowledge of the optimal balance of metabolic processes and the corresponding microbial taxa required to provide a stable, balanced ecosystem will enable a more holistic understanding of the rumen. Future studies should aim to identify key ecosystem processes and components within the rumen, including microbial taxa, metabolites and plant-based traits amenable to breeding-based modification. As well as gaining valuable insights into the biology of the rumen ecosystem, this will deliver realistic and appropriate novel targets for beneficial manipulation of rumen function.     

Oral Samples as Non-Invasive Proxies for Assessing the Composition of the Rumen Microbial Community

Microbial community analysis was carried out on ruminal digesta obtained directly via rumen fistula and buccal fluid, regurgitated digesta (bolus) and faeces of dairy cattle to assess if non-invasive samples could be used as proxies for ruminal digesta. Samples were collected from five cows receiving grass silage based diets containing no additional lipid or four different lipid supplements in a 5 x 5 Latin square design. Extracted DNA was analysed by qPCR and by sequencing 16S and 18S rRNA genes or the fungal ITS1 amplicons. Faeces contained few protozoa, and bacterial, fungal and archaeal communities were substantially different to ruminal digesta. Buccal and bolus samples gave much more similar profiles to ruminal digesta, although fewer archaea were detected in buccal and bolus samples. Bolus samples overall were most similar to ruminal samples. The differences between both buccal and bolus samples and ruminal digesta were consistent across all treatments. It can be concluded that either proxy sample type could be used as a predictor of the rumen microbial community, thereby enabling more convenient large-scale animal sampling for phenotyping and possible use in future animal breeding programs aimed at selecting cattle with a lower environmental footprint.     

Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle

Feed-efficient animals have lower production costs and reduced environmental impact. Given that rumen microbial fermentation plays a pivotal role in host nutrition, the premise that rumen microbiota may contribute to host feed efficiency is gaining momentum. Since diet is a major factor in determining rumen community structure and fermentation patterns, we investigated the effect of divergence in phenotypic residual feed intake (RFI) on ruminal community structure of beef cattle across two contrasting diets. PCR-denaturing gradient gel electrophoresis (DGGE) and quantitative PCR (qPCR) were performed to profile the rumen bacterial population and to quantify the ruminal populations of Entodinium spp., protozoa, Fibrobacter succinogenes, Ruminococcus flavefaciens, Ruminococcus albus, Prevotella brevis, the genus Prevotella, and fungi in 14 low (efficient)- and 14 high (inefficient)-RFI animals offered a low-energy, high-forage diet, followed by a high-energy, low-forage diet. Canonical correspondence and Spearman correlation analyses were used to investigate associations between physiological variables and rumen microbial structure and specific microbial populations, respectively. The effect of RFI on bacterial profiles was influenced by diet, with the association between RFI group and PCR-DGGE profiles stronger for the higher forage diet. qPCR showed that Prevotella abundance was higher (P < 0.0001) in inefficient animals. A higher (P < 0.0001) abundance of Entodinium and Prevotella spp. and a lower (P < 0.0001) abundance of Fibrobacter succinogenes were observed when animals were offered the low-forage diet. Thus, differences in the ruminal microflora may contribute to host feed efficiency, although this effect may also be modulated by the diet offered.     

Coupling 16S-ITS rDNA clone libraries and automated ribosomal intergenic spacer analysis to show marine microbial diversity: development and application to a time series

We outline an approach to simultaneously assess multilevel microbial diversity patterns utilizing 16S-ITS rDNA clone libraries coupled with automated ribosomal intergenic spacer analysis (ARISA). Sequence data from 512 clones allowed estimation of ARISA fragment lengths associated with bacteria in a coastal marine environment. We matched 92% of ARISA peaks (each comprising >1% total amplified product) with corresponding lengths from clone libraries. These peaks with putative identification accounted for an average of 83% of total amplified community DNA. At 16S rDNA similarities <98%, most taxa displayed differences in ARISA fragment lengths >10 bp, readily detectable and suggesting ARISA resolution is near the 'species' level. Prochlorococcus abundance profiles from ARISA were strongly correlated (r2=0.86) to Prochlorococcus cell counts, indicating ARISA data are roughly proportional to actual cell abundance within a defined taxon. Analysis of ARISA profiles for 42 months elucidated patterns of microbial presence and abundance providing insights into community shifts and ecological niches for specific organisms, including a coupling of ecological patterns for taxa within the Prochlorococcus, the Gamma Proteobacteria and Actinobacteria. Clade-specific ARISA protocols were developed for the SAR11 and marine cyanobacteria to resolve ambiguous identifications and to perform focused studies. 16S-ITS data allowed high-resolution identification of organisms by ITS sequence analysis, and examination of microdiversity.     

Effect of the tropical forage calliandra on microbial protein synthesis and ecology in the rumen

C.S. MCSWEENEY, B. PALMER, R. BUNCH AND D.O. KRAUSE. 2001 .
Aims: To determine the effect of condensed tannins in Calliandra calothyrsus (calliandra) on rumen microbial function.
Methods and Results: Microbial populations, ruminal protein synthesis and fermentation end-products were measured in sheep fed roughage hay supplemented with calliandra (30%), with and without inclusions of polyethylene glycol (PEG) to counteract the effect of tannin. Molecular and conventional enumeration techniques were used to quantify rumen bacteria, fungi and protozoa, and protein synthesis was predicted from estimates of urinary purine excretion. The total number of cellulolytic bacteria, including populations of Fibrobacter succinogenes and Ruminococcus spp., was significantly lower in sheep supplemented with calliandra and these populations increased when animals were treated with PEG. By contrast, protozoa and fungi and the microbial group containing Bacteroides - Porphyromonas - Prevotella bacteria appeared to be less affected. The efficiency of microbial protein synthesis in the rumen was not altered significantly.
Conclusions: Calliandra caused significant shifts in rumen microbial populations without changing the efficiency of protein synthesis.
Significance and Impact of the Study: The effect of calliandra tannins on rumen digestion may result more from complexing with nutrients than direct inhibition of micro-organisms.     

Distribution of High Bacterial Taxa Across the Chronosequence of Two Alpine Glacier Forelands

Little is known about the changes in abundance of microbial taxa in relation to the chronosequence of receding glaciers. This study investigated how the abundances of ten bacterial phyla or classes varied along successional gradients in two glaciers, Ödenwinkelkees and Rotmoosferner, in the central Alps. Quantitative PCR was used to estimate the abundance of the different bacterial taxa in extended glacier chronosequences, including 10- to 160-year-old successional stages, the surface of the glacier, and a fully established soil. Actinobacteria (15–30%) was the dominant group within the chronosequences. Several taxa showed significant differences in the number of taxa-specific 16S rRNA gene copies per nanogram of DNA and/or in the ratio of taxa-specific to the total bacterial 16S rRNA gene copies (i.e., the relative abundance of the different taxa within the bacterial community) between the established soils or the glacier surface and the 10- to 160-year-old successional stages. A significantly higher proportion of Βetaproteobacteria (20%) was observed on the surface of both glaciers. However, no differences were observed between the 10- to 160-year-old successional stages in the number of taxa-specific 16S rRNA gene copies per nanogram of DNA or in the ratio of taxa-specific to the total bacterial 16S rRNA gene copies for the different taxa. Nevertheless, when the relative abundance data from all the studied taxa were combined and analyzed altogether, most of the sites could be distinguished from one other. This indicates that the overall composition of the bacterial community was more affected than the abundance of the targeted taxa by changes in environmental conditions along the chronosequences.     

Evaluation of group-specific, 16S rRNA-targeted scissor probes for quantitative detection of predominant bacterial populations in dairy cattle rumen

AIMS: To develop a suite of group-specific, rRNA-targeted oligonucleotide scissor probes for the quantitative detection of the predominant bacterial groups within the ruminal microbial community with the rRNA cleavage reaction-mediated microbial quantification method. METHODS AND RESULTS: Oligonucleotides that complement the conserved sites of the 16S rRNA of phylogenetically defined groups of bacteria that significantly contribute to the anaerobic fermentation of carbohydrates in ruminal ecosystems were selected from among published probes or were newly designed. For each probe, target-specific rRNA cleavage was achieved by optimizing the formamide concentration in the reaction mixture. The set of scissor probes was then used to analyse the bacterial community in the rumen fluids of four healthy dairy cows. In the rumen fluid samples, the genera Bacteroides/Prevotella and Fibrobacter and the Clostridium coccoides-Eubacterium rectale group were detected in abundance, accounting for 44-48%, 2.9-10%, and 9.1-10% of the total 16S rRNA, respectively. The coverage with the probe set was 71-78% of the total bacterial 16S rRNA. CONCLUSIONS: The probe set coupled with the sequence-specific small-subunit rRNA cleavage method can be used to analyse the structure of a ruminal bacterial community. SIGNIFICANCE AND IMPACT OF THE STUDY: The probe set developed in this study provides a tool for comprehensive rRNA-based monitoring of the community members that dominate ruminal ecosystems. As the ruminal microbial community can be perturbed, it is important to track its dynamics by analysing microbiological profiles under specific conditions. The method described here will provide a convenient approach for such tracking.