Sandra Kittelmann

Simultaneous Amplicon Sequencing to Explore Co-Occurrence Patterns of Bacterial, Archaeal and Eukaryotic Microorganisms in Rumen Microbial Communities

Ruminants rely on a complex rumen microbial community to convert dietary plant material to energy-yielding products. Here we developed a method to simultaneously analyze the communitys bacterial and archaeal 16S rRNA genes, ciliate 18S rRNA genes and anaerobic fungal internal transcribed spacer 1 genes using 12 DNA samples derived from 11 different rumen samples from three host species (Ovis aries, Bos taurus, Cervus elephas) and multiplex 454 Titanium pyrosequencing. We show that the mixing ratio of the group-specific DNA templates before emulsion PCR is crucial to compensate for differences in amplicon length. This method, in contrast to using a non-specific universal primer pair, avoids sequencing non-targeted DNA, such as plant- or endophyte-derived rRNA genes, and allows increased or decreased levels of community structure resolution for each microbial group as needed. Communities analyzed with different primers always grouped by sample origin rather than by the primers used. However, primer choice had a greater impact on apparent archaeal community structure than on bacterial community structure, and biases for certain methanogen groups were detected. Co-occurrence analysis of microbial taxa from all three domains of life suggested strong within- and between-domain correlations between different groups of microorganisms within the rumen. The approach used to simultaneously characterize bacterial, archaeal and eukaryotic components of a microbiota should be applicable to other communities occupying diverse habitats.

Effect of DNA extraction methods and sampling techniques on the apparent structure of cow and sheep rumen microbial communities.

Molecular microbial ecology techniques are widely used to study the composition of the rumen microbiota and to increase understanding of the roles they play. Therefore, sampling and DNA extraction methods that result in adequate yields of microbial DNA that also accurately represents the microbial community are crucial. Fifteen different methods were used to extract DNA from cow and sheep rumen samples. The DNA yield and quality, and its suitability for downstream PCR amplifications varied considerably, depending on the DNA extraction method used. DNA extracts from nine extraction methods that passed these first quality criteria were evaluated further by quantitative PCR enumeration of microbial marker loci. Absolute microbial numbers, determined on the same rumen samples, differed by more than 100-fold, depending on the DNA extraction method used. The apparent compositions of the archaeal, bacterial, ciliate protozoal, and fungal communities in identical rumen samples were assessed using 454 Titanium pyrosequencing. Significant differences in microbial community composition were observed between extraction methods, for example in the relative abundances of members of the phyla Bacteroidetes and Firmicutes. Microbial communities in parallel samples collected from cows by oral stomach-tubing or through a rumen fistula, and in liquid and solid rumen digesta fractions, were compared using one of the DNA extraction methods. Community representations were generally similar, regardless of the rumen sampling technique used, but significant differences in the abundances of some microbial taxa such as the Clostridiales and the Methanobrevibacter ruminantium clade were observed. The apparent microbial community composition differed between rumen sample fractions, and Prevotellaceae were most abundant in the liquid fraction. DNA extraction methods that involved phenol-chloroform extraction and mechanical lysis steps tended to be more comparable. However, comparison of data from studies in which different sampling techniques, different rumen sample fractions or different DNA extraction methods were used should be avoided.

Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH4 mitigation at the levels of microbiota composition and transcriptional regulation.

Two different bacterial community types are linked with the low-methane emission trait in sheep.

The potent greenhouse gas methane (CH4) is produced in the rumens of ruminant animals from hydrogen produced during microbial degradation of ingested feed. The natural animal-to-animal variation in the amount of CH4 emitted and the heritability of this trait offer a means for reducing CH4 emissions by selecting low-CH4 emitting animals for breeding. We demonstrate that differences in rumen microbial community structure are linked to high and low CH4 emissions in sheep. Bacterial community structures in 236 rumen samples from 118 high- and low-CH4 emitting sheep formed gradual transitions between three ruminotypes. Two of these (Q and S) were linked to significantly lower CH4 yields (14.4 and 13.6 g CH4/kg dry matter intake [DMI], respectively) than the third type (H; 15.9 g CH4/kg DMI; p<0.001). Low-CH4 ruminotype Q was associated with a significantly lower ruminal acetate to propionate ratio (3.7±0.4) than S (4.4±0.7; p<0.001) and H (4.3±0.5; p<0.001), and harbored high relative abundances of the propionate-producing Quinella ovalis. Low-CH4 ruminotype S was characterized by lactate- and succinate-producing Fibrobacter spp., Kandleria vitulina, Olsenella spp., Prevotella bryantii, and Sharpea azabuensis. High-CH4 ruminotype H had higher relative abundances of species belonging to Ruminococcus, other Ruminococcaceae, Lachnospiraceae, Catabacteriaceae, Coprococcus, other Clostridiales, Prevotella, other Bacteroidales, and Alphaproteobacteria, many of which are known to form significant amounts of hydrogen. We hypothesize that lower CH4 yields are the result of bacterial communities that ferment ingested feed to relatively less hydrogen, which results in less CH4 being formed.

Characterization of rumen ciliate community composition in domestic sheep, deer, and cattle, feeding on varying diets, by means of PCR-DGGE and clone libraries

The structure and variability of ciliate protozoal communities in the rumens of domestic New Zealand ruminants feeding on different diets was investigated. The relative abundance of ciliates compared with bacteria was similar across all samples. However, molecular fingerprinting of communities showed ruminant-specific differences in species composition. Community compositions of cattle were significantly influenced by diet. In contrast, diet effects in deer and sheep were weaker than the animal-to-animal variation. Cloning and sequencing of almost-full-length 18S rRNA genes from representative samples revealed that New Zealand ruminants were colonized by at least nine genera of ciliates and allowed the assignment of samples to two distinct community types. Cattle contained A-type communities, with most sequences closely related to those of the genera Polyplastron and Ostracodinium. Deer and sheep (with one exception) harboured B-type communities, with the majority of sequences belonging to the genera Epidinium and Eudiplodinium. It has been suggested that species composition of ciliate communities may impact methane formation in ruminants, with the B-type producing more methane. Therefore, manipulation of ciliate communities may be a means of mitigating methane emissions from grazing sheep and deer in New Zealand.

Nitrogen metabolism and rumen microbial enumeration in lactating cows with divergent residual feed intake fed high-digestibility pasture

Dairy cattle selected for negative residual feed intake (n-RFI; efficient) should maintain production while reducing dry matter intake over a lactation because of improvements in feed digestion and efficient use of nutrients. The objective of this study was to measure nitrogen (N) digestibility and rumen microbial community composition over a short period during early lactation in lactating Holstein-Friesian cows selected previously for divergent RFI. It was proposed that n-RFI cows would have greater apparent digestibility of N than the positive RFI (p-RFI; inefficient) animals, to compensate for the lower dry matter intake determined during selection for divergence. Sixteen 3-yr-old rumen-cannulated, lactating cows (56 ± 10d in milk) selected for n-RFI (n = 8) and p-RFI (n = 8) were housed in metabolism stalls and fed fresh vegetative ryegrass (Lolium perenne L.) pasture ad libitum as a sole diet during an 8-d digestibility study. Intake of nutrients and outputs of milk, feces, and urine were determined. Rumen parameters were determined by removing, weighing, and sampling digesta, and by cobalt-EDTA dilution. Intakes of N, dry matter, organic matter, or its components did not differ with RFI. Compared with p-RFI cows, n-RFI cows had a greater apparent N digestibility (77.2 vs. 75.5%), and a tendency toward greater dry matter and organic matter digestibilities. The n-RFI cows had a lower fecal N output (126 vs. 138 g/d) and a lower partition of feed N to fecal N (23.1 vs. 24.7%) compared with p-RFI animals. We found no differences between phenotypes in the partition of N to urinary N or milk crude protein but did observe a trend for n-RFI cows to partition less N to milk casein (16.8 vs. 17.9%). Rumen digesta mass was similar for both groups, despite differences in calculated fractional liquid outflow rates, and most bacterial, archaeal, protozoal, and fungal communities were similar for both phenotype groups. In conclusion, dry matter intake and rumen function were similar for both phenotypes when the animals were fed highly digestible fresh ryegrass, but apparent digestibility of dietary N was higher in the efficient (n-RFI) cows. Future research should measure digestion parameters in cows with divergent RFI when fed diets differing in chemical composition (e.g., divergent crude protein contents).

RIM-DB: a taxonomic framework for community structure analysis of methanogenic archaea from the rumen and other intestinal environments

Methane is formed by methanogenic archaea in the rumen as one of the end products of feed fermentation in the ruminant digestive tract. To develop strategies to mitigate anthropogenic methane emissions due to ruminant farming, and to understand rumen microbial differences in animal feed conversion efficiency, it is essential that methanogens can be identified and taxonomically classified with high accuracy. Currently available taxonomic frameworks offer only limited resolution beyond the genus level for taxonomic assignments of sequence data stemming from high throughput sequencing technologies. Therefore, we have developed a QIIME-compatible database (DB) designed for species-level taxonomic assignment of 16S rRNA gene amplicon data targeting methanogenic archaea from the rumen, and from animal and human intestinal tracts. Called RIM-DB (Rumen and Intestinal Methanogen-DB), it contains a set of 2,379 almost full-length chimera-checked 16S rRNA gene sequences, including 20 previously unpublished sequences from isolates from three different orders. The taxonomy encompasses the recently-proposed seventh order of methanogens, the Methanomassiliicoccales, and allows differentiation between defined groups within this order. Sequence reads from rumen contents from a range of ruminant-diet combinations were taxonomically assigned using RIM-DB, Greengenes and SILVA. This comparison clearly showed that taxonomic assignments with RIM-DB resulted in the most detailed assignment, and only RIM-DB taxonomic assignments allowed methanogens to be distinguished taxonomically at the species level. RIM-DB complements the use of comprehensive databases such as Greengenes and SILVA for community structure analysis of methanogens from the rumen and other intestinal environments, and allows identification of target species for methane mitigation strategies.

A proposed taxonomy of anaerobic fungi (class neocallimastigomycetes) suitable for large-scale sequence-based community structure analysis.

Anaerobic fungi are key players in the breakdown of fibrous plant material in the rumen, but not much is known about the composition and stability of fungal communities in ruminants. We analyzed anaerobic fungi in 53 rumen samples from farmed sheep (4 different flocks), cattle, and deer feeding on a variety of diets. Denaturing gradient gel electrophoresis fingerprinting of the internal transcribed spacer 1 (ITS1) region of the rrn operon revealed a high diversity of anaerobic fungal phylotypes across all samples. Clone libraries of the ITS1 region were constructed from DNA from 11 rumen samples that had distinctly different fungal communities. A total of 417 new sequences were generated to expand the number and diversity of ITS1 sequences available. Major phylogenetic groups of anaerobic fungi in New Zealand ruminants belonged to the genera Piromyces, Neocallimastix, Caecomyces and Orpinomyces. In addition, sequences forming four novel clades were obtained, which may represent so far undetected genera or species of anaerobic fungi. We propose a revised phylogeny and pragmatic taxonomy for anaerobic fungi, which was tested and proved suitable for analysis of datasets stemming from high-throughput next-generation sequencing methods. Comparing our revised taxonomy to the taxonomic assignment of sequences deposited in the GenBank database, we believe that >29% of ITS1 sequences derived from anaerobic fungal isolates or clones are misnamed at the genus level.

Isolation of previously uncultured rumen bacteria by dilution to extinction using a new liquid culture medium

A new anaerobic medium that mimics the salts composition of rumen fluid was used in conjunction with a dilution method of liquid culture to isolate fermentative bacteria from the rumen of a grass-fed sheep. The aim was to inoculate a large number of culture tubes each with a mean of <1 culturable cell, which should maximize the number of cultures that develop from a single bacterium. This minimizes the effort that has to be put into purifying the resultant cultures. Of 1000 tubes, 139 were growth positive. Of the 93 that were able to be subcultured, 54 (58%) appeared to be pure cultures. The phylogenetic placements of these 54 cultures, together with another 6 cultures obtained from a preliminary study, were determined. Using a criterion of <93% 16S rRNA gene sequence identity to a previously named bacterium as a proxy for defining a new genus, 27 (45%) of the 60 cultures belonged to 14 potentially novel genera. Many of these had 16S rRNA genes that shared >97% sequence identity to genes of uncultured bacteria detected in various gastrointestinal environments. This strategy has therefore allowed us to cultivate many novel rumen bacteria, opening the way to overcoming the lack of cultures of many of the groups detected using cultivation-independent methods.

Faecal Microbiota of Forage-Fed Horses in New Zealand and the Population Dynamics of Microbial Communities following Dietary Change

The effects of abrupt dietary transition on the faecal microbiota of forage-fed horses over a 3-week period were investigated. Yearling Thoroughbred fillies reared as a cohort were exclusively fed on either an ensiled conserved forage-grain diet (“Group A”; n = 6) or pasture (“Group B”; n = 6) for three weeks prior to the study. After the Day 0 faecal samples were collected, horses of Group A were abruptly transitioned to pasture. Both groups continued to graze similar pasture for three weeks, with faecal samples collected at 4-day intervals. DNA was isolated from the faeces and microbial 16S and 18S rRNA gene amplicons were generated and analysed by pyrosequencing. The faecal bacterial communities of both groups of horses were highly diverse (Simpson’s index of diversity >0.8), with differences between the two groups on Day 0 (P<0.017 adjusted for multiple comparisons). There were differences between Groups A and B in the relative abundances of four genera, BF311 (family Bacteroidaceae; P = 0.003), CF231 (family Paraprevotellaceae; P = 0.004), and currently unclassified members within the order Clostridiales (P = 0.003) and within the family Lachnospiraceae (P = 0.006). The bacterial community of Group A horses became similar to Group B within four days of feeding on pasture, whereas the structure of the archaeal community remained constant pre- and post-dietary change. The community structure of the faecal microbiota (bacteria, archaea and ciliate protozoa) of pasture-fed horses was also identified. The initial differences observed appeared to be linked to recent dietary history, with the bacterial community of the forage-fed horses responding rapidly to abrupt dietary change.