Ninety-nine de novo assembled genomes from the moose (Alces alces) rumen microbiome provide new insights into microbial plant biomass degradation

Multi-omic Directed Discovery of Cellulosomes, Polysaccharide Utilization Loci, and Lignocellulases from an Enriched Rumen Anaerobic Consortium

The lignocellulolytic ERAC displays a unique set of plant polysaccharide-degrading enzymes (with multimodular characteristics), cellulosomal complexes, and PULs. The MAGs described here represent an expansion of the genetic content of rumen bacterial genomes dedicated to plant polysaccharide degradation, therefore providing a valuable resource for the development of biocatalytic toolbox strategies to be applied to lignocellulose-based biorefineries.

Lignocellulose is one of the most abundant renewable carbon sources, representing an alternative to petroleum for the production of fuel and chemicals. Nonetheless, the lignocellulose saccharification process, to release sugars for downstream applications, is one of the most crucial factors economically challenging to its use. The synergism required among the various carbohydrate-active enzymes (CAZymes) for efficient lignocellulose breakdown is often not satisfactorily achieved with an enzyme mixture from a single strain. To overcome this challenge, enrichment strategies can be applied to develop microbial communities with an efficient CAZyme arsenal, incorporating complementary and synergistic properties, to improve lignocellulose deconstruction. We report a comprehensive and deep analysis of an enriched rumen anaerobic consortium (ERAC) established on sugarcane bagasse (SB). The lignocellulolytic abilities of the ERAC were confirmed by analyzing the depolymerization of bagasse by scanning electron microscopy, enzymatic assays, and mass spectrometry. Taxonomic analysis based on 16S rRNA sequencing elucidated the community enrichment process, which was marked by a higher abundance of Firmicutes and Synergistetes species. Shotgun metagenomic sequencing of the ERAC disclosed 41 metagenome-assembled genomes (MAGs) harboring cellulosomes and polysaccharide utilization loci (PULs), along with a high diversity of CAZymes. The amino acid sequences of the majority of the predicted CAZymes (60% of the total) shared less than 90% identity with the sequences found in public databases. Additionally, a clostridial MAG identified in this study produced proteins during consortium development with scaffoldin domains and CAZymes appended to dockerin modules, thus representing a novel cellulosome-producing microorganism.

IMPORTANCE The lignocellulolytic ERAC displays a unique set of plant polysaccharide-degrading enzymes (with multimodular characteristics), cellulosomal complexes, and PULs. The MAGs described here represent an expansion of the genetic content of rumen bacterial genomes dedicated to plant polysaccharide degradation, therefore providing a valuable resource for the development of biocatalytic toolbox strategies to be applied to lignocellulose-based biorefineries.

The effect of a diet based on rice straw co-fermented with probiotics and enzymes versus a fresh corn Stover-based diet on the rumen bacterial community and metabolites of beef cattle

Improvement of the food value of rice straw is urgently required in rice crop growing areas to mitigate pollution caused by rice straw burning and enhance the supply of high-quality forages for ruminants. The aims of the present study were to compare the effects of fresh corn Stover and rice straw co-fermented with probiotics and enzymes on rumen fermentation and establish the feasibility of increasing the rice straw content in ruminant diets and, by extension, reducing air pollution caused by burning rice straw. Twenty Simmental hybrid beef cattle were randomly allotted to two groups with ten cattle per group. They were fed diets based either on rice straw co-fermented with probiotics and enzymes or fresh corn Stover for 90 days. Rumen fluid was sampled with an esophageal tube vacuum pump device from each animal on the mornings of days 30, 60, and 90. Bacterial diversity was evaluated by sequencing the V4–V5 region of the 16S rRNA gene. Metabolomes were analyzed by gas chromatography/time-of-flight mass spectrometry (GC–TOF/MS). Compared to cattle fed fresh corn Stover, those fed rice straw co-fermented with probiotics and enzymes had higher (P < 0.05) levels of acetic acid and propionate in rumen liquid at d 60 and d 90 respectively, higher (P < 0.05) abundances of the phyla Bacteroidetes and Fibrobacteres and the genera Ruminococcus, Saccharofermentans, Pseudobutyrivibrio, Treponema, Lachnoclostridium, and Ruminobacter, and higher (P < 0.05) concentrations of metabolites involved in metabolisms of amino acid, carbohydrate, and cofactors and vitamins. Relative to fresh corn Stover, rice straw co-fermented with probiotics and enzymes resulted in higher VFA concentrations, numbers of complex carbohydrate-decomposing and H₂-utilizing bacteria, and feed energy conversion efficiency in the rumen.

A catalog of microbial genes from the bovine rumen unveils a specialized and diverse biomass-degrading environment

Background

The rumen microbiota provides essential services to its host and, through its role in ruminant production, contributes to human nutrition and food security. A thorough knowledge of the genetic potential of rumen microbes will provide opportunities for improving the sustainability of ruminant production systems. The availability of gene reference catalogs from gut microbiomes has advanced the understanding of the role of the microbiota in health and disease in humans and other mammals. In this work, we established a catalog of reference prokaryote genes from the bovine rumen.

Results

Using deep metagenome sequencing we identified 13,825,880 non-redundant prokaryote genes from the bovine rumen. Compared to human, pig, and mouse gut metagenome catalogs, the rumen is larger and richer in functions and microbial species associated with the degradation of plant cell wall material and production of methane. Genes encoding enzymes catalyzing the breakdown of plant polysaccharides showed a particularly high richness that is otherwise impossible to infer from available genomes or shallow metagenomics sequencing. The catalog expands the dataset of carbohydrate-degrading enzymes described in the rumen. Using an independent dataset from a group of 77 cattle fed 4 common dietary regimes, we found that only <0.1% of genes were shared by all animals, which contrast with a large overlap for functions, i.e., 63% for KEGG functions. Different diets induced differences in the relative abundance rather than the presence or absence of genes, which explains the great adaptability of cattle to rapidly adjust to dietary changes.

Conclusions

These data bring new insights into functions, carbohydrate-degrading enzymes, and microbes of the rumen to complement the available information on microbial genomes. The catalog is a significant biological resource enabling deeper understanding of phenotypes and biological processes and will be expanded as new data are made available.

Identification of Complex Rumen Microbiome Interaction Within Diverse Functional Niches as Mechanisms Affecting the Variation of Methane Emissions in Bovine

A network analysis including relative abundances of all ruminal microbial genera (archaea, bacteria, fungi, and protists) and their genes was performed to improve our understanding of how the interactions within the ruminal microbiome affects methane emissions (CH4). Metagenomics and CH4 data were available from 63 bovines of a two-breed rotational cross, offered two basal diets. Co-abundance network analysis revealed 10 clusters of functional niches. The most abundant hydrogenotrophic Methanobacteriales with key microbial genes involved in methanogenesis occupied a different functional niche (i.e., "methanogenesis" cluster) than methylotrophic Methanomassiliicoccales (Candidatus Methanomethylophylus) and acetogens (Blautia). Fungi and protists clustered together and other plant fiber degraders like Fibrobacter occupied a seperate cluster. A Partial Least Squares analysis approach to predict CH4 variation in each cluster showed the methanogenesis cluster had the best prediction ability (57.3%). However, the most important explanatory variables in this cluster were genes involved in complex carbohydrate degradation, metabolism of sugars and amino acids and Candidatus Azobacteroides carrying nitrogen fixation genes, but not methanogenic archaea and their genes. The cluster containing Fibrobacter, isolated from other microorganisms, was positively associated with CH4 and explained 49.8% of its variability, showing fermentative advantages compared to other bacteria and fungi in providing substrates (e.g., formate) for methanogenesis. In other clusters, genes with enhancing effect on CH4 were related to lactate and butyrate (Butyrivibrio and Pseudobutyrivibrio) production and simple amino acids metabolism. In comparison, ruminal genes negatively related to CH4 were involved in carbohydrate degradation via lactate and succinate and synthesis of more complex amino acids by γ-Proteobacteria. When analyzing low- and high-methane emitters data in separate networks, competition between methanogens in the methanogenesis cluster was uncovered by a broader diversity of methanogens involved in the three methanogenesis pathways and larger interactions within and between communities in low compared to high emitters. Generally, our results suggest that differences in CH4 are mainly explained by other microbial communities and their activities rather than being only methanogens-driven. Our study provides insight into the interactions of the rumen microbial communities and their genes by uncovering functional niches affecting CH4, which will benefit the development of efficient CH4 mitigation strategies.

Spatially remote motifs cooperatively affect substrate preference of a ruminal GH26-type endo-β-1,4-mannanase

β-Mannanases from the glycoside hydrolase 26 (GH26) family are retaining hydrolases that are active on complex heteromannans and whose genes are abundant in rumen metagenomes and metatranscriptomes. These enzymes can exhibit distinct modes of substrate recognition and are often fused to carbohydrate-binding modules (CBMs), resulting in a molecular puzzle of mechanisms governing substrate preference and mode of action that has not yet been pieced together. In this study, we recovered a novel GH26 enzyme with a CBM35 module linked to its N terminus (CrMan26) from a cattle rumen metatranscriptome. CrMan26 exhibited a preference for galactomannan as substrate and the crystal structure of the full-length protein at 1.85 Å resolution revealed a unique orientation of the ancillary domain relative to the catalytic interface, strategically positioning a surface aromatic cluster of the ancillary domain as an extension of the substrate-binding cleft, contributing to galactomannan preference. Moreover, systematic investigation of nonconserved residues in the catalytic interface unveiled that residues Tyr195 (-3 subsite) and Trp234 (-5 subsite) from distal negative subsites have a key role in galactomannan preference. These results indicate a novel and complex mechanism for substrate recognition involving spatially remote motifs, distal negative subsites from the catalytic domain, and a surface-associated aromatic cluster from the ancillary domain. These findings expand our molecular understanding of the mechanisms of substrate binding and recognition in the GH26 family and shed light on how some CBMs and their respective orientation can contribute to substrate preference.

Phylogenomic analysis of 589 metagenome-assembled genomes encompassing all major prokaryotic lineages from the gut of higher termites

"Higher" termites have been able to colonize all tropical and subtropical regions because of their ability to digest lignocellulose with the aid of their prokaryotic gut microbiota. Over the last decade, numerous studies based on 16S rRNA gene amplicon libraries have largely described both the taxonomy and structure of the prokaryotic communities associated with termite guts. Host diet and microenvironmental conditions have emerged as the main factors structuring the microbial assemblages in the different gut compartments. Additionally, these molecular inventories have revealed the existence of termite-specific clusters that indicate coevolutionary processes in numerous prokaryotic lineages. However, for lack of representative isolates, the functional role of most lineages remains unclear. We reconstructed 589 metagenome-assembled genomes (MAGs) from the different gut compartments of eight higher termite species that encompass 17 prokaryotic phyla. By iteratively building genome trees for each clade, we significantly improved the initial automated assignment, frequently up to the genus level. We recovered MAGs from most of the termite-specific clusters in the radiation of, for example, Planctomycetes, Fibrobacteres, Bacteroidetes, Euryarchaeota, Bathyarchaeota, Spirochaetes, Saccharibacteria, and Firmicutes, which to date contained only few or no representative genomes. Moreover, the MAGs included abundant members of the termite gut microbiota. This dataset represents the largest genomic resource for arthropod-associated microorganisms available to date and contributes substantially to populating the tree of life. More importantly, it provides a backbone for studying the metabolic potential of the termite gut microbiota, including the key members involved in carbon and nitrogen biogeochemical cycles, and important clues that may help cultivating representatives of these understudied clades.

Fibre digestion by rumen microbiota — a review of recent metagenomic and metatranscriptomic studies

Plant biomass is the most abundant renewable resource on the planet, and the biopolymers of lignocellulose are the foundation of ruminant production systems. Optimizing the saccharification of lignocellulosic feeds is a crucial step in their bioconversion to ruminant protein. Plant cell walls are chemically heterogeneous structures that have evolved to provide structural support and protection to the plant. Ruminants are the most efficient digesters of lignocellulose due to a rich array of bacteria, archaea, fungi, and protozoa within the rumen and lower digestive tract. Metagenomic and metatranscriptomic studies have enhanced the current understanding of the composition, diversity, and function of the rumen microbiome. There is particular interest in identifying the carbohydrate-active enzymes responsible for the ruminal degradation of plant biomass. Understanding the roles of cellulosomes- and polysaccharide-utilising loci in ruminal fibre degradation could provide insight into strategies to enhance forage utilisation by ruminants. Despite advancements in “omics” technology, the majority of rumen microorganisms are still uncharacterised, and their ability to act synergistically is still not understood. By advancing our current knowledge of rumen fibre digestion, there may be opportunity to further improve the productive performance of ruminants fed forage diets.

Top-Down Enrichment Guides in Formation of Synthetic Microbial Consortia for Biomass Degradation

Consortium-based approaches are a promising avenue toward efficient bioprocessing. However, many complex microbial interactions dictate community dynamics and stability that must be replicated in synthetic systems. The rumen and/or hindguts of large mammalian herbivores harbor complex communities of biomass-degrading fungi and bacteria, as well as archaea and protozoa that work collectively to degrade lignocellulose, yet the microbial interactions responsible for stability, resilience, and activity of the community remain largely uncharacterized. In this work, we demonstrate a "top-down" enrichment-based methodology for selecting a minimal but effective lignocellulose-degrading community that produces methane-rich fermentation gas (biogas). The resulting enrichment consortium produced 0.75-1.9-fold more fermentation gas at 1.4-2.1 times the rate compared to a monoculture of fungi from the enrichment. Metagenomic sequencing of the top-down enriched consortium revealed genomes encoding for functional compartmentalization of the community, spread across an anaerobic fungus (Piromyces), a bacterium (Sphaerochaeta), and two methanogenic archaea (Methanosphaera and Methanocorpusculum). Guided by the composition of the top-down enrichment, several synthetic cocultures were formed from the "bottom-up" using previously isolated fungi, Neocallimastix californiae and Anaeromyces robustus paired with the methanogen Methanobacterium bryantii. While cross-feeding occurred in synthetic co-cultures, removal of fungal metabolites by methanogens did not increase the rate of gas production or the rate of substrate deconstruction by the synthetic community relative to fungal monocultures. Metabolomic characterization verified that syntrophy was established within synthetic co-cultures, which generated methane at similar concentrations compared to the enriched consortium but lacked the temporal stability (resilience) seen in the native system. Taken together, deciphering the membership and metabolic potential of an enriched gut consortium enables the design of methanogenic synthetic co-cultures. However, differences in the growth rate and stability of enriched versus synthetic consortia underscore the difficulties in mimicking naturally occurring syntrophy in synthetic systems.

Compendium of 4,941 rumen metagenome-assembled genomes for rumen microbiome biology and enzyme discovery

Ruminants provide essential nutrition for billions of people worldwide. The rumen is a specialized stomach that is adapted to the breakdown of plant-derived complex polysaccharides. The genomes of the rumen microbiota encode thousands of enzymes adapted to digestion of the plant matter that dominates the ruminant diet. We assembled 4,941 rumen microbial metagenome-assembled genomes (MAGs) using approximately 6.5 terabases of short- and long-read sequence data from 283 ruminant cattle. We present a genome-resolved metagenomics workflow that enabled assembly of bacterial and archaeal genomes that were at least 80% complete. Of note, we obtained three single-contig, whole-chromosome assemblies of rumen bacteria, two of which represent previously unknown rumen species, assembled from long-read data. Using our rumen genome collection we predicted and annotated a large set of rumen proteins. Our set of rumen MAGs increases the rate of mapping of rumen metagenomic sequencing reads from 15% to 50–70%. These genomic and protein resources will enable a better understanding of the structure and functions of the rumen microbiota.

Comprehensive resource of cattle rumen genomes and a database of predicted proteins.

Carbohydrate Hydrolytic Potential and Redundancy of an Anaerobic Digestion Microbiome Exposed to Acidosis, as Uncovered by Metagenomics

The enzymatic hydrolysis of lignocellulosic biomass is mainly driven by the action of carbohydrate-active enzymes. By characterizing the gene profiles at the different stages of the anaerobic digestion experiment, we showed that the microbiome retains its hydrolytic functional redundancy even during severe acidosis, despite significant changes in taxonomic composition. By analyzing reconstructed bacterial genomes, we demonstrate that Bacteroidetes hydrolytic gene diversity likely favors the abundance of this phylum in some anaerobic digestion systems. Further, we observe genetic redundancy within the Bacteroidetes group, which accounts for the preserved hydrolytic potential during acidosis. This work also uncovers new polysaccharide utilization loci involved in the deconstruction of various biomasses and proposes the model of acetylated glucomannan degradation by Bacteroidetes. Acetylated glucomannan-enriched biomass is a common substrate for many industries, including pulp and paper production. Using naturally evolved cocktails of enzymes for biomass pretreatment could be an interesting alternative to the commonly used chemical pretreatments.

Increased hydrolysis of easily digestible biomass may lead to acidosis of anaerobic reactors and decreased methane production. Previously, it was shown that the structure of microbial communities changed during acidosis; however, once the conditions are back to optimal, biogas (initially CO₂) production quickly restarts. This suggests the retention of the community functional redundancy during the process failure. In this study, with the use of metagenomics and downstream bioinformatics analyses, we characterize the carbohydrate hydrolytic potential of the microbial community, with a special focus on acidosis. To that purpose, carbohydrate-active enzymes were identified, and to further link the community hydrolytic potential with key microbes, bacterial genomes were reconstructed. In addition, we characterized biochemically the specificity and activity of selected enzymes, thus verifying the accuracy of the in silico predictions. The results confirm the retention of the community hydrolytic potential during acidosis and indicate Bacteroidetes to be largely involved in biomass degradation. Bacteroidetes showed higher diversity and genomic content of carbohydrate hydrolytic enzymes that might favor the dominance of this phylum over other bacteria in some anaerobic reactors. The combination of bioinformatic analyses and activity tests enabled us to propose a model of acetylated glucomannan degradation by Bacteroidetes.

IMPORTANCE The enzymatic hydrolysis of lignocellulosic biomass is mainly driven by the action of carbohydrate-active enzymes. By characterizing the gene profiles at the different stages of the anaerobic digestion experiment, we showed that the microbiome retains its hydrolytic functional redundancy even during severe acidosis, despite significant changes in taxonomic composition. By analyzing reconstructed bacterial genomes, we demonstrate that Bacteroidetes hydrolytic gene diversity likely favors the abundance of this phylum in some anaerobic digestion systems. Further, we observe genetic redundancy within the Bacteroidetes group, which accounts for the preserved hydrolytic potential during acidosis. This work also uncovers new polysaccharide utilization loci involved in the deconstruction of various biomasses and proposes the model of acetylated glucomannan degradation by Bacteroidetes. Acetylated glucomannan-enriched biomass is a common substrate for many industries, including pulp and paper production. Using naturally evolved cocktails of enzymes for biomass pretreatment could be an interesting alternative to the commonly used chemical pretreatments.

In marine Bacteroidetes the bulk of glycan degradation during algae blooms is mediated by few clades using a restricted set of genes

We investigated Bacteroidetes during spring algae blooms in the southern North Sea in 2010–2012 using a time series of 38 deeply sequenced metagenomes. Initial partitioning yielded 6455 bins, from which we extracted 3101 metagenome-assembled genomes (MAGs) including 1286 Bacteroidetes MAGs covering ~120 mostly uncultivated species. We identified 13 dominant, recurrent Bacteroidetes clades carrying a restricted set of conserved polysaccharide utilization loci (PULs) that likely mediate the bulk of bacteroidetal algal polysaccharide degradation. The majority of PULs were predicted to target the diatom storage polysaccharide laminarin, alpha-glucans, alpha-mannose-rich substrates, and sulfated xylans. Metaproteomics at 14 selected points in time revealed expression of SusC-like proteins from PULs targeting all of these substrates. Analyses of abundant key players and their PUL repertoires over time furthermore suggested that fewer and simpler polysaccharides dominated early bloom stages, and that more complex polysaccharides became available as blooms progressed.

Mining of camel rumen metagenome to identify novel alkali-thermostable xylanase capable of enhancing the recalcitrant lignocellulosic biomass conversion

The aim of this study was to isolate and characterize novel alkali-thermostable xylanase genes from the mixed genome DNA of camel rumen metagenome. In this study, a five-stage computational screening procedure was utilized to find the primary candidate enzyme with superior properties from the camel rumen metagenome. This enzyme was subjected to cloning, purification, and structural and functional characterization. It showed high thermal stability, high activity in a broad range of pH (6-11) and temperature (30-90 °C) and effectivity in recalcitrant lignocellulosic biomass degradation. Our results demonstrated the power of in silico analysis to discover novel alkali-thermostable xylanases, effective for the bioconversion of lignocellulosic biomass.

Invited review: Application of meta-omics to understand the dynamic nature of the rumen microbiome and how it responds to diet in ruminants

Ruminants are unique among livestock due to their ability to efficiently convert plant cell wall carbohydrates into meat and milk. This ability is a result of the evolution of an essential symbiotic association with a complex microbial community in the rumen that includes vast numbers of bacteria, methanogenic archaea, anaerobic fungi and protozoa. These microbes produce a diverse array of enzymes that convert ingested feedstuffs into volatile fatty acids and microbial protein which are used by the animal for growth. Recent advances in high-throughput sequencing and bioinformatic analyses have helped to reveal how the composition of the rumen microbiome varies significantly during the development of the ruminant host, and with changes in diet. These sequencing efforts are also beginning to explain how shifts in the microbiome affect feed efficiency. In this review, we provide an overview of how meta-omics technologies have been applied to understanding the rumen microbiome, and the impact that diet has on the rumen microbial community.

Molecular hydrogen produced by elemental magnesium inhibits rumen fermentation and enhances methanogenesis in dairy cows

Hydrogen is a key metabolite that connects microbial fermentation and methanogenesis in the rumen. This study was to investigate the effects of elevated H₂ produced by elemental Mg on rumen fermentation and methanogenesis in dairy cows. Four nonlactating Chinese Holstein dairy cows were employed for this experiment in a replicated crossover design. The 2 dietary treatments included a basal diet supplemented with Mg(OH)₂ (14.5 g/kg of feed dry matter) or elemental Mg (6.00 g/kg of feed dry matter). When compared with Mg(OH)₂ treatment, cows fed diet with elemental Mg had similar rumen Mg²⁺ concentration, but higher rumen dissolved H₂ and methane concentrations at 2.5 h after morning feeding. Also, elemental Mg supplementation decreased feed digestibility, rumen volatile fatty acid concentration, and relative abundance of group Ruminococcaceae_UCG-014, genus Bifidobacterium, and group Mollicutes_RF9, increased acetate to propionate ratio, succinate concentration, and abundance of family Christensenellaceae. Elemental Mg supplementation increased enteric CH₄ emission, altered methanogen community with increased abundance of order Methanomassiliicoccales, 16S ribosomal RNA gene copies of methanogens, and order Methanobacteriales. In summary, the pulse of elevated dissolved H₂ after feeding produced by elemental Mg inhibited rumen fermentation and feed digestibility by decreasing the abundance of carbohydrate-degrading bacteria, promoted H₂ incorporation into succinate by increasing family Christensenellaceae and genus Bacteroidales_BS11, and increased H₂ utilization for methanogenesis by favoring growth of methanogens.

Bacterial colonisation of reeds and cottonseed hulls in the rumen of Tarim red deer (Cervus elaphus yarkandensis)

The rumen microbiome contributes greatly to the degradation of plant fibres to volatile fatty acids and microbial products, affecting the health and productivity of ruminants. In this study, we investigated the dynamics of colonisation by bacterial communities attached to reeds and cottonseed hulls in the rumen of Tarim red deer, a native species distributed in the desert of the Tarim Basin. The reed and cottonseed hull samples incubated in nylon bags for 1, 6, 12, and 48 h were collected and used to examine the bacterial communities by next-generation sequencing of the bacterial 16S rRNA gene. Prevotella1 and Rikenellaceae RC9 were the most abundant taxa in both the reed and cottonseed hull groups at various times, indicating a key role of these organisms in rumen fermentation in Tarim red deer. The relative abundances of cellulolytic bacteria, such as members of Fibrobacter, Treponema 2, Ruminococcaceae NK4A214 and Succiniclasticum increased, while that of the genus Prevotella 1 decreased, with increasing incubation time in both reeds and cottonseed hulls. Moreover, the temporal changes in bacterial diversity between reeds and cottonseed hulls were different, as demonstrated by the variations in the taxa Ruminococcaceae UCG 010 and Papillibacter in the reed group and Sphaerochaeta and Erysipelotrichaceae UCG 004 in the cottonseed hull group; the abundances of these bacteria first decreased and then increased. In conclusion, our results reveal the dynamics of bacterial colonisation of reeds and cottonseed hulls in the rumen of Tarim red deer.

Metagenomic reconstructions of gut microbial metabolism in weanling pigs

BACKGROUND

The piglets' transition from milk to solid feed induces a succession of bacterial communities, enhancing the hosts' ability to harvest energy from dietary carbohydrates. To reconstruct microbial carbohydrate metabolism in weanling pigs, this study combined 16S rRNA gene sequencing (n = 191) and shotgun metagenomics (n = 72).

RESULTS

Time and wheat content in feed explained most of the variation of the microbiota as assessed by 16S rRNA gene sequencing in weanling pigs. De novo metagenomic binning reconstructed 360 high-quality genomes that represented 11 prokaryotic and 1 archaeal phylum. Analysis of carbohydrate metabolism in these genomes revealed that starch fermentation is carried out by a consortium of Firmicutes expressing extracellular α-(1 → 4)-glucan branching enzyme (GH13) and Bacteroidetes expressing periplasmic neopullulanase (GH13) and α-glucosidase (GH97). Fructans were degraded by extracellular GH32 enzymes from Bacteriodetes and Lactobacillus. Lactose fermentation by β-galactosidases (GH2 and GH42) was identified in Firmicutes. In conclusion, the assembly of 360 high-quality genomes as the first metagenomic reference for swine intestinal microbiota allowed identification of key microbial contributors to degradation of starch, fructans, and lactose.

CONCLUSIONS

Microbial consortia that are responsible for degradation of these glycans differ substantially from the microbial consortia that degrade the same glycans in humans. Our study thus enables improvement of feeding models with higher feed efficiency and better pathogen control for weanling pigs.

Functional metagenomics reveals abundant polysaccharide-degrading gene clusters and cellobiose utilization pathways within gut microbiota of a wood-feeding higher termite

Plant cell-wall polysaccharides constitute the most abundant but recalcitrant organic carbon source in nature. Microbes residing in the digestive tract of herbivorous bilaterians are particularly efficient at depolymerizing polysaccharides into fermentable sugars and play a significant support role towards their host's lifestyle. Here, we combine large-scale functional screening of fosmid libraries, shotgun sequencing, and biochemical assays to interrogate the gut microbiota of the wood-feeding "higher" termite Globitermes brachycerastes. A number of putative polysaccharide utilization gene clusters were identified with multiple fibrolytic genes. Our large-scale functional screening of 50,000 fosmid clones resulted in 464 clones demonstrating plant polysaccharide-degrading activities, including 267 endoglucanase-, 24 exoglucanase-, 72 β-glucosidase-, and 101 endoxylanase-positive clones. We sequenced 173 functionally active clones and identified ~219 genes encoding putative carbohydrate-active enzymes (CAZymes) targeting cellulose, hemicellulose and pectin. Further analyses revealed that 68 of 154 contigs encode one or more CAZyme, which includes 35 examples of putative saccharolytic operons, suggesting that clustering of CAZymes is common in termite gut microbial inhabitants. Biochemical characterization of a representative xylanase cluster demonstrated that constituent enzymes exhibited complementary physicochemical properties and saccharolytic capabilities. Furthermore, diverse cellobiose-metabolizing enzymes include β-glucosidases, cellobiose phosphorylases, and phopho-6-β-glucosidases were identified and functionally verified, indicating that the termite gut micro-ecosystem utilizes diverse metabolic pathways to interconnect hydrolysis and central metabolism. Collectively, these results provide an in-depth view of the adaptation and digestive strategies employed by gut microbiota within this tiny-yet-efficient host-associated ecosystem.

Fiber-associated spirochetes are major agents of hemicellulose degradation in the hindgut of wood-feeding higher termites

Symbiotic digestion of lignocellulose in wood-feeding higher termites (family Termitidae) is a two-step process that involves endogenous host cellulases secreted in the midgut and a dense bacterial community in the hindgut compartment. The genomes of the bacterial gut microbiota encode diverse cellulolytic and hemicellulolytic enzymes, but the contributions of host and bacterial symbionts to lignocellulose degradation remain ambiguous. Our previous studies of Nasutitermes spp. documented that the wood fibers in the hindgut paunch are consistently colonized not only by uncultured members of Fibrobacteres, which have been implicated in cellulose degradation, but also by unique lineages of Spirochaetes. Here, we demonstrate that the degradation of xylan, the major component of hemicellulose, is restricted to the hindgut compartment, where it is preferentially hydrolyzed over cellulose. Metatranscriptomic analysis documented that the majority of glycoside hydrolase (GH) transcripts expressed by the fiber-associated bacterial community belong to family GH11, which consists exclusively of xylanases. The substrate specificity was further confirmed by heterologous expression of the gene encoding the predominant homolog. Although the most abundant transcripts of GH11 in Nasutitermes takasagoensis were phylogenetically placed among their homologs of Firmicutes, immunofluorescence microscopy, compositional binning of metagenomics contigs, and the genomic context of the homologs indicated that they are encoded by Spirochaetes and were most likely obtained by horizontal gene transfer among the intestinal microbiota. The major role of spirochetes in xylan degradation is unprecedented and assigns the fiber-associated Treponema clades in the hindgut of wood-feeding higher termites a prominent part in the breakdown of hemicelluloses.

Interspecies cross-feeding orchestrates carbon degradation in the rumen ecosystem

Because of their agricultural value, there is a great body of research dedicated to understanding the microorganisms responsible for rumen carbon degradation. However, we lack a holistic view of the microbial food web responsible for carbon processing in this ecosystem. Here, we sampled rumen-fistulated moose, allowing access to rumen microbial communities actively degrading woody plant biomass in real time. We resolved 1,193 viral contigs and 77 unique, near-complete microbial metagenome-assembled genomes, many of which lacked previous metabolic insights. Plant-derived metabolites were measured with NMR and carbohydrate microarrays to quantify the carbon nutrient landscape. Network analyses directly linked measured metabolites to expressed proteins from these unique metagenome-assembled genomes, revealing a genome-resolved three-tiered carbohydrate-fuelled trophic system. This provided a glimpse into microbial specialization into functional guilds defined by specific metabolites. To validate our proteomic inferences, the catalytic activity of a polysaccharide utilization locus from a highly connected metabolic hub genome was confirmed using heterologous gene expression. Viral detected proteins and linkages to microbial hosts demonstrated that phage are active controllers of rumen ecosystem function. Our findings elucidate the microbial and viral members, as well as their metabolic interdependencies, that support in situ carbon degradation in the rumen ecosystem.

A combination of proteomics, metagenome-assembled genomes and heterologous gene expression experiments reveals a trophic system for carbon utilization in the moose rumen microbiome and provides insights into phage dynamics in this ecosystem.

Review: The application of omics to rumen microbiota function

Rumen microbiome profiling uses 16S rRNA (18S rRNA, internal transcribed spacer) gene sequencing, a method that usually sequences a small portion of a single gene and is often biased and varies between different laboratories. Functional information can be inferred from this data, but only for those that are closely related to known annotated species, and even then may not truly reflect the function performed within the environment being studied. Genome sequencing of isolates and metagenome-assembled genomes has now reached a stage where representation of the majority of rumen bacterial genera are covered, but this still only represents a portion of rumen microbial species. The creation of a microbial genome (bins) database with associated functional annotations will provide a consistent reference to allow mapping of RNA-Seq reads for functional gene analysis from within the rumen microbiome. The integration of multiple omic analytics is linking functional gene activity, metabolic pathways and rumen metabolites with the responsible microbiota, supporting our biological understanding of the rumen system. The application of these techniques has advanced our understanding of the major microbial populations and functional pathways that are used in relation to lower methane emissions, higher feed efficiencies and responses to different feeding regimes. Continued and more precise use of these tools will lead to a detailed and comprehensive understanding of compositional and functional capacity and design of techniques for the directed intervention and manipulation of the rumen microbiota towards a desired state.

Review: Using physiologically based models to predict population responses to phytochemicals by wild vertebrate herbivores

To understand how foraging decisions impact individual fitness of herbivores, nutritional ecologists must consider the complex in vivo dynamics of nutrient-nutrient interactions and nutrient-toxin interactions associated with foraging. Mathematical modeling has long been used to make foraging predictions (e.g. optimal foraging theory) but has largely been restricted to a single currency (e.g. energy) or using simple indices of nutrition (e.g. fecal nitrogen) without full consideration of physiologically based interactions among numerous co-ingested phytochemicals. Here, we describe a physiologically based model (PBM) that provides a mechanistic link between foraging decisions and demographic consequences. Including physiological mechanisms of absorption, digestion and metabolism of phytochemicals in PBMs allows us to estimate concentrations of ingested and interacting phytochemicals in the body. Estimated phytochemical concentrations more accurately link intake of phytochemicals to changes in individual fitness than measures of intake alone. Further, we illustrate how estimated physiological parameters can be integrated with the geometric framework of nutrition and into integral projection models and agent-based models to predict fitness and population responses of vertebrate herbivores to ingested phytochemicals. The PBMs will improve our ability to understand the foraging decisions of vertebrate herbivores and consequences of those decisions and may help identify key physiological mechanisms that underlie diet-based ecological adaptations.

Addressing Global Ruminant Agricultural Challenges Through Understanding the Rumen Microbiome: Past, Present, and Future

The rumen is a complex ecosystem composed of anaerobic bacteria, protozoa, fungi, methanogenic archaea and phages. These microbes interact closely to breakdown plant material that cannot be digested by humans, whilst providing metabolic energy to the host and, in the case of archaea, producing methane. Consequently, ruminants produce meat and milk, which are rich in high-quality protein, vitamins and minerals, and therefore contribute to food security. As the world population is predicted to reach approximately 9.7 billion by 2050, an increase in ruminant production to satisfy global protein demand is necessary, despite limited land availability, and whilst ensuring environmental impact is minimized. Although challenging, these goals can be met, but depend on our understanding of the rumen microbiome. Attempts to manipulate the rumen microbiome to benefit global agricultural challenges have been ongoing for decades with limited success, mostly due to the lack of a detailed understanding of this microbiome and our limited ability to culture most of these microbes outside the rumen. The potential to manipulate the rumen microbiome and meet global livestock challenges through animal breeding and introduction of dietary interventions during early life have recently emerged as promising new technologies. Our inability to phenotype ruminants in a high-throughput manner has also hampered progress, although the recent increase in “omic” data may allow further development of mathematical models and rumen microbial gene biomarkers as proxies. Advances in computational tools, high-throughput sequencing technologies and cultivation-independent “omics” approaches continue to revolutionize our understanding of the rumen microbiome. This will ultimately provide the knowledge framework needed to solve current and future ruminant livestock challenges.

The Planktonic Core Microbiome and Core Functions in the Cattle Rumen by Next Generation Sequencing

The cow rumen harbors a great variety of diverse microbes, which form a complex, organized community. Understanding the behavior of this multifarious network is crucial in improving ruminant nutrient use efficiency. The aim of this study was to expand our knowledge by examining 10 Holstein dairy cow rumen fluid fraction whole metagenome and transcriptome datasets. DNA and mRNA sequence data, generated by Ion Torrent, was subjected to quality control and filtering before analysis for core elements. The taxonomic core microbiome consisted of 48 genera belonging to Bacteria (47) and Archaea (1). The genus Prevotella predominated the planktonic core community. Core functional groups were identified using co-occurrence analysis and resulted in 587 genes, from which 62 could be assigned to metabolic functions. Although this was a minimal functional core, it revealed key enzymes participating in various metabolic processes. A diverse and rich collection of enzymes involved in carbohydrate metabolism and other functions were identified. Transcripts coding for enzymes active in methanogenesis made up 1% of the core functions. The genera associated with the core enzyme functions were also identified. Linking genera to functions showed that the main metabolic pathways are primarily provided by Bacteria and several genera may serve as a “back-up” team for the central functions. The key actors in most essential metabolic routes belong to the genus Prevotella. Confirming earlier studies, the genus Methanobrevibacter carries out the overwhelming majority of rumen methanogenesis and therefore methane emission mitigation seems conceivable via targeting the hydrogenotrophic methanogenesis.

The rumen microbiome: a crucial consideration when optimising milk and meat production and nitrogen utilisation efficiency

Methane is generated in the foregut of all ruminant animals by the microorganisms present. Dietary manipulation is regarded as the most effective and most convenient way to reduce methane emissions (and in turn energy loss in the animal) and increase nitrogen utilization efficiency. This review examines the impact of diet on bovine rumen function and outlines what is known about the rumen microbiome. Our understanding of this area has increased significantly in recent years due to the application of omics technologies to determine microbial composition and functionality patterns in the rumen. This information can be combined with data on nutrition, rumen physiology, nitrogen excretion and/or methane emission to provide comprehensive insights into the relationship between rumen microbial activity, nitrogen utilisation efficiency and methane emission, with an ultimate view to the development of new and improved intervention strategies.

Development and characterization of stable anaerobic thermophilic methanogenic microbiomes fermenting switchgrass at decreasing residence times

BACKGROUND

Anaerobic fermentation of lignocellulose occurs in both natural and managed environments, and is an essential part of the carbon cycle as well as a promising route to sustainable production of fuels and chemicals. Lignocellulose solubilization by mixed microbiomes is important in these contexts.

RESULTS

Here, we report the development of stable switchgrass-fermenting enrichment cultures maintained at various residence times and moderately high (55 °C) temperatures. Anaerobic microbiomes derived from a digester inoculum were incubated at 55 °C and fed semi-continuously with medium containing 30 g/L mid-season harvested switchgrass to achieve residence times (RT) of 20, 10, 5, and 3.3 days. Stable, time-invariant cellulolytic methanogenic cultures with minimal accumulation of organic acids were achieved for all RTs. Fractional carbohydrate solubilization was 0.711, 0.654, 0.581 and 0.538 at RT = 20, 10, 5 and 3.3 days, respectively, and glucan solubilization was proportional to xylan solubilization at all RTs. The rate of solubilization was described well by the equation r = k(C - C₀f_r), where C represents the concentration of unutilized carbohydrate, C₀ is the concentration of carbohydrate (cellulose and hemicellulose) entering the bioreactor and f_r is the extrapolated fraction of entering carbohydrate that is recalcitrant at infinite residence time. The 3.3 day RT is among the shortest RT reported for stable thermophilic, methanogenic digestion of a lignocellulosic feedstock. 16S rDNA phylotyping and metagenomic analyses were conducted to characterize the effect of RT on community dynamics and to infer functional roles in the switchgrass to biogas conversion to the various microbial taxa. Firmicutes were the dominant phylum, increasing in relative abundance from 54 to 96% as RT decreased. A Clostridium clariflavum strain with genetic markers for xylose metabolism was the most abundant lignocellulose-solubilizing bacterium. A Thermotogae (Defluviitoga tunisiensis) was the most abundant bacterium in switchgrass digesters at RT = 20 days but decreased in abundance at lower RTs as did multiple Chloroflexi. Synergistetes and Euryarchaeota were present at roughly constant levels over the range of RTs examined.

CONCLUSIONS

A system was developed in which stable methanogenic steady-states were readily obtained with a particulate biomass feedstock, mid-season switchgrass, at laboratory (1 L) scale. Characterization of the extent and rate of carbohydrate solubilization in combination with 16S rDNA and metagenomic sequencing provides a multi-dimensional view of performance, species composition, glycoside hydrolases, and metabolic function with varying residence time. These results provide a point of reference and guidance for future studies and organism development efforts involving defined cultures.

Metagenomic binning reveals the functional roles of core abundant microorganisms in twelve full-scale biogas plants

The aim of this work was to elucidate the microbial ecology in twelve mesophilic and thermophilic full-scale biogas plants using a genome-centric metagenomic approach. In this study both biogas plants treating manure and those treating sludge from waste water treatment plants were considered. The identification of 132 Metagenome-Assembled Genomes (MAGs) and analysis of their abundance profile in different samples allowed the identification of the most abundant core members of the anaerobic digestion microbiome. Canonical correspondence analysis was used to determine the influence of biotic and environmental factors on MAGs abundance and to investigate the methanogenic performance of the biogas plants. Prediction of the functional properties of MAGs was obtained analyzing their KEGG pathways and their carbohydrate active domains. Network analysis allowed investigation of species-species associations and shed light on syntrophic interactions between members belonging to the anaerobic digestion dark matter (phylum Fermentibacteria). By stratifying and comparing different levels of information, it was predicted that some MAGs have a crucial role in the manure-supplemented thermophilic biogas plants and it was highlighted the importance of the glycine cleavage system in complementing the "truncated" Wood-Ljungdahl pathway.

Pseudomonas PB1-Like Phages: Whole Genomes from Metagenomes Offer Insight into an Abundant Group of Bacteriophages

Despite the abundance, ubiquity and impact of environmental viruses, their inherent genomic plasticity and extreme diversity pose significant challenges for the examination of bacteriophages on Earth. Viral metagenomic studies have offered insight into broader aspects of phage ecology and repeatedly uncover genes to which we are currently unable to assign function. A combined effort of phage isolation and metagenomic survey of Chicago&rsquo;s nearshore waters of Lake Michigan revealed the presence of Pbunaviruses, relatives of the Pseudomonas phage PB1. This prompted our expansive investigation of PB1-like phages. Genomic signatures of PB1-like phages and Pbunaviruses were identified, permitting the unambiguous distinction between the presence/absence of these phages in soils, freshwater and wastewater samples, as well as publicly available viral metagenomic datasets. This bioinformatic analysis led to the de novo assembly of nine novel PB1-like phage genomes from a metagenomic survey of samples collected from Lake Michigan. While this study finds that Pbunaviruses are abundant in various environments of Northern Illinois, genomic variation also exists to a considerable extent within individual communities.

Glucuronoyl esterases: diversity, properties and biotechnological potential. A review

Glucuronoyl esterases (GEs) belonging to the carbohydrate esterase family 15 (CE15) are involved in microbial degradation of lignocellulosic plant materials. GEs are capable of degrading complex polymers of lignin and hemicellulose cleaving ester bonds between glucuronic acid residues in xylan and lignin alcohols. GEs promote separation of lignin, hemicellulose and cellulose which is crucial for efficient utilization of biomass as an energy source and feedstock for further processing into products or chemicals. Genes encoding GEs are found in both fungi and bacteria, but, so far, bacterial GEs are essentially unexplored, and despite being discovered >10 years ago, only a limited number of GEs have been characterized. The first laboratory scale example of improved xylose and glucuronic acid release by the synergistic action of GE with cellulolytic enzymes was only reported recently (improved C5 sugar and glucuronic acid yields) and, until now, not much is known about their biotechnology potential. In this review, we discuss the diversity, structure and properties of microbial GEs and consider the status of their action on natural substrates and in biological systems in relation to their future industrial use.

Higher order scaffoldin assembly in Ruminococcus flavefaciens cellulosome is coordinated by a discrete cohesin-dockerin interaction

Cellulosomes are highly sophisticated molecular nanomachines that participate in the deconstruction of complex polysaccharides, notably cellulose and hemicellulose. Cellulosomal assembly is orchestrated by the interaction of enzyme-borne dockerin (Doc) modules to tandem cohesin (Coh) modules of a non-catalytic primary scaffoldin. In some cases, as exemplified by the cellulosome of the major cellulolytic ruminal bacterium Ruminococcus flavefaciens, primary scaffoldins bind to adaptor scaffoldins that further interact with the cell surface via anchoring scaffoldins, thereby increasing cellulosome complexity. Here we elucidate the structure of the unique Doc of R. flavefaciens FD-1 primary scaffoldin ScaA, bound to Coh 5 of the adaptor scaffoldin ScaB. The RfCohScaB5-DocScaA complex has an elliptical architecture similar to previously described complexes from a variety of ecological niches. ScaA Doc presents a single-binding mode, analogous to that described for the other two Coh-Doc specificities required for cellulosome assembly in R. flavefaciens. The exclusive reliance on a single-mode of Coh recognition contrasts with the majority of cellulosomes from other bacterial species described to date, where Docs contain two similar Coh-binding interfaces promoting a dual-binding mode. The discrete Coh-Doc interactions observed in ruminal cellulosomes suggest an adaptation to the exquisite properties of the rumen environment.

Assembly of 913 microbial genomes from metagenomic sequencing of the cow rumen

The cow rumen is adapted for the breakdown of plant material into energy and nutrients, a task largely performed by enzymes encoded by the rumen microbiome. Here we present 913 draft bacterial and archaeal genomes assembled from over 800 Gb of rumen metagenomic sequence data derived from 43 Scottish cattle, using both metagenomic binning and Hi-C-based proximity-guided assembly. Most of these genomes represent previously unsequenced strains and species. The draft genomes contain over 69,000 proteins predicted to be involved in carbohydrate metabolism, over 90% of which do not have a good match in public databases. Inclusion of the 913 genomes presented here improves metagenomic read classification by sevenfold against our own data, and by fivefold against other publicly available rumen datasets. Thus, our dataset substantially improves the coverage of rumen microbial genomes in the public databases and represents a valuable resource for biomass-degrading enzyme discovery and studies of the rumen microbiome.

A metagenomic analysis of the camel rumen's microbiome identifies the major microbes responsible for lignocellulose degradation and fermentation

BACKGROUND

The diverse microbiome present in the rumen of ruminant animals facilitates the digestion of plant-based fiber. In this study, a shotgun metagenomic analysis of the microbes adhering to plant fiber in the camel rumen was undertaken to identify the key species contributing to lignocellulose degradation and short chain volatile fatty acids (VFA) fermentation.

RESULTS

The density of genes in the metagenome encoding glycoside hydrolases was estimated to be 25 per Mbp of assembled DNA, which is significantly greater than what has been reported in other sourced metagenomes, including cow rumen. There was also a substantial representation of sequences encoding scaffoldins, dockerins and cohesins, indicating the potential for cellulosome-mediated lignocellulose degradation. Binning of the assembled metagenome has enabled the definition of 65 high-quality genome bins which showed high diversity for lignocellulose degrading enzymes. Species associated to Bacteroidetes showed a high proportion of genes for debranching and oligosaccharide degrading enzymes, while those belonging to Firmicutes and Fibrobacteres were rich in cellulases and hemicellulases and thus these lineages were probably the key for ensuring the degradation of lignocellulose. The presence of many "polysaccharide utilization loci" (PULs) in Bacteroidetes genomes indicates their broad substrate specificity and high potential carbohydrate degradation ability. An analysis of VFA biosynthesis pathways showed that genes required for the synthesis of acetate were present in a range of species, except for Elusimicrobiota and Euryarchaeota. The production of propionate, exclusively via the succinate pathway, was carried out by species belonging to the phyla Bacteroidetes, Firmicutes, Spirochaetes and Fibrobacteres. Butyrate was generated via the butyrylCoA: acetate CoA-transferase pathway by Bacteroidetes and Lentisphaerae species, but generally via the butyrate kinase pathway by Firmicutes species.

CONCLUSION

The analysis confirmed the camel rumen's microbiome as a dense and yet largely untapped source of enzymes with the potential to be used in a range of biotechnological processes including biofuel, fine chemicals and food processing industries.

Comparative Metagenomics of Cellulose- and Poplar Hydrolysate-Degrading Microcosms from Gut Microflora of the Canadian Beaver (Castor canadensis) and North American Moose (Alces americanus) after Long-Term Enrichment

To identify carbohydrate-active enzymes (CAZymes) that might be particularly relevant for wood fiber processing, we performed a comparative metagenomic analysis of digestive systems from Canadian beaver (Castor canadensis) and North American moose (Alces americanus) following 3 years of enrichment on either microcrystalline cellulose or poplar hydrolysate. In total, 9,386 genes encoding CAZymes and carbohydrate-binding modules (CBMs) were identified, with up to half predicted to originate from Firmicutes, Bacteroidetes, Chloroflexi, and Proteobacteria phyla, and up to 17% from unknown phyla. Both PCA and hierarchical cluster analysis distinguished the annotated glycoside hydrolase (GH) distributions identified herein, from those previously reported for grass-feeding mammals and herbivorous foragers. The CAZyme profile of moose rumen enrichments also differed from a recently reported moose rumen metagenome, most notably by the absence of GH13-appended dockerins. Consistent with substrate-driven convergence, CAZyme profiles from both poplar hydrolysate-fed cultures differed from cellulose-fed cultures, most notably by increased numbers of unique sequences belonging to families GH3, GH5, GH43, GH53, and CE1. Moreover, pairwise comparisons of moose rumen enrichments further revealed higher counts of GH127 and CE15 families in cultures fed with poplar hydrolysate. To expand our scope to lesser known carbohydrate-active proteins, we identified and compared multi-domain proteins comprising both a CBM and domain of unknown function (DUF) as well as proteins with unknown function within the 416 predicted polysaccharide utilization loci (PULs). Interestingly, DUF362, identified in iron-sulfur proteins, was consistently appended to CBM9; on the other hand, proteins with unknown function from PULs shared little identity unless from identical PULs. Overall, this study sheds new light on the lignocellulose degrading capabilities of microbes originating from digestive systems of mammals known for fiber-rich diets, and highlights the value of enrichment to select new CAZymes from metagenome sequences for future biochemical characterization.