Metabolic influence of core ciliates within the rumen microbiome

News

HomeHome / News / Metabolic influence of core ciliates within the rumen microbiome

Jan 13, 2024

Metabolic influence of core ciliates within the rumen microbiome

The ISME Journal (2023)Cite

The ISME Journal (2023)Cite this article

2034 Accesses

32 Altmetric

Metrics details

Protozoa comprise a major fraction of the microbial biomass in the rumen microbiome, of which the entodiniomorphs (order: Entodiniomorphida) and holotrichs (order: Vestibuliferida) are consistently observed to be dominant across a diverse genetic and geographical range of ruminant hosts. Despite the apparent core role that protozoal species exert, their major biological and metabolic contributions to rumen function remain largely undescribed in vivo. Here, we have leveraged (meta)genome-centric metaproteomes from rumen fluid samples originating from both cattle and goats fed diets with varying inclusion levels of lipids and starch, to detail the specific metabolic niches that protozoa occupy in the context of their microbial co-habitants. Initial proteome estimations via total protein counts and label-free quantification highlight that entodiniomorph species Entodinium and Epidinium as well as the holotrichs Dasytricha and Isotricha comprise an extensive fraction of the total rumen metaproteome. Proteomic detection of protozoal metabolism such as hydrogenases (Dasytricha, Isotricha, Epidinium, Enoploplastron), carbohydrate-active enzymes (Epidinium, Diplodinium, Enoploplastron, Polyplastron), microbial predation (Entodinium) and volatile fatty acid production (Entodinium and Epidinium) was observed at increased levels in high methane-emitting animals. Despite certain protozoal species having well-established reputations for digesting starch, they were unexpectedly less detectable in low methane emitting-animals fed high starch diets, which were instead dominated by propionate/succinate-producing bacterial populations suspected of being resistant to predation irrespective of host. Finally, we reaffirmed our abovementioned observations in geographically independent datasets, thus illuminating the substantial metabolic influence that under-explored eukaryotic populations have in the rumen, with greater implications for both digestion and methane metabolism.

Ruminants operate in symbiosis with their intrinsic rumen microbiome, which is responsible for the degradation of forage into nutrients, in the form of volatile fatty acids (VFAs), supplying ~70% of net energy for the host [1]. The rumen microbiome itself is a complex assemblage of bacterial, fungal, archaeal, viral and protozoal microorganisms whose intricate composition and function is connected to host productivity traits, such as feed efficiency, milk yield, animal health and greenhouse gas (GHG) emissions [2,3,4,5]. Large collaborative research efforts have been made to identify and characterize the core rumen microbiome including creating a publicly available catalog for cultivated and sequenced genomes [6,7,8]. In the rumen, bacteria are estimated to constitute 50–90%, protozoa 10–50%, fungi 5–10% and archaea <4% of the total microbial biomass [9, 10]. Due to the difficulties of axenically culturing rumen eukaryotic populations and their complex genomic features that are obstinate to current metagenomic assembly and binning technologies, the reconstruction of the rumen microbiome has been heavily biased towards bacterial and archaeal members, whereas the fungal and protozoal contributions of the rumen currently remain poorly characterized. While anaerobic fungi have a reputable role as fiber degraders in the rumen, only 18 anaerobic gut fungi from herbivores are currently described, with only 11 genomes available [11,12,13]. Similarly, to date few rumen protozoal genomes are sequenced and publicly available, chief among them, the rumen ciliate protozoa Entodinium caudatum [14]. More recently, single-cell amplified genomes (SAGs) have been recovered from rumen microbiome samples representative of 5 holotrich and 14 entodiniomorph species spanning 13 genera including Isotricha, Dasytricha, Diplodinum, Enoploplastron, Metadinium, Eremoplastron, Ostracodinium, Polyplastron, Ophryoscolex, Epidinum and Entodinium [15].

Entodinium and Epidinium are entodiniomorphs of the order Entodiniomorphida and represent two of the most dominant genera of rumen protozoa, previously being detected in more than 99% of 592 rumen samples at a mean protozoal relative abundance of ~38 and 16% respectively (2015 rumen census: 32 animal species, 35 countries) [16]. Entodiniomorphs have been previously observed to stimulate methane production [17], but while some possess hydrogenosomes such as Epidinium caudatum [18], it is believed the numerically dominant Entodinium caudatum does not, and instead encodes mitosomes and iron hydrogenases that may indicate hydrogen production, beneficial for methanogenic endosymbionts [19]. Holotrich species of the order Vestibuliferida such as Dasytricha and Isotricha are generally found at lower abundances compared to entodiniomorphs but are renowned for their fermentation of soluble carbohydrates and their production of VFAs and hydrogen [20, 21]. Here, we present a genome-centric metaproteomics analysis of the rumen microbiome from two different host species; Holstein dairy cattle (Bos taurus) and alpine goats (Capra hircus), that were fed diets of first-cut grassland hay with a 45:55 forage:concentrate ratio, with concentrates supplemented with either no additional lipid (CTL), or corn oil and wheat starch (COS) [5, 22, 23]. Previous microscopy analysis revealed high counts of entodiniomorphs and to a lesser extent holotrichs, while metadata revealed that animals fed COS displayed reduced methane emissions, irrespective of host. To describe how these diets affect digestion and production of methane and VFAs, we sought to investigate changes in function and composition in the complex rumen microbiome. By using shotgun metagenomic sequencing we recovered in total 244 prokaryote metagenome-assembled genomes (MAGs) that together with selected isolate- and SAG-derived eukaryote genomes [14, 15, 24,25,26,27,28,29] formed the database for our genome-centric metaproteomic analysis of rumen protozoal species. Key protozoal metabolisms were detected related to plant fiber and starch degradation, bacterial predation as well as hydrogen and VFA production. Despite numerous protozoal species having the genetic ability to degrade starch, our analysis showed contrasting data that suggests rumen protozoa were less metabolically active in the rumen microbiome of animals fed a starch-rich diet. In such a scenario other starch-degrading and/or propionate and succinate producing bacterial genera, as Prevotella and Fibrobacter and members of the families Succinivibrionaceae and Aminobacteriaceae appeared to be more prevalent. In concert, our analysis showed that reduced methane production in both cattle and goats eating feeds supplemented with COS is likely caused via a redirection of hydrogen to succinate and propionate production instead of methanogenesis. Finally, by analysing a secondary, geographically independent dataset, we reaffirmed our primary observations of protozoal dominance and patterns concerning hydrogen, plant fiber and starch-related metabolism, thus supporting our hypothesis that this protozoal species plays a core role in rumen microbiome function.

The experimental procedures were approved by the Auvergne-Rhône-Alpes Ethics Committee for Experiments on Animals (France; DGRI agreement APAFIS#3277–2015121411432527 v5) and complied with the European Union Directive 2010/63/EU guidelines. Experiments were performed at the animal experimental facilities of HerbiPôle site de Theix at the Institut National de la Recherche pour l’Agriculture, l’Alimentation l’Environnement (INRAE, Saint-Genès-Champanelle, France) from February to July 2016. Experimental design, animals and diets were as described elsewhere [5, 23]. Briefly, four Holstein cows and four Alpine goats, all lactating, were enrolled in respectively two 4 × 4 Latin square design trials to study the effects of 4 diets over four 28-d experimental periods (Supplementary Table S1). The original study included four different experimental grassland hay basal diets with concentrates supplemented with various lipid sources; control diet with no added lipids (CTL), diet supplemented with corn oil and wheat starch (COS), diet supplemented with marine algae powder (MAP) and diet supplemented with hydrogenated palm oil (HPO) (Supplementary Table S2) [23]. In the present study, we focused on the CTL and COS diets, which were associated with the most extreme methane (CH4) emission profiles in both ruminant species. The CTL diet composed of grass hay ad libitum with concentrates containing no additional lipid, whereas COS contained corn oil (5.0% total dry matter intake (DMI)) and wheat starch −5.0% of total DMI (COS) (Table 1). Corn oil (Olvea, Saint Léonard, France) was added to the concentrate, at 5% of total DMI and contained (g/kg of total FA): 16:0 (114), 18:0 (16.4), cis-9 18:1 (297), cis-11 18:1 (6.30), 18:2n-6 (535), 18:3n-3 (7.57), 20:0 (3.48), 22:0 (1.0), 24:0 (1.5) and total FA (1000 g/kg). Detailed diet composition is available in Martin et al. [5]. Each experimental period lasted for 28 days. Diets were offered as two equal meals at 0830 and 1600h. Animals had access to a constant supply of freshwater ad libitum. Concentrate and hay refusals were weighed daily. The amounts of feed offered the following day was adjusted regarding to refusals to maintain the targeted the dietary 45% dry matter (DM) forage and 55% DM concentrate ratio. We acknowledge that by only including one sampling point for rumen fluid we are reducing the power that the original experimental Latin square design gives regarding minimizing individual error between the animals.

Rumen fluid was collected through stomach-tubing before the morning feeding on day 27 of each experimental period. The stomach tube consisted of a flexible 23 mm external diameter PVC hose fitted to a 10 cm-strainer at the head of the probe for cattle, and a flexible 15 mm PVC hose with a 12 cm-strainer for goats. The first 200 ml of rumen fluid was discarded from to minimize contamination from saliva. Samples were filtered through a polyester monofilament fabric (280 μm pore size), dispatched in 2-ml screw-cap tubes, centrifuged at 15,000 × g for 10 min and the pellet snap-frozen in liquid nitrogen. Samples were stored at −80 °C until DNA extraction using the Yu and Morrison bead-beating procedure [30]. In total, 32 rumen fluid samples (four cows and four goats fed four diets included in the original study [23]) were sent to the Norwegian University of Life Sciences (NMBU) for metagenomic and metaproteomic analysis. Respiration chambers were used to measure methane emissions over a 5-day period, while VFA and NH3 concentrations were determined by gas chromatography using a flame ionization detector [31]. Protozoa were counted by microscopy and categorized as either small entodiniomorphs (<100 µm), large entodiniomorphs (>100 µm) or as holotrichs Dasytricha and Isotricha [9]. Further specifics about diets and measurements can be found in Martin et al. [5] and VFA and methane measurements are summarized in Table 1.

Metagenomic shotgun sequencing was performed at the Norwegian Sequencing Center on two lanes of the HiSeq 3/4000 (Illumina) generating 150 bp paired-end reads in both lanes. Sequencing libraries were prepared using the TruSeq DNA PCR-Free High Throughput Library Prep Kit (Illumina) prior to sequencing. All 32 samples (four cows and four goats fed four diets included in the original study [23]) were run on both lanes to prevent potential lane-to-lane sequencing bias. FASTQ files were quality filtered and Illumina adapters removed using Trimmomatic [32] (v. 0.36) with parameters -phred33 for base quality encoding, leading and trailing base threshold set to 20. Sequences with an average quality score below 15 in a 4-base sliding window were trimmed and the minimum length of reads was set to 36 bp. MEGAHIT [33] (v.1.2.9) was used to co-assemble reads originating from samples collected from cattle and goats separately, with options –kmin-1pass, --k-list 27,37,47,57,67,77,87, --min-contig-len 1000 in accordance with [7]. Bowtie2 [34] (v. 2.3.4.1) was used to map reads back to the assemblies and SAMtools [35] (v. 1.3.1) was used to convert SAM files to BAM format and index sorted BAM files.

The two co-assemblies (one from the samples originating from cattle and the other originating from the samples of goats) were binned using Maxbin2 [36], MetaBAT2 [37] and CONCOCT [38]. MetaBAT2 (v. 2.12.1) was run using parameters –minContig 2000 and –numThreads 4, Maxbin2 (v. 2.2.7) ran with default parameters and -thread 4, min_contig_length 2000, and CONCOCT (v. 1.1.0) ran with default parameters and –length_threshold 2000. Further, bins were filtered, dereplicated and aggregated using DASTool [39] (v. 1.1.2) with the parameters –write_bins 1, --threads 2 and BLAST [40] as search engine. This resulted in a total of 244 dereplicated MAGs across the two host species (104 originating from cattle and 140 from goats). CheckM [41] (v. 1.1.3) lineage workflow was used to calculate completeness and contamination of each MAG, with parameters –threads 8, --extension fa, and CoverM (v. 0.5.0) (https://github.com/wwood/CoverM) was used to calculate relative abundance of each MAG, while GTDB-tk [42, 43] (v. 1.3.0) was used for taxonomic annotation. Approximately 90% (219 of 244) of the recovered MAGs were considered high or medium quality MAGs according to MIMAGs threshold for completeness and contamination for genome reporting standards [44] (Supplementary Table S3). Gene calling and functional annotation of the final MAGs were performed using the DRAM [45] pipeline with the databases dbCAN [46], Pfam [47], Uniref90 [48], Merops [49], VOGdb and KOfam [50]. The translated amino acid sequences from the publicly available drafted En. caudatum macronucleus genome were annotated with the KEGG metabolic pathway database using BlastKOALA [51] by Park et al. [14]. Proteins with resulting KEGG Orthology (KO) numbers were functionally assigned to metabolic pathways using KEGG Mapper Reconstruct Tool [52].

The resulting protein sequences for each MAG, as well as those from the host genomes of goat (Capra hircus, NCBI ID: 10731) and cattle (Bos taurus, NCBI ID: 82) were compiled into two databases, from now on referred to as sample specific RUmen DataBase for Goat (RUDB-G) and sample specific RUmen DataBase for Cattle (RUDB-C). In addition, both databases were supplemented with the genome of the protozoa Entodinium caudatum [14], as well as 18 SAGs from Li et al. [15] including Dasytricha ruminantium SAG3, Diplodinium dentatum SAGT1, Diplodinium flabellum SAG1, Enoploplastron triloricatum SAGT1, Entodinium bursa SAG3, Entodinium longinucleatum SAG4, Epidinium cattanei SAG3, Epidinium caudatum SAG1, Eremoplastron rostratum SAG2, Isotricha intestinalis SAGT2, Isotricha prostoma SAG3, Isotricha sp YL-2021a SAG1, Isotricha sp YL-2021b SAG3, Metadinium minomm SAG1, Ophryoscolex caudatus SAGT3, Ostracodinium dentatum SAG1, Ostracodinium gracile SAG1 and Polyplastron multivesiculatum SAGT3. The genome of Fibrobacter succinogenes S85 (NCBI ID: 932) was also added to both rumen databases since it is well recognized as a primary cellulolytic bacterium in the rumen microbiome and has previously been observed as an active microorganism in similar studies yet was not thoroughly binned as a unique MAG in this study. Nine available fungal genomes downloaded from Joint Genome Institute (JGI) Mycocosm [53] were added to RUDB-C (Anaeromyces sp. S4 [25], Caecomyces churrovis [28], Neocallimastix californiae [25], Neocallimastix lanati [29], Piromyces finnis [25], Piromyces sp. E2 [25], Piromyces UH3-1 [24, 25, 27], Piromyces eukMAG [54], Orpinomyces sp [26]). In total the complete databases consisted of 1 291 083 and 1 190 938 protein entries for RUDB-G and RUDB-C, respectively.

To 300 μL of rumen fluid sample (in total 32 rumen fluid samples originating from four cows and four goats fed four diets included in the original study [23]) 150 μL lysis buffer (30 mM DTT, 150 mM Tris-HCl (pH = 8), 0.3% Triton X-100, 12% SDS) was added together with 4 mm glass beads (≤160 μm), followed by brief vortexing and resting on ice for 30 mins. Lysis was performed using FastPrep-24 Classic Grinder (MP Biomedical, Ohio, USA) for 3 × 60 s at 4.0 m/s [55]. Samples were centrifuged at 16,000 × g for 15 min at 4 °C and lysate was carefully removed. Protein concentration was measured using the Bio-Rad DC Protein Assay (Bio-Rad, California USA) with bovine serum albumin as standard. Absorbance of sample lysates was measured at A750 on BioTek Synergy H4 Hybrid Microplate Reader (Thermo Fisher Scientific Inc., Massaschusetts, USA). 40–50 μg of protein was prepared in SDS-buffer, heated in water bath for 5 min at 99 °C and analysed by SDS-PAGE using Any-kD Mini-PROTEAN TGX Stain-Free gels (Bio-Rad, California, USA) in a 2-minute run for sample clean-up purposes, before staining with Coomassie Blue R-250. The visible bands were carefully excised from the gel and divided as 1 × 1 mm pieces before reduction, alkylation and digestion with trypsin. Peptides were concentrated and eluted using C18 ZipTips (Merck Millipore, Darmstadt, Germany) according to manufacturer's instructions, dried and analysed by nano-LC-MS/MS using a Q-Exactive hybrid quadrapole Orbitrap MS (Thermo Fisher Scientific Inc., Massaschusetts, USA) as previously described [56].

Mass spectrometry (MS) raw data were analysed with FragPipe version 19 and searched against the sample-specific protein sequence database (1 291 083 and 1 190 938 protein entries for RUDB-G and RUDB-C, respectively) with MSFragger [57] (Supplementary Table S4). The database was supplemented with contaminant protein entries, such as human keratin, trypsin and bovine serum albumin, in addition to reversed sequences of all protein entries for estimation of false discovery rates (FDR) [58]. Oxidation of methionine and protein N-terminal acetylation were used as variable modifications, while carbomidomethylation of cysteine residues were used as fixed modification. Trypsin was chosen as digestive enzyme, maximum missed cleavages allowed was one and matching tolerance levels for both MS and MS/MS were 20 ppm. The results were filtered to 1% FDR and quantification was done using IonQuant [59] including normalization between samples and the feature ‘match between runs’ to maximize protein identifications. Perseus [60] version 1.6.2.3 was used for further analysis. A protein group was considered valid if it was quantified in at least 50% of the replicates in at least one condition. Protein groups identified as potential contaminants were removed. Calculations of individual genome/SAG/MAG contributions were calculated using the workflow outlined by Kleiner et al. [61], which sums protein abundances (label free quantification (LFQ) values) for all proteins assigned to each genome/SAG/MAG, and differential abundance between diets were detected by an unpaired two-sided Student's t-test (p value < 0.05).

After filtration, we resolved 3657 unique protein groups across the 16 samples from cattle and 3818 unique protein groups across 15 samples from goats (Supplementary Table S5 and Supplementary Table S6). Dot plots for Fig. 1 were made with ggplot2 [62] in R (v. 4.2.0) [63]. To determine which expressed metabolic pathways En. caudatum were significantly enriched for in each diet/animal, we used the hypeR() function from the hyperR package [64] in R which employs the hypergeometric test for significance. The ‘geneset’ for hyperR was generated by using the keggGet() functionof KEGGREST R package to retrieve entries from the KEGG database and determine which pathways the En. caudatum KOs belong to. The geneset was then manually curated to only include metabolic pathways of interest (i.e., we remove pathways such as "Huntington disease"). For the "background" setting in hyperR(), to be conservative, we used the total number of unique KOs (7592) in the En. caudatum genome that could possibly be expressed. Significantly enriched KOs with p value <0.05 were manually filtered. Both p values and FDR-adjusted p values are available in Supplementary Table S7, with FDR- adjusted p values <0.05 highlighted in bold.

Dot plots for total proteins identified (a), and average recovered metaproteomic expression (b presented as proportion of summed LFQ intensities) belonging to protozoal, bacterial, archaeal, or fungal species across the control diet (CTL) and diets supplemented with corn oil and wheat starch (COS) for dairy cattle (n = 4) and goats (n = 4). Detected protein abundances for protozoal and bacterial populations can be found in Supplementary Table S6.

Samples were also analysed from previously performed feeding experiments with Holstein Friesian bulls [65]. In brief, these bulls were subjected to either ad libitum or restricted feeding regime in a compensatory growth model detailed in Keogh et al. [65]. Both feeding groups received the same ratio of concentrate and grass silage, respectively 70 and 30%, of which the concentrate was mainly composed of starch-rich rolled barley (72.5%) and soya (22.5%). Rumen samples were collected at slaughter and stored at −80 °C prior to metagenomics and metaproteomic analysis in this study.

Sample preparation, cell lysis and extraction of DNA was carried out as previously described by McCabe et al. [66]. Quality check of fastq files and removal of low-quality reads was performed using fastp (V.0.19.5). Sequence reads were mapped against the bovine genome (ARS-UCD1.3) using minimap2 (V.2.16), and host sequences were removed. Reads were co-assembled using Megahit (V1.2.6) with "−meta-large" pre-set option as the metagenome was complex. Metagenomic binning was applied to the co-assembly using MetaBAT2 using standard parameters (V.2.12.1). MAGs were then dereplicated using dRep (V.1.4.3), and the resulting MAGs were taxonomically assigned using Bin Annotation Tool (BAT), available on (https://github.com/dutilh/CAT). This tool internally uses prodigal (V.2.6.3) for gene prediction and DIAMOND (V.0.9.14) for the alignment against the non-redundant (nr) protein database (as of Feb 2020).

Sample preparation for metaproteomics was done by lysing cells with bead beating with two glass bead sizes (≤106 µm and 0.5 mm), in 100 mM Tris, pH8, 5% SDS and 10 mM DTT. A FastPrep 24 instrument was operated for 3 × 45 s at a speed of 6.5 m/s. The samples were centrifuged for 15 min at 20.000 × g and the protein extracts were cleaned by Wessel-Flügge precipitation [67]; pellets were dissolved in 5% SDS, 100 mM Tris-Cl, pH8, 10 mM DTT and kept at −20 °C until further processing. Protein digestion was performed using suspension trapping (STrap) [68], dried in a SpeedVac (Eppendorf Concentrator Plus) and re-dissolved in 0.05 % trifluoroacetic acid, 2% acetonitrile for peptide concentration estimation using a Nanodrop One instrument and subsequent MS/MS-analysis. The samples were analyzed using an Ultimate3000 RSLCnano UHPLC coupled to a QExactive hybrid quadrupole-orbitrap mass spectrometer (Thermo Fisher, Bremen, Germany) as described previously [56]. MS raw data were analysed with FragPipe and searched against the Holstein Friesian bulls sample-specific protein sequence database (2 533 362 protein entries) using the same workflow described above.

Because of their large size protozoal species can comprise a significant fraction of the microbial biomass in the rumen [9]. While the total count and diversity of protozoal species are lesser than their bacterial counterparts in the rumen, their genome size and total number of genes encoded within are considerably larger. In addition, due to alternative splicing and post-translational modifications, the protein representation, or proteome, of protozoal populations will be larger than the number of genes in their protozoal genome [69]. Thus, the amount of protein in a protozoa species can be expected to far exceed the amount of protein in a bacterial species that can be identified and quantified in proteomic studies. In this context, our metaproteomic data showed an extensive fraction of detectable proteins affiliated to protozoal species, and other closely related species, in both dairy cattle and goats in proportion to the combined bacterial species that were represented in our genome databases (Fig. 1a). Protein detection intensity of proteins assigned to protozoa were also proportionally greater than that of the bacterial fraction of the rumen microbiome in animals fed the control (CTL) diet (Fig. 1b) further supporting the significant levels of protozoal detection in our samples. Proteomes from entodiniomorph species were detected at high levels in both animal hosts, with Ep. cattanei and En. caudatum dominating the rumen metaproteome of cattle and goats respectively (Fig. 2). In cattle we also detected entodinomorphs D. dentatum, En. bursa and Ep. caudatum (Fig. 2a, Supplementary Table S6) at high levels whereas in goats Entodinium species En. caudatum, En. longinucleatum and En. bursa where among the highest detected entodinomorphs (Fig. 2c). More than three times as many proteins from En. caudatum were detected in goats than cattle, however this was somewhat expected given the 7x higher counts of entodiniomorph concentration (cells/mL) previously observed in the goat samples compared to cattle (Table 1) [5]. In agreement with microscopic cell counts of holotrichs, we detected proteins from both Dasytricha spp. and Isotricha spp (Table 1) in both cattle and goats (Fig. 2a, c).

The figure displays metabolically active populations (as genomes, SAGs or MAGs), with selected expressed proteins (presented as Enzyme Commission (EC) number with short enzyme descriptions) active in fiber/starch degradation, glycolysis and production of pyruvate, butyrate, acetate, and succinate in cattle (a, b) and goats (c, d) fed CTL or COS diets. a, c depict protozoal proteomes that were detected in cattle and goats respectively are presented separately to bacteria (b, d) as the scale of their protein quantification values were ~10× larger. Protein quantification values (y-axis) were calculated by considering both the number of proteins detected per MAG/SAG/genome and their LFQ intensity: we averaged LFQ intensities for each detected protein across biological replicates for each dietary condition (CTL: green or COS: orange), which were subsequently summed for all detected proteins per MAG/SAG/genome. Heatmaps show selected MAG/SAG/genome with metabolically active proteins, presented as EC numbers, recovered from cattle (RUDB-C) and goats (RUDB-G) fed either CTL or COS diets. Where average protein abundance level for a given MAG/SAG/genome was significantly different (determined by paired Wilcoxon test, p value <0.05) between the COS and CTL diet, the given MAG/SAG/genome has been marked with a red star. MAGs are presented with their MAG ID and taxonomic annotation from GTDB-tk. Genome annotations and LFQ intensities used to create heatmaps can be found in Supplementary Table S6.

Rumen protozoa have long been characterized as degraders of plant fibers as well as non-structural storage polysaccharides such as starch [70, 71]. Closer examination of the proteomes from both highly detected Epidinium species Ep. cattanei and Ep. caudatum in cattle highlighted a plethora of detected carbohydrate active enzymes (CAZymes) targeting lignocellulosic polysaccharides such as cellulose (GH5, GH9), xylan (GH10, GH11, GH30, CE1), mannan (GH26) and fiber-derived oligosaccharides (GH2, GH3, GH94) (Fig. 2a, Supplementary Table S6). In goats similar CAZy families were detected but from different entodinomorphs species, such as E. triloricatum and P. multivesiculatum (Fig. 2c, Supplementary Table S6). Our in vivo CAZyme observations were in strong agreement to long standing in vitro data, which have shown entodinomorphs such as Epidinium, Polyplastron and Enoploplastron have greater endoglucanase and xylanase activity, while Entodinium spp. have only weak activity [70]. Similarly, it has previously been shown in vitro that holotrichs such as Dasytricha spp. have glucosidase and cellobiosidase activity but negligible fibrolytic activity [72], which our in vivo data also supported (Fig. 2a, c). The composition of protozoal species metabolizing starch (both degrading and synthesizing) also varied between the two animal hosts with goats dominated by Entodinium spp. and D. ruminantium while protozoal starch metabolism in cattle was Epidinium spp., D. ruminantium and Isotricha spp (Fig. 2, Supplementary Table S6). While bacterial CAZymes targeting lignocellulose were also detected (e.g., GH48 cellulases from Ruminococcaceae populations in cattle), CAZymes targeting lignocellulose and starch were detected at higher frequency and protein quantification levels in protozoal species suggesting they have a central role in ruminal fiber digestion.

En. caudatum is renowned for its predatory activity and is acknowledged as the most abundant protozoa in the rumen, whereby it has been estimated that 0.1% of rumen prokaryotes are digested by the rumen protozoal population every minute [73]. The dominance of En. caudatum in our data presented an opportunity to better understand the metabolic influence of this universal species. While previous efforts have investigated the genome and transcriptome of En. caudatum grown in monoculture [14, 74], our metaproteomic analysis sought to reveal in vivo metabolism and functions of En. caudatum within the rumen microbiome. In accordance with Wang et al. [74], our metaproteomic analysis revealed expressed proteins attributed to metabolic pathways such as carbon metabolism, glycolysis/gluconeogenesis, starch and sucrose (and glycogen) metabolism, pyruvate metabolism, oxidative phosphorylation and production of alcohol (Supplementary Table S6). Wang et al. found that En. caudatum uses carbohydrates such as starch as its primary substrate, as well as cellulose and hemicellulose to a certain degree [74], and their transcript analyses showed that En. caudatum had high levels of expression of amylases and low-level expression of hemicellulases, cellulases and pectinases. Similarly, our metaproteomic analysis reveals expression of amylases by En. caudatum that are predicted to enable En. caudatum to engulf and degrade starch granules to simpler sugars and to produce glycogen, its most important storage carbohydrate [75]. However, no detection of En. caudatum CAZymes related to hemicellulose or pectin were observed in any of our metaproteomes, suggesting that it is not engaging in the deconstruction of these carbohydrates at the time our samples were collected for analysis (before feeding). It should be noted that ruminal fermentation activity as well as production of VFAs and methane will be at its highest after feeding, as a result of an increased availability of fermentable substrate [76]. While sampling time can influence the recovered microbial composition and hence function, any differences in metabolic parameters or species abundance in this study is comparable across both diets given the consistent sampling times.

While monoculture cultures of En. caudatum have not been established to verify the VFAs it can produce, Wang et al. found transcripts of enzymes involved in fermentative formation of acetate and butyrate [74]. Similarly, we detected proteins inferred in metabolism of acetate, butyrate, and alcohol in En. caudatum. Goats fed the CTL diet had a significantly higher proportion of proteins from En. caudatum and concurrently had increased relative levels of acetate and butyrate compared to animals fed the corn oil and wheat starch diet (COS), which had fewer En. caudatum proteins and lower acetate/butyrate levels (Fig. 2 and Table 1). As En. caudatum populations were seemingly most abundant in goats fed the CTL diet, we used these metaproteomes to reconstruct metabolic features (Fig. 3). Of the 538 En. caudatum proteins identified in CTL-fed goats, 244 had unique KO numbers assigned, from which KEGG Mapper reconstructions [52] enabled functional assignment of 217 proteins to metabolic pathways. Our metabolic reconstructions showed expressed proteins involved in endocytosis, phagosome and lysosome processes for predatory activity, engulfment and digestion of bacteria (Fig. 3, Supplementary Table S6). For the rumen samples used in this study Martin et al. previously observed higher NH3 concentrations in goats compared to cattle and hypothesized that it might have resulted from increased bacterial protein breakdown and feed protein degradability due to higher density of entodiniomorphs known for their predatory activity [5]. In support of these observations, we performed metaproteomic pathway enrichment analysis of En. caudatum (Fig. 3, Supplementary Table S7) proteins detected in goats, which revealed significantly enriched endocytosis metabolism. Other biological processes such as glycolysis and starch and sucrose metabolism were significantly enriched pathways in En. caudatum (Supplementary Table S7). Although suspected of having metabolic interactions with methanogenic archaea, several protozoal populations such as Entodinium are hypothesized as having associations with certain members of the Gram-negative Gammaproteobacteria, which multiple studies have speculated are resistant to protozoal engulfment [77,78,79]. In contrast, Gutierrez and Davis previously demonstrated that Entodinium-species engulf Gram positive starch-degraders [79]. In the context of our data, we speculate that CTL fed animals provided En. caudatum-like populations optimal conditions for predation, whereas increased starch levels in the COS diets facilitated increases in Gram-negative species that are possibly resistant to protozoal engulfment and/or reduced pH levels, leading to sub-optimal conditions for Entodinium and a decrease in its levels (Fig. 2).

KO identifiers for identified proteins were analysed via KEGG mapper to reconstruct expressed key features in the metabolism of En. caudatum. Dashed arrows represent proteins or pathways that were not detected in our metaproteomes but are key steps in their respective pathways. Detailed information connecting KO identifiers to their respective gene ID, LFQ and animal/diet can be found in Supplementary Table S6.

The changes in VFA and methane levels in animals fed the high starch COS diet, previously measured by Martin et al. [5], suggested significant alterations in composition and thus functions of the rumen microbiome irrespective of host species. In particular, a decrease in proportions of acetate and butyrate, decrease in the acetate:propionate ratio and an increase in relative propionate levels were observed in animals fed the COS diet, compared to the CTL diet (Table 1). Diets that are high in starch content or with low forage:concentrate ratios have previously been shown to result in higher production of propionate and succinate, as they are easily fermented in the rumen and accordingly have high passage rates [80, 81]. We therefore leveraged our genome-centric metaproteomic data from both cattle (Fig. 2a, b) and goats (Fig. 2c, d) fed either the COS or CTL diet to gain an overview of protein expression from individual populations. We specifically focused on pathways involved in the degradation of starch (CTL: corn starch, COS: corn + wheat starch) to pyruvate through glycolysis and finally formation of acetate, butyrate and propionate (via succinate). Irrespective of host, and despite their starch-metabolizing reputation [74, 75], Entodinium, Epidinium, Isotricha and Dasytricha spp had lower abundance and less proteins involved in starch degradation in animals fed the COS diet compared to those fed the CTL diet (Fig. 2a, c). Further, we observed opposing patterns for proteins identified as Entodinium and Epidinium-spp. involved in glycolysis, and production of pyruvate, acetate and butyrate, which were detected in higher levels in both cattle and goats fed the CTL diet compared to the COS diet.

While several putative protozoal amylases were detected across all animals and diets, their quantification levels (i.e., protein abundances) did not increase as expected when higher levels of starch were available (Fig. 2a, c). We therefore hypothesized that the observed shift in VFA profiles in response to increased starch was additionally influenced by the bacterial fraction of the rumen microbiome. In contrast to lower levels of protozoa in the animals fed the COS diet, we observed an increase in bacterial proteome detection irrespective of host (Fig. 4), including suspected starch-degrading and succinate- and propionate-producing bacterial species (Fig. 2b, d). For example, starch fermentation pathways from population genomes affiliated with the Succinivibrionaceae family, Prevotella species and, additionally for goats, members of the Aminobacteriaceae families, were detected at higher proteomic levels in the animals fed the COS diet compared to those fed the CTL diet (Fig. 2b, d).

Volcano plots indicating different rumen microbiome proteins from dairy cattle (a) and goats (b) fed either the control or COS diets and which displayed both large magnitude of fold-changes in LFQ intensities (x axis) and high statistical significance (−log10 of nominal p values using an unpaired t-test, y axis). Dashed horizontal line denotes a p value <0.05 cut-off. For both animal hosts we observed an increase in protein detection for bacterial populations (greenblue) in animals fed high-starch diets compared to protozoal populations (purple), which were detected at higher LFQ intensities in animals fed the control diet.

For the animals sampled in this study, Martin et al. demonstrated a ~25–30% reduction in methane emissions in both cattle and goats fed the COS diet compared to the control (Table 1) [5]. Previous comparisons between defaunated and faunated animals have shown decrease in methane production in protozoa-free ruminants, suggesting symbiotic interactions between methanogenic archaea and protozoal species [77]. Methanogen's epi- and endo-symbiotic relationships with protozoa have also been suggested to contribute to 9–37% of rumen methanogenesis [77, 82,83,84]. Moreover, studying microcosms with the presence and absence of protozoal species Solomon et al. reported higher levels of acetate and butyrate in microcosms with protozoa present in addition to increased methane emissions [77], which supports the main findings of animals fed the CTL diet in this study (Table 1). Several of the protozoal species detected in our metaproteomic data are known to have hydrogenosomes and closer examination showed that [FeFe]-hydrogenases from both holotrich and entodinomorph genera, such as Dasytricha, Isotricha, Diplodinium, Epidinium, Entodinium and Enoploplastron, were detected at higher quantification levels in animals fed the CTL compared to those fed COS (Fig. 2a, c, Supplementary Table S6). We also observed partial evidence in our proteomic data that the hydrogenosome-lacking En. caudatum makes contributions to ruminal hydrogen production, supporting previous reports associating its abundance with higher methane levels [77, 84, 85]. Only one of its eight [FeFe]-hydrogenases were detected in goats (absent in cattle), though it showed no changes in protein abundance in either the high (CTL) or low (COS) methane yielding animals. Given our observations, we speculate that predicted protozoal-methanogen relationships that lead to increased methane levels in this study are centered on hydrogen transfer.

The addition of lipids in diets for ruminants have been largely studied as a methane mitigation strategy [86, 87]. However, it is widely believed that lipids are not fermented in the rumen and rather modified/hydrogenated and hence do not contribute to hydrogen production or methane production in the rumen. While biohydrogenation of unsaturated fatty acids such as corn oil can serve as a hydrogen sink, it has been observed that very small amounts of metabolic hydrogen (1–2%) are used for biohydrogenation [86, 87]. Several studies have also shown that the CH4 mitigation effects of lipids are dependent on both dose and fatty acid composition [82, 88]. For example, medium chain fatty acids (MCFAs), such as lauric acid, have been proven to be a more effective inhibitor of protozoa than long chain polyunsaturated fatty acids (PUFAs), such as corn oil, which mainly consists of linoleic acid [82, 87,88,89]. In a study conducted by Zhang et al. [88], goats fed corn oil as a supplement decreased ruminal hydrogen concentrations and total methane emissions. Nevertheless, there was seemingly no effect on rumen protozoal populations, which suggests that the corn oil dosage used in the Zhang study did not act as an anti-protozoal agent [88]. In addition, high levels of PUFAs in the diet will also likely impact the metabolism of various bacterial and archaeal species [82, 87, 88], particularly keystone fiber degrading Gram-positive bacterial species that produce various levels of VFAs and hydrogen [90, 91]. Collectively, these previous findings were in agreement with the decreased CH4 production in cattle and goats in this study fed the COS diet, which was observed to additionally impact other bacterial metabolism and ruminal fermentation parameters, such as increased propionate and decreased butyrate and acetate levels (Table 1) [5]. Thus, while we cannot unequivocally rule out that diets supplemented with lipids at higher levels could be directly affecting the abundance of protozoa, we believe it is also possible that corn oil is influencing other fiber-degrading and hydrogen-producing species in the rumen, hence effecting ruminal fermentation and methane emissions [87, 92]. Nonetheless, validation of these results under different conditions and extended sampling is necessary to verify our findings.

Increases in dietary starch for ruminants is known to stimulate the propionate and succinate pathways of starch-degrading bacteria, which due to their net incorporation of metabolic hydrogen [H] represent a [H] sink in rumen fermentation besides hydrogenotrophic methanogenesis [82, 93]. Diets rich in starch are more fermentable in the rumen, which can decrease the ruminal pH to levels that can inhibit methanogenic archaea and fiber-degrading bacterial species [94, 95]. Yet, lowered pH levels in the rumen can also lead to clinical (or sub-clinical in most production scenarios) ruminal acidosis [96, 97]. Hence, high concentrate diets, which increase production of propionate at the expense of methane, does not necessarily opt for a viable methane mitigation strategy in the long term. Our results suggests that decreased methanogenesis in COS-fed animals is likely due to a decrease in available hydrogen and/or decrease in pH levels, which we predict is caused by the lower levels of hydrogen-producing protozoa and bacteria as well as the metabolism of dominant bacterial populations that likely do not produce exogenous hydrogen due to their own [H]-utilizing succinate and propionate metabolism. This prediction was partially supported by the detection of several putative hydrogenases believed to use molecular hydrogen in putative fumarate-reducing Aminobacteriaceae populations that were dominant in the COS-fed goat rumen. In particular, two [FeFe]-hydrogenases were detected in COS-fed animals from RUDB_G023 (ORF_4456741_13) and RUDB_G016 (ORF_1143158_19) that were annotated as NADP-reducing hydrogenase subunits (EC:1.12.1.3) (Fig. 2d, Supplementary Table S6).

To further test our hypothesis that protozoal species play a central role in the rumen ecosystem, we explored additional metagenome-centric metaproteomic datasets originating from an independent feeding experiment performed in Ireland on 60 Holstein Friesian bulls [65]. In brief, these bulls were subjected to the same ratio of concentrate and grass silage at either an ad libitum or restricted feeding regime in a compensatory growth model detailed in Keogh et al. [65]. We applied the same strategy as for the described Holstein dairy cattle and alpine goats to resolve the metaproteomic dataset for a subset of 15 animals (7 restricted and 8 ad libitum) against 781 reconstructed sample-specific MAGs (RUDB-HF), which were supplemented with the genome of En. caudatum, 18 protozoal SAGs [15] as well as genomes of available anaerobic fungi. This collection of microbial prokaryote and eukaryote genomes was then used as a sequence database for the generated protein spectra. Consistent with our previous observation in dairy cattle, a substantial proportion of the detected proteins were affiliated to Epidinium (Ep. cattanei and Ep. caudatum), Entodinium (En. caudatum and En. bursa) and Isotricha spp. providing further support that protozoa are an important and metabolically active contributor to the rumen microbiome (Fig. 5, Supplementary Table S8). The protein quantification (measured as sum of LFQ intensities affiliated to each MAG/genome, averaged for each diet) was also higher in the rumen sample from bulls on the restricted diet, which likely had less starch available compared to the ad libitum group and a higher retention time (Fig. 5a, Supplementary Table S8). A previously published 16S rRNA amplicon investigation of the phylogenetic differences between the rumen microbiomes of these two diet groups highlighted an increase in Succinivibrionaceae in the starch-rich ad libitum diet [66]. Our metaproteomic analysis confirmed a significantly higher (unpaired t-test, p value < 0.05) proteomic detection of several Succinivibrionaceae-MAGs under the ad libitum group (Supplementary Table S8), accompanied with a reduced acetate:propionate ratio in the rumen, which is often associated with increased feed efficiency and reduced production of methane [66]. These observations largely mirror the dominance of Succinivibrionaceae-MAGs in the dairy cattle and goats fed the COS diet, further strengthening our hypothesis that Entodinium and Epidinium spp do not metabolically respond to increases in available starch in the host animals’ diets and have other roles than being a primary starch degrader.

A total of 60 beef cattle were subjected to two dietary contrasting condition: 30 animals with ad libitum feeding and 30 subjected to 125 days of feed restriction. Dietary components in both treatments consisted of 70% concentrate, and 30% grass silage, with the concentrate containing rolled barley 72.5%, soya 22.5%, molasses 3% and calf mineral 2%. Rolled barley is high in energy and starch content (~50%). Metaproteomes for a subset of 15 animals (7 restricted and 8 ad libitum) were analysed against a database that combined 19 protozoal SAGs/genome (a) and 781 MAGs and isolate genomes (b), including fungal representatives. Protozoal proteomes are presented separately as the scale of their protein quantification values were ~10× larger. Protein quantification values were calculated by considering both the number of proteins detected per MAG/SAG/genome and their abundance; we averaged LFQ intensities for each detected protein across biological replicates for each dietary condition (ad libitum: green or restricted: orange), which were subsequently summed for all detected proteins per MAG/SAG/genome. Similar to our observations in Holstein dairy heifers and Alpine goats (Fig. 2), the proteomes of protozoal species were a major fraction of the total rumen metaproteome, however they were substantially reduced in dietary conditions where starch content was higher. Further supporting our results in dairy cattle, we observed in high-starch conditions an increase in protein detection for bacterial populations affiliated to the Succinivibrionaceae. c Volcano plot indicating different rumen microbiome proteins that displayed both large magnitude of fold-changes in LFQ intensities (x axis) and high statistical significance (−log10 of nominal p values using an unpaired t-test, y axis). Dashed horizontal line denotes a p value <0.05 cut-off. Volcano plots showed an increase in protein detection for bacterial populations (greenblue) in animals fed ad libitum compared to protozoal populations (purple), which were detected at higher LFQ intensities in animals under restricted conditions. LFQ intensities used to create (a, b) can be found in Supplementary Table S8.

To assess whether similar patterns of diet-induced microbiome shifts were consistent across the two geographically distinct datasets used in this study, we compared ratios of protein abundance level for commonly observed taxonomic families. Spearman's rank correlation coefficients were calculated for the comparisons, where two tests were performed: correlation between common taxonomic families (rarely observed families with <2 MAGs with detected proteins across both conditions were excluded) in the rumen of beef (RUDB_HF) and dairy (RUDB_C) cattle as well as beef cattle (RUDB_HF) and alpine goats (RUDB_G), respectively. Changes in protein abundance levels from a total of 16 taxonomic families were found in common between the rumen microbiomes of beef and dairy cattle (significant correlation: p value < 0.05 and correlation coefficient rho of 0.68), further strengthening our hypothesis that similar structural rumen microbiome shifts occur across independent datasets of cattle subjected to diets of high and low starch (Fig. S1, Supplementary Table S9). However, unsurprisingly no significant correlations were observed between the rumen microbiomes of goats (RUDB_G) and beef cattle (RUDB_HF) (p value >0.05 and correlation coefficient rho of 0.19, Supplementary Table S9, Supplementary Fig. S1), despite similar trends in protein abundance levels across taxa being observed (Figs. 2, 5, and Supplementary Fig. S1).

In conclusion, by using a (meta)genome-centric metaproteomics approach we primarily investigated the role of the rumen protozoa in the rumen microbiome of beef and dairy cattle as well as dairy goats that were subjected to varying dietary conditions. We showed that the proteomes of core entodinomorph and holotrich genera such as Entodinium, Epidinium, Dasytricha and Isotricha constitute a substantial fraction of the recovered rumen microbial proteome, which supports previous 16S/18S rRNA gene-based rumen census data that have highlighted their global dominance across a plethora of ruminant species. In both animal hosts, protozoal CAZymes targeting lignocellulose were detected though more frequently and at higher quantification levels in animals fed the CTL diet compared to the COS supplemented diet, presumably due to lower cell density in the later (Table 1). Proteins identified as En. caudatum, Epidinium, Isotricha and Dasytricha spp. were surprisingly detected at lower levels in animals that were fed increased levels of wheat starch, despite their reputable starch-metabolizing capabilities (Figs. 2–5). We hypothesize that protozoa were being out competed by Gram-negative bacterial species (e.g., Succinivibrionaceae in cattle and/or Aminobacteriaceae in goats) in the COS-fed animals, which were possibly resistant to protozoal engulfment [77,78,79] and/or lower pH levels, creating sub-optimal conditions for starch-degrading protozoal populations. We also observed lower detection of protozoal [FeFe]-hydrogenases in COS-fed animals that measured a ~25–30% lower CH4 yield at the time of sampling in this study (prior to feeding), further suggesting the specific conditions that enable succinate- and propionate-producing populations to flourish subsequently exert an impact on hydrogen and methane metabolisms in the rumen microbiome. For animals that were additionally fed corn oil, it is acknowledged that diets supplemented with lipids can have effects on the decreased abundance of protozoa, however our study did not provide evidence of a direct effect. While much work is still needed to confirm our abovementioned hypotheses, our integrated metaproteomics approaches have demonstrated the future importance of including eukaryote populations for accurate and meaningful analyses of the rumen microbiome and its impact on GHG mitigation strategies and host productivity traits.

Raw shotgun metagenomic data has been deposited in the National Center for Biotechnology Sequence Read Archive (NCBI-SRA) under accessions numbers SRR19524239 to SRR19524270 with links to BioProject accession number PRJNA844951. All annotated prokaryote MAGs are available publicly at https://doi.org/10.6084/m9.figshare.20079461 and https://doi.org/10.6084/m9.figshare.20066972. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [98] partner repository with the dataset identifiers PXD040467, PXD040454 and PXD040349. NCBI accession numbers for the protozal SAGs used in this study are as follows: Dasytricha ruminantium SAG3 (JAJJKK000000000), Diplodinium dentatum SAGT1 (JAJJLX000000000), Diplodinium flabellum SAG1 (JAJJLB000000000), Enoploplastron triloricatum SAGT1 (JAJNAB000000000), Entodinium bursa SAG3 (JAJJKR000000000), Epidinium cattanei SAG3 (JAJJKG000000000), Epidinium caudatum SAG1 (JAJMZR000000000), Eremoplastron rostratum SAG2 (JAJMZO000000000), Isotricha intestinalis SAGT2 (JAJJLT000000000), Isotricha prostoma SAG3 (JAJMZS000000000), Isotricha sp. YL-2021a SAG1 (JAJJLH000000000), Isotricha sp. YL-2021b SAG3 (JAJJKW000000000), Metadinium minomm SAG1 (JAJMZT000000000), Ophryoscolex caudatus SAGT3 (JAJJLZ000000000), Ostracodinium dentatum SAG1 (JAJJKF000000000), Ostracodinium gracile SAG1 (JAJJKZ000000000), Polyplastron multivesiculatum SAGT3 (JAJJLS000000000).

Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70:567–90.

Article CAS PubMed Google Scholar

Wallace RJ, Rooke JA, McKain N, Duthie CA, Hyslop JJ, Ross DW, et al. The rumen microbial metagenome associated with high methane production in cattle. BMC Genom. 2015;16:839.

Article Google Scholar

Jami E, White BA, Mizrahi I. Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS ONE. 2014;9:e85423.

Article PubMed PubMed Central Google Scholar

McCann JC, Luan S, Cardoso FC, Derakhshani H, Khafipour E, Loor JJ. Induction of Subacute Ruminal Acidosis Affects the Ruminal Microbiome and Epithelium. Front Microbiol. 2016;7:701.

Article PubMed PubMed Central Google Scholar

Martin C, Coppa M, Fougère H, Bougouin A, Baumont R, Eugène M, et al. Diets supplemented with corn oil and wheat starch, marine algae, or hydrogenated palm oil modulate methane emissions similarly in dairy goats and cows, but not feeding behavior. Anim Feed Sci Technol. 2021;272:114783.

Article CAS Google Scholar

Seshadri R, Leahy SC, Attwood GT, Teh KH, Lambie SC, Cookson AL, et al. Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection. Nat Biotechnol. 2018;36:359–67.

Article CAS PubMed PubMed Central Google Scholar

Stewart RD, Auffret MD, Warr A, Walker AW, Roehe R, Watson M. Compendium of 4,941 rumen metagenome-assembled genomes for rumen microbiome biology and enzyme discovery. Nat Biotechnol. 2019;37:953–61.

Article CAS PubMed PubMed Central Google Scholar

Xie F, Jin W, Si H, Yuan Y, Tao Y, Liu J, et al. An integrated gene catalog and over 10,000 metagenome-assembled genomes from the gastrointestinal microbiome of ruminants. Microbiome. 2021;9:137.

Article CAS PubMed PubMed Central Google Scholar

Williams AG, Coleman GS. The rumen microbial ecosystem. 2nd ed. Dordrecht: Springer; 1997.

Google Scholar

Lin C, Raskin L, Stahl DA. Microbial community structure in gastrointestinal tracts of domestic animals: comparative analyses using rRNA-targeted oligonucleotide probes. FEMS Microbiol Ecol. 1997;22:281–94.

Article CAS Google Scholar

Hagen LH, Brooke CG, Shaw CA, Norbeck AD, Piao H, Arntzen MO, et al. Proteome specialization of anaerobic fungi during ruminal degradation of recalcitrant plant fiber. ISME J. 2021;15:421–34.

Article CAS PubMed Google Scholar

Saye LMG, Navaratna TA, Chong JPJ, O’Malley MA, Theodorou MK, Reilly M. The anaerobic fungi: challenges and opportunities for industrial lignocellulosic biofuel production. Microorganisms. 2021;9:694.

Article CAS PubMed PubMed Central Google Scholar

Stabel M, Hanafy RA, Schweitzer T, Greif M, Aliyu H, Flad V, et al. Aestipascuomyces dupliciliberans gen. nov, sp. nov., the first cultured representative of the uncultured SK4 clade from aoudad sheep and alpaca. Microorganisms. 2020;8:1734.

Article CAS PubMed PubMed Central Google Scholar

Park T, Wijeratne S, Meulia T, Firkins JL, Yu Z. The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment. Genomics. 2021;113:1416–27.

Article CAS PubMed Google Scholar

Li Z, Wang X, Zhang Y, Yu Z, Zhang T, Dai X, et al. Genomic insights into the phylogeny and biomass-degrading enzymes of rumen ciliates. ISME J. 2022;16:2775–87.

Article CAS PubMed Google Scholar

Henderson G, Cox F, Ganesh S, Jonker A, Young W, Global Rumen Census C, et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5:14567.

Article CAS PubMed PubMed Central Google Scholar

Ranilla MJ, Jouany JP, Morgavi DP. Methane production and substrate degradation by rumen microbial communities containing single protozoal species in vitro. Lett Appl Microbiol. 2007;45:675–80.

Article CAS PubMed Google Scholar

Yarlett N, Coleman GS, Williams AG, Lloyd D. Hydrogenosomes in known species of rumen entodiniomorphid protozoa. FEMS Microbiol Lett. 1984;21:15–9.

Article CAS Google Scholar

Park T, Meulia T, Firkins JL, Yu Z. Inhibition of the Rumen Ciliate Entodinium caudatum by Antibiotics. Front Microbiol. 2017;8:1189.

Article PubMed PubMed Central Google Scholar

Belanche A, de la Fuente G, Newbold CJ. Study of methanogen communities associated with different rumen protozoal populations. FEMS Microbiol Ecol. 2014;90:663–77.

Article CAS PubMed Google Scholar

Hillman D, Lloyd K, Yarlett N, Williams AG. Hydrogen production by rumen holotrich protozoa: effects of oxygen and implications for metabolic control by in situ conditions. J Protozool. 1989;36:205–13.

Article PubMed Google Scholar

Fougere H, Bernard L. Effect of diets supplemented with starch and corn oil, marine algae, or hydrogenated palm oil on mammary lipogenic gene expression in cows and goats: A comparative study. J Dairy Sci. 2019;102:768–79.

Article CAS PubMed Google Scholar

Fougere H, Delavaud C, Bernard L. Diets supplemented with starch and corn oil, marine algae, or hydrogenated palm oil differentially modulate milk fat secretion and composition in cows and goats: A comparative study. J Dairy Sci. 2018;101:8429–45.

Article CAS PubMed Google Scholar

Solomon KV, Haitjema CH, Henske JK, Gilmore SP, Borges-Rivera D, Lipzen A, et al. Early-branching gut fungi possess a large, comprehensive array of biomass-degrading enzymes. Science. 2016;351:1192–5.

Article CAS PubMed PubMed Central Google Scholar

Haitjema CH, Gilmore SP, Henske JK, Solomon KV, de Groot R, Kuo A, et al. A parts list for fungal cellulosomes revealed by comparative genomics. Nat Microbiol. 2017;2:17087.

Article CAS PubMed Google Scholar

Youssef NH, Couger MB, Struchtemeyer CG, Liggenstoffer AS, Prade RA, Najar FZ, et al. The genome of the anaerobic fungus Orpinomyces sp. strain C1A reveals the unique evolutionary history of a remarkable plant biomass degrader. Appl Environ Microbiol. 2013;79:4620–34.

Article CAS PubMed PubMed Central Google Scholar

Hooker CA, Hillman ET, Overton JC, Ortiz-Velez A, Schacht M, Hunnicutt A, et al. Hydrolysis of untreated lignocellulosic feedstock is independent of S-lignin composition in newly classified anaerobic fungal isolate, Piromyces sp. UH3-1. Biotechnol Biofuels. 2018;11:293.

Article CAS PubMed PubMed Central Google Scholar

Brown JL, Swift CL, Mondo SJ, Seppala S, Salamov A, Singan V, et al. Cocultivation of the anaerobic fungus Caecomyces churrovis with Methanobacterium bryantii enhances transcription of carbohydrate binding modules, dockerins, and pyruvate formate lyases on specific substrates. Biotechnol Biofuels. 2021;14:234.

Article CAS PubMed PubMed Central Google Scholar

Wilken SE, Monk JM, Leggieri PA, Lawson CE, Lankiewicz TS, Seppälä S, et al. Experimentally validated reconstruction and analysis of a genome-scale metabolic model of an anaerobic Neocallimastigomycota fungus. mSystems. 2021;6:e00002–21.

Article CAS PubMed PubMed Central Google Scholar

Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques. 2004;36:808–12.

Article CAS PubMed Google Scholar

Morgavi DP, Jouany J-P, Martin C. Changes in methane emission and rumen fermentation parameters induced by refaunation in sheep. Aust J Exp Agric. 2008;48:69–72.

Article CAS Google Scholar

Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

Article CAS PubMed PubMed Central Google Scholar

Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.

Article CAS PubMed Google Scholar

Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.

Article CAS PubMed PubMed Central Google Scholar

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.

Article PubMed PubMed Central Google Scholar

Wu YW, Simmons BA, Singer SW. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics. 2016;32:605–7.

Article CAS PubMed Google Scholar

Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7:e7359.

Article PubMed PubMed Central Google Scholar

Alneberg J, Bjarnason BS, de Bruijn I, Schirmer M, Quick J, Ijaz UZ, et al. Binning metagenomic contigs by coverage and composition. Nat Methods. 2014;11:1144–6.

Article CAS PubMed Google Scholar

Sieber CMK, Probst AJ, Sharrar A, Thomas BC, Hess M, Tringe SG, et al. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat Microbiol. 2018;3:836–43.

Article CAS PubMed PubMed Central Google Scholar

Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36:W5–9.

Article CAS PubMed PubMed Central Google Scholar

Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.

Article CAS PubMed PubMed Central Google Scholar

Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil PA, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996–1004.

Article CAS PubMed Google Scholar

Parks DH, Chuvochina M, Chaumeil P-A, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol. 2020;38:1079–86.

Article CAS PubMed Google Scholar

Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35:725–31.

Article CAS PubMed PubMed Central Google Scholar

Shaffer M, Borton MA, McGivern BB, Zayed AA, La Rosa SL, Solden LM, et al. DRAM for distilling microbial metabolism to automate the curation of microbiome function. Nucleic Acids Res. 2020;48:8883–900.

Article CAS PubMed PubMed Central Google Scholar

Zhang H, Yohe T, Huang L, Entwistle S, Wu P, Yang Z, et al. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018;46:W95–101.

Article CAS PubMed PubMed Central Google Scholar

Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021;49:D412–9.

Article CAS PubMed Google Scholar

Suzek BE, Wang Y, Huang H, McGarvey PB, Wu CH, UniProt C. UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. Bioinformatics. 2015;31:926–32.

Article CAS PubMed Google Scholar

Rawlings ND, Barrett AJ, Thomas PD, Huang X, Bateman A, Finn RD. The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res. 2018;46:D624–32.

Article CAS PubMed Google Scholar

Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36:2251–2.

Article CAS PubMed Google Scholar

Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428:726–31.

Article CAS PubMed Google Scholar

Kanehisa M, Sato Y, Kawashima M. KEGG mapping tools for uncovering hidden features in biological data. Protein Sci. 2022;31:47–53.

Article CAS PubMed Google Scholar

Nordberg H, Cantor M, Dusheyko S, Hua S, Poliakov A, Shabalov I, et al. The genome portal of the Department of Energy Joint Genome Institute: 2014 updates. Nucleic Acids Res. 2014;42:D26–31.

Article CAS PubMed Google Scholar

Peng X, Wilken SE, Lankiewicz TS, Gilmore SP, Brown JL, Henske JK, et al. Genomic and functional analyses of fungal and bacterial consortia that enable lignocellulose breakdown in goat gut microbiomes. Nat Microbiol. 2021;6:499–511.

Article CAS PubMed PubMed Central Google Scholar

Mićić M, Whyte JD, Karsten V. in Sample Preparation Techniques for Soil, Plant, and Animal Samples. Micic M, editor. New Yok: Springer New York; 2016. p. 99–116.

Michalak L, Gaby JC, Lagos L, La Rosa SL, Hvidsten TR, Tetard-Jones C, et al. Microbiota-directed fibre activates both targeted and secondary metabolic shifts in the distal gut. Nat Commun. 2020;11:5773.

Article CAS PubMed PubMed Central Google Scholar

Kong AT, Leprevost FV, Avtonomov DM, Mellacheruvu D, Nesvizhskii AI. MSFragger: ultrafast and comprehensive peptide identification in mass spectrometry–based proteomics. Nat Methods. 2017;14:513–20.

Article CAS PubMed PubMed Central Google Scholar

da Veiga Leprevost F, Haynes SE, Avtonomov DM, Chang H-Y, Shanmugam AK, Mellacheruvu D, et al. Philosopher: a versatile toolkit for shotgun proteomics data analysis. Nat Methods. 2020;17:869–70.

Article PubMed PubMed Central Google Scholar

Yu F, Haynes SE, Nesvizhskii AI. IonQuant enables accurate and sensitive label-free quantification with FDR-controlled match-between-runs. Mol Cell Proteom. 2021;20:100077.

Article CAS Google Scholar

Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods. 2016;13:731–40.

Article CAS PubMed Google Scholar

Kleiner M, Thorson E, Sharp CE, Dong X, Liu D, Li C, et al. Assessing species biomass contributions in microbial communities via metaproteomics. Nat Commun. 2017;8:1558.

Article PubMed PubMed Central Google Scholar

Wickham, H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. 2016. https://ggplot2.tidyverse.org.

R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2021. https://www.R-project.org/.

Federico A, Monti S. hypeR: an R package for geneset enrichment workflows. Bioinformatics. 2020;36:1307–8.

Article CAS PubMed Google Scholar

Keogh K, Waters SM, Kelly AK, Kenny DA. Feed restriction and subsequent realimentation in Holstein Friesian bulls: I. Effect on animal performance; muscle, fat, and linear body measurements; and slaughter characteristics. J Anim Sci. 2015;93:3578–89.

Article CAS PubMed Google Scholar

McCabe MS, Cormican P, Keogh K, O’Connor A, O’Hara E, Palladino RA, et al. Illumina MiSeq phylogenetic amplicon sequencing shows a large reduction of an uncharacterised Succinivibrionaceae and an increase of the Methanobrevibacter gottschalkii clade in feed restricted cattle. PLOS ONE. 2015;10:e0133234.

Article PubMed PubMed Central Google Scholar

Wessel D, Flugge UI. A method for the quantitative recovery of protein in dilute solution in the presence of detergents and lipids. Anal Biochem. 1984;138:141–3.

Article CAS PubMed Google Scholar

Zougman A, Selby PJ, Banks RE. Suspension trapping (STrap) sample preparation method for bottom-up proteomics analysis. Proteomics. 2014;14:1006–000.

Article CAS PubMed Google Scholar

Harper JW, Bennett EJ. Proteome complexity and the forces that drive proteome imbalance. Nature. 2016;537:328–38.

Article CAS PubMed PubMed Central Google Scholar

Williams AG, Coleman GS. The Rumen Protozoa. New York: Springer; 1992.

Book Google Scholar

Williams CL, Thomas BJ, McEwan NR, Stevens PR, Creevey CJ, Huws SA. Rumen protozoa play a significant role in fungal predation and plant carbohydrate breakdown. Front Microbiol. 2020;11:720.

Takenaka A, Tajima K, Mitsumori M, Kajikawa H. Fiber digestion by rumen ciliate protozoa. Microbes Environ. 2004;19:203–10.

Article Google Scholar

Coleman GS, Sandford DC. The engulfment and digestion of mixed rumen bacteria and individual bacterial species by single and mixed species of rumen ciliate protozoa grown in vivo. J Agric Sci. 1979;92:729–42.

Article Google Scholar

Wang L, Abu-Doleh A, Plank J, Catalyurek UV, Firkins JL, Yu Z. The transcriptome of the rumen ciliate Entodinium caudatum reveals some of its metabolic features. BMC Genom. 2019;20:1008.

Article CAS Google Scholar

Belzecki G, McEwan NR, Kowalik B, Michalowski T, Miltko R. Effect of Entodinium caudatum on starch intake and glycogen formation by Eudiplodinium maggii in the rumen and reticulum. Eur J Protistol. 2017;57:38–49.

Article PubMed Google Scholar

Allen MS. Relationship between fermentation acid production in the rumen and the requirement for physically effective fiber. J Dairy Sci. 1997;80:1447–62.

Article CAS PubMed Google Scholar

Solomon R, Wein T, Levy B, Eshed S, Dror R, Reiss V, et al. Protozoa populations are ecosystem engineers that shape prokaryotic community structure and function of the rumen microbial ecosystem. ISME J. 2022;16:1187–97.

Article PubMed Google Scholar

Park T, Yu Z. Do ruminal ciliates select their preys and prokaryotic symbionts? Front Microbiol. 2018;9:1710.

Article PubMed PubMed Central Google Scholar

Gutierrez J, Davis RE. Bacterial ingestion by the rumen ciliates Entodinium and Diplodinium. J Protozool. 1959;6:222–6.

Article Google Scholar

Bauman DE, Davis CL, Bucholtz HF. Propionate production in the rumen of cows fed either a control or high-grain, low-fiber diet. J Dairy Sci. 1971;54:1282–7.

Article CAS PubMed Google Scholar

Jiao HP, Dale AJ, Carson AF, Murray S, Gordon AW, Ferris CP. Effect of concentrate feed level on methane emissions from grazing dairy cows. J Dairy Sci. 2014;97:7043–53.

Article CAS PubMed Google Scholar

Martin C, Morgavi DP, Doreau M. Methane mitigation in ruminants: from microbe to the farm scale. Animal. 2010;4:351–65.

Article CAS PubMed Google Scholar

Finlay BJ, Esteban G, Clarke KJ, Williams AG, Embley TM, Hirt RP. Some rumen ciliates have endosymbiotic methanogens. FEMS Microbiol Lett. 1994;117:157–61.

Article CAS PubMed Google Scholar

Newbold CJ, Lassalas B, Jouany JP. The importance of methanogens associated with ciliate protozoa in ruminal methane production in vitro. Lett Appl Microbiol. 1995;21:230–4.

Article CAS PubMed Google Scholar

Newbold CJ, de la Fuente G, Belanche A, Ramos-Morales E, McEwan NR. The role of ciliate protozoa in the rumen. Front Microbiol. 2015;6:1313.

Article PubMed PubMed Central Google Scholar

Hristov AN, Oh J, Firkins JL, Dijkstra J, Kebreab E, Waghorn G, et al. Special topics–Mitigation of methane and nitrous oxide emissions from animal operations: I. A review of enteric methane mitigation options. J Anim Sci. 2013;91:5045–69.

Article CAS PubMed Google Scholar

Jenkins TC, Wallace RJ, Moate PJ, Mosley EE. Board-invited review: Recent advances in biohydrogenation of unsaturated fatty acids within the rumen microbial ecosystem. J Anim Sci. 2008;86:397–412.

Article CAS PubMed Google Scholar

Zhang XM, Medrano RF, Wang M, Beauchemin KA, Ma ZY, Wang R, et al. Corn oil supplementation enhances hydrogen use for biohydrogenation, inhibits methanogenesis, and alters fermentation pathways and the microbial community in the rumen of goats. J Anim Sci. 2019;97:4999–5008.

Article PubMed PubMed Central Google Scholar

Guyader J, Eugène M, Nozière P, Morgavi DP, Doreau M, Martin C. Influence of rumen protozoa on methane emission in ruminants: a meta-analysis approach. Animal. 2014;8:1816–25.

Article CAS PubMed Google Scholar

Niu P, Schwarm A, Bonesmo H, Kidane A, Aspeholen Aby B, Storlien TM, et al. A basic model to predict enteric methane emission from dairy cows and Its application to update operational models for the national inventory in Norway. Animals. 2021;11:1891.

Article PubMed PubMed Central Google Scholar

McAllister TA, Cheng K-J, Okine EK, Mathison GW. Dietary, environmental and microbiological aspects of methane production in ruminants. Can J Anim Sci. 1996;76:231–43.

Article CAS Google Scholar

Gruninger RJ, Zhang XM, Smith ML, Kung L Jr., Vyas D, McGinn SM, et al. Application of 3-nitrooxypropanol and canola oil to mitigate enteric methane emissions of beef cattle results in distinctly different effects on the rumen microbial community. Anim Microbiome. 2022;4:35.

Article CAS PubMed PubMed Central Google Scholar

Popova M, Martin C, Eugène M, Mialon MM, Doreau M, Morgavi DP. Effect of fibre- and starch-rich finishing diets on methanogenic Archaea diversity and activity in the rumen of feedlot bulls. Anim Feed Sci Technol. 2011;166-167:113–21.

Article CAS Google Scholar

Van Kessel JAS, Russell JB. The effect of pH on ruminal methanogenesis. FEMS Microbiol Ecol. 1996;20:205–10.

Article Google Scholar

Russell JB, Muck RE, Weimer PJ. Quantitative analysis of cellulose degradation and growth of cellulolytic bacteria in the rumen. FEMS Microbiol Ecol. 2009;67:183–97.

Article CAS PubMed Google Scholar

Russell JB. The importance of pH in the regulation of ruminal acetate to propionate ratio and methane production in vitro. J Dairy Sci. 1998;81:3222–30.

Article CAS PubMed Google Scholar

Kleen JL, Hooijer GA, Rehage J, Noordhuizen JP. Subacute ruminal acidosis (SARA): a review. J Vet Med A Physiol Pathol Clin Med. 2003;50:406–14.

Article CAS PubMed Google Scholar

Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47:D442–50.

Article CAS PubMed Google Scholar

Download references

PBP and TOA are grateful for support from The Research Council of Norway (FRIPRO program, PBP: 250479), as well as the European Research Commission Starting Grant Fellowship (awarded to PBP: 336355—MicroDE), and the Novo Nordisk Foundation (awarded to PBP: 0054575—SuPAcow). LHH was supported by The Research Council of Norway (FRIPRO program, LHH: 302639—SeaCow), while MØA was supported by the Novo Nordisk Foundation, project No. NNF20OC006131. The experimental trial was financed by APIS-GENE (Paris, France) as part of the NutriLip project. Gas emission measurements were funded by UMR 1213 Herbivores (INRAE, Saint-Genès-Champanelle, France). The sequencing service was provided by the Norwegian Sequencing Center (www.sequencing.uio.no), a national technology platform hosted by the University of Oslo and supported by the "Functional Genomics" and "Infrastructure" programs of the Research Council of Norway and the Southeastern Regional Health Authorities. The authors acknowledge the Orion High Performance Computing Center at the Norwegian University of Life Sciences and Sigma2—the National Infrastructure for High Performance Computing and Data Storage in Norway for providing computational resources that have contributed to meta-omics computations reported in this paper. Mass spectrometry-based proteomic analyses were performed by The MS and Proteomics Core Facility, Norwegian University of Life Sciences (NMBU). This facility is a member of the National Network of Advanced Proteomics Infrastructure (NAPI), which is funded by the Research Council of Norway INFRASTRUKTUR-program (project number: 295910).

Department of Animal and Aquacultural Sciences, Faculty of Biosciences, Norwegian University of Life Sciences, Ås, Norway

Thea O. Andersen, Ianina Altshuler & Phillip B. Pope

Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway

Arturo Vera-Ponce de León, Juline M. Walter, Torgeir R. Hvidsten, Magnus Ø. Arntzen, Live H. Hagen & Phillip B. Pope

Teagasc, Animal and Bioscience Research Department, Animal and Grassland Research and Innovation Centre, Teagasc, Grange, Dunsany, County, Meath, Ireland

Emily McGovern, Kate Keogh & Sinead M. Waters

INRAE, VetAgro Sup, UMR Herbivores, Université Clermont Auvergne, Saint-Genes-Champanelle, France

Cécile Martin, Laurence Bernard, Diego P. Morgavi & Milka Popova

Department of Animal Science and Technology, Chung-Ang University, Anseong-si, Gyeonggi-do, Republic of Korea

Tansol Park

Department of Animal Sciences, The Ohio State University, Columbus, OH, USA

Tansol Park, Jeffrey L. Firkins & Zhongtang Yu

Center for Ruminant Genetics and Evolution, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang, China

Zongjun Li & Yu Jiang

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

TOA, LHH, MØA and PBP conceived the study, performed the primary analysis of the data and wrote the paper (with input from all authors). MP, CM, LB and DPM designed, performed and analyzed the INRAE animal experiments. SMW, KK and EM designed, performed and analyzed the TEAGASC animal experiments. TOA, LHH, MØA, AVPDL, IA, JMW and TRH generated the metaproteomic data and/or contributed to the data analyses. TP, ZY, JLF, YJ and ZL generated and analyzed protozoal genomic data.

Correspondence to Phillip B. Pope.

The authors declare no competing interests.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Andersen, T.O., Altshuler, I., Vera-Ponce de León, A. et al. Metabolic influence of core ciliates within the rumen microbiome. ISME J (2023). https://doi.org/10.1038/s41396-023-01407-y

Download citation

Received: 29 June 2022

Revised: 29 March 2023

Accepted: 30 March 2023

Published: 11 May 2023

DOI: https://doi.org/10.1038/s41396-023-01407-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative