Susanne Künzel1,2†
Timur Yergaliyev1,2†
Katharina J. Wild1,2Hanna Philippi1,2Asta H. Petursdottir3Helga Gunnlaugsdottir3,4Chris K. Reynolds5David J. Humphries5
Amélia Camarinha-Silva1,2
Markus Rodehutscord1,2*- 1Institute of Animal Science, University of Hohenheim, Stuttgart, Germany
- 2Hohenheim Center for Livestock Microbiome Research, University of Hohenheim, Stuttgart, Germany
- 3Matís, Reykjavík, Iceland
- 4Faculty of Food Science and Nutrition, School of Health Sciences, University of Iceland, Reykjavík, Iceland
- 5School of Agriculture, Policy and Development, University of Reading, Reading, United Kingdom
This study aimed to investigate the effects of two brown Icelandic seaweed samples (Ascophyllum nodosum and Fucus vesiculosus) on in vitro methane production, nutrient degradation, and microbiota composition. A total mixed ration (TMR) was incubated alone as control or together with each seaweed at two inclusion levels (2.5 and 5.0% on a dry matter basis) in a long-term rumen simulation technique (Rusitec) experiment. The incubation period lasted 14 days, with 7 days of adaptation and sampling. The methane concentration of total gas produced was decreased at the 5% inclusion level of A. nodosum and F. vesiculosus by 8.9 and 3.6%, respectively (P < 0.001). The total gas production was reduced by all seaweeds, with a greater reduction for the 5% seaweed inclusion level (P < 0.001). Feed nutrient degradation and the production of volatile fatty acids and ammonia in the effluent were also reduced, mostly with a bigger effect for the 5% inclusion level of both seaweeds, indicating a reduced overall fermentation (all P ≤ 0.001). Microbiota composition was analyzed by sequencing 16S rRNA amplicons from the rumen content of the donor cows, fermenter liquid and effluent at days 7 and 13, and feed residues at day 13. Relative abundances of the most abundant methanogens varied between the rumen fluid used for the start of incubation and the samples taken at day 7, as well as between days 7 and 13 in both fermenter liquid and effluent (P < 0.05). According to the differential abundance analysis with q2-ALDEx2, in effluent and fermenter liquid samples, archaeal and bacterial amplicon sequence variants were separated into two groups (P < 0.05). One was more abundant in samples taken from the treatment without seaweed supplementation, while the other one prevailed in seaweed supplemented treatments. This group also showed a dose-dependent response to seaweed inclusion, with a greater number of differentially abundant members between a 5% inclusion level and unsupplemented samples than between a 2.5% inclusion level and TMR. Although supplementation of both seaweeds at a 5% inclusion level decreased methane concentration in the total gas due to the high iodine content in the seaweeds tested, the application of practical feeding should be done with caution.
Introduction
The rapidly growing global population brings serious challenges to the food industry. While ruminant livestock are vital in sustaining food security by converting inedible plant matter into meat and dairy products, they significantly contribute to global methane (CH4) emissions, a potent greenhouse gas (Lassey, 2007). These CH4 emissions are mainly related to fermentation processes orchestrated by the rumen microbiome. Diet composition is one of the major factors influencing rumen microbial communities (de Menezes et al., 2011; Henderson et al., 2015) and is therefore linked to CH4 production by ruminants.
Feed production does not need arable land, which may help us cope with the increasing demand for animal products. Seaweed could be an alternative animal feed material that has already been used for thousands of years in coastal areas (Evans and Critchley, 2014). Some seaweed species can also affect ruminal CH4 production even at a low inclusion level in the feed. Pronounced effects were reported for Asparagopsis taxiformis grown in the Pacific Ocean, with CH4 formation reduced by up to 99% with seaweed inclusion ≤5% in vitro (Kinley et al., 2016; Machado et al., 2018; Roque et al., 2019). Numerous studies have demonstrated that the effect of seaweed supplementation on methanogenesis was associated with a modified rumen microbiome (Molina-Alcaide et al., 2017; Roque et al., 2019; Abbott et al., 2020). In the rumen, hydrogenotrophic methanogens produce CH4 by using hydrogen (H2) for carbon dioxide (CO2) reduction through the Wolfe cycle (Leahy et al., 2010; Thauer, 2012). Accordingly, the composition of archaeal methanogens and abundances of hydrogen-producing bacteria were identified as key factors associated with levels of CH4 emissions in ruminants (Tapio et al., 2017). Considering the relevance of H2 for CH4 production, manipulating H2 production or its utilization pathways through diet and microbiome composition is an approach that may provide new insights into the development of CH4 mitigation strategies.
Seaweeds are a very heterogeneous group of feeding substances, and their application in animal feeding has been restricted due to a lack of information about species-specific nutritive value. Several factors influence nutrient and bioactive compounds in seaweed, such as species, season, and site of harvesting (Paiva et al., 2018; Britton et al., 2021). To the best of our knowledge, most studies concerning seaweed in ruminant nutrition were conducted with species harvested in the Pacific Ocean. Since there is also a high variation of Atlantic seaweeds, there is a need for information about these seaweeds. The current study aimed to investigate the effect of two Icelandic seaweed samples on in vitro gas and methane production, nutrient degradation, and microbiota composition using a continuous long-term rumen simulation technique (Rusitec). Both seaweeds are endemic to Iceland and are also available in large quantities in other European countries. Our hypothesis was that two common and abundant North Atlantic seaweed species could reduce CH4 formation to the point where their abundance and ease of access would make them a viable option as a ruminant feed substance. A secondary objective was to investigate changes in microbiota composition over time in the in vitro system fed with or without seaweed.
Materials and Methods
Treatments
Two samples of seaweed naturally occurring in Iceland were studied for their effect on TMR formulated for cattle in a Rusitec system (Czerkawski and Breckenridge, 1977). The TMR was composed of 20% corn grain, 20% soybean meal, 40% corn silage, and 20% grass silage and served as the control treatment. One seaweed was Ascophyllum nodosum (AN), harvested in August 2018, and the other was Fucus vesiculosus (FV), harvested in June 2019. Both seaweeds were used at 2.5% (AN2.5 and FV2.5) and 5% (AN5 and FV5) inclusion in the TMR on a dry matter (DM) basis in exchange for TMR. The seaweeds and all ingredients of the TMR were dried and ground to pass a 1 mm screen. The analyzed nutrient composition, bromoform, and total phenolic content (TPC) of the seaweeds and TMR and the calculated nutrient composition of the experimental treatments are shown in Table 1.
Table 1. Nutrient composition, bromoform, total phenolic content of the used seaweeds, and five experimental diets.
Rumen Content and Synthetic Saliva
Rumen content was collected from three rumen-fistulated non-lactating Jersey cows before the morning feeding. Animals had free access to water and a diet composed of 33% corn silage, 33% grass silage, 23% grass hay, 10% barley straw, and 1% mineral mixture (on a DM basis). During the daytime, the cows were kept on pasture. Two liters of rumen fluid were taken from each cow into prewarmed thermos flasks, comprised of 1 L pumped from the liquid phase and 1 L squeezed out from the solid phase. Additionally, 200 g of squeezed solid phase from each cow was transferred into prewarmed, isolated containers. Rumen fluid was strained through two layers of cheesecloth and mixed with equal parts from the donor animals. Subsequently, rumen fluid was mixed with a buffer solution (1:1), flushed with CO2, and stirred at 39°C until fermenters of the Rusitec were inoculated. The buffer solution was prepared according to the suggested composition of artificial saliva by McDougall (1948) with the addition of 15N enriched NH4Cl (0.0378 g/L; 104 mg 15N/g N) used for the calculation of microbial protein synthesis.
Rusitec System
Each run lasted 14 days (days 0–13), with 7 days of adaptation, followed by 7 days of sampling. Ten fermenters were arranged side by side, with five fermenters sharing one circulation thermostat (Lauda ECO E 4 S, Lauda-Königshofen, Germany) and one buffer pump (Ismatec IPC ISM 931, Wertheim, Germany) (Figure 1). Each of the five treatments was replicated in two fermenters for two runs each, resulting in four replications per treatment. Fermenters were randomized in each run. The circulation thermostat was used as a block design, and each treatment was used once in each block. Due to technical arrangements, an additional empty fermenter was connected to each circulation thermostat. The glass fermenters were kept at a constant temperature of 39°C with a heating jacket. The buffer pump continuously infused buffer solution into each fermenter at a daily rate of 75% of the fermenter’s capacity (950 ml) to stimulate salivation. Each fermenter contained a feed container with continuous vertical movement ensured by a lift motor to simulate rumen motility (10–12 strokes/min; Figure 2). The fluid effluent was separated from the fermentation gas in a glass cylinder and then collected in 1 L glass bottles. The glass cylinder and effluent bottles were placed in a 4°C tempered water bath. Gaseous effluent was passed through gas counters (BlueVCount, BlueSens gas sensor GmbH, Herten, Germany) to measure gas production via gas-tight tubes and was subsequently collected quantitatively in gas-tight five-layered plastic-aluminum bags (Dr.-Ing. Ritter Apparatebau GmbH & Co., KG, Bochum, Germany). The methane concentration of total gas production was measured using an infrared CH4 analyzer from the gas collected in the plastic-aluminum bags (PRONOVA Analysentechnik GmbH & Co., KG, Berlin, Germany).
Figure 1. Schematic diagram of the Rusitec set up with the 12 fermenters (a) connected to a lift motor (b), a buffer pump (c), and a gas-tight tube for the effluent (d). The effluent is separated in glass cylinders in a cooling water bath (e) into fluid effluent (f), which is collected in glass bottles in a cooling water bath (g), and gaseous effluent (h). The gaseous effluent is passed through a cold trap with an H2S absorber granulate (i), measured with gas counters (j), and sampled for methane analysis in plastic bags (k).
Figure 2. Schematic diagram of a Rusitec fermenter with a feed container (a) containing feed bags (b). The feed container is connected via a bar (c) to a lift motor. Buffer is consistently infused through a tube (d) connected to the buffer pump. The connection to the effluent (e) is closable to make it possible to remove the fermenter. The heating jacket (f) is connected via an inflow (g) and outflow (h) to a circulation thermostat.
Feed and solid rumen content were weighed into nylon bags (120 mm × 70 mm) with a 100 μm pore size and closed with cable ties. To start the system, one feed bag containing 15 g of the respective treatment and one feed bag containing 60 g of solid rumen content were put into the container of each fermenter. After 24 h, the bag with rumen content was replaced by another feed bag. Subsequently, each feed bag was removed after 48 h, rinsed with 50 (adaptation period) or 100 ml (sampling period) of buffer solution, squeezed moderately, and replaced by a new bag containing the specific treatment. The resulting liquid was returned to the respective fermenter.
Sampling and Chemical Analyses
The total gas production and CH4 concentration of total gas were measured in 24 h intervals during the sampling period. The temperature, pH, and redox potential of the fermenter fluid were measured daily in the fermenter before the feed bags were changed (SenTix ORP, WTW, Weilheim in Oberbayern, Germany; BlueLine 14 pH IDS, SI Analytics, Mainz, Germany). From days 7 to 13, a sample from the fermenter liquid phase (FL; 30 ml/d) was taken daily and pooled by the fermenter to obtain the fraction of liquid-associated microbes (LAM) by centrifugation according to Boguhn et al. (2006), with minor modifications. In brief, the suspension was centrifuged twice at 2,000 × g for 5 min at 4°C. Then, the supernatant was centrifuged three times at 15,000 × g for 15 min at 4°C. Afterward, samples were freeze-dried and pulverized with a ball mill (MM 400; Retsch GmbH, Haan, Germany). On day 13, solid-associated microbes (SAM) were separated from the feed residues (FR), as described by Ranilla and Carro (2003). Ammonia nitrogen (NH3-N), 15N abundance, and volatile fatty acids (VFA) were analyzed in the effluent (E), which was weighed and sampled daily (70 ml/d), pooled by the fermenter, and stored at −20°C until it was centrifuged for 15 min at 24,000 × g. The analysis of VFA was performed by vacuum distillation and gas chromatography (Hewlett-Packard 6890; Agilent, Waldbronn, Germany) measurement, as described by Wischer et al. (2013). NH3-N was analyzed using Kjeldahl steam distillation with phosphate buffer (Vapodest 50, Gerhardt, Königswinter, Germany). The daily production of VFA and NH3-N was calculated by relating the analyzed concentrations to the daily measured amount of E. The removed feed bags were dried for 24 h at 65°C, weighed, and pooled by the fermenter to determine nutrient degradation from days 7 to 12. Feed and FR were analyzed according to official methods in Germany (Verband Deutscher Landwirtschaftlicher Untersuchungs- und Forschungsanstalten, 1976) for DM (method 3.1), crude protein (CP; method 4.1.1), neutral detergent fiber on an ash-free basis (aNDFom; method 6.5.1), and acid detergent fiber on an ash-free basis (ADFom; method 6.5.2). Seaweeds and TMR were additionally analyzed for crude ash (method 8.1) and ether extract (method 5.1.1). Degradation of nutrients was calculated as the difference between the input and output of each fermenter in relation to input and expressed as a percentage.
Bromoform was analyzed in seaweeds according to method 8260B (EPA, 1996) by ALS Global (Miami, FL, United States). For the analysis of TPC, dried seaweeds were extracted in Milli-Q water (1:7 sample:water) for 1 h at room temperature (400 rpm). The TPC was determined on the sample supernatant according to the method by Singleton and Rossi (1965) adapted to the microplate format. For the analysis of heavy metals and minerals in the seaweeds, an ultraWAVE acid digestion system (Milestone Inc., Italy) and an Agilent 7900 quadrupole inductively coupled plasma mass spectrometer (Agilent Technologies, Singapore) were used according to NMKL (2007). Analyzed heavy metals and minerals are shown in Supplementary Table 1. Feed, FR, freeze-dried particle-free E, LAM, SAM, and NH4Cl were analyzed for 15N abundance using an elemental analyzer (EA, 1108; Carlo Erba Instruments, Biberach, Germany) combined with an isotope mass spectrometer (MS Finnigan MAT; ThermoQuest Italia S.p.A., Milan, Italy). SAM was additionally analyzed for total N.
Microbiota composition was analyzed in rumen solid phase (RSP) and rumen fluid (RF) from the liquid phase of each cow in both experimental runs (n = 6), and the mixture of rumen fluid from all cows and buffer solution, treated as RF replicates (n = 2) on day 0. On days 7 and 13, 1 ml of FL and E were taken from each fermenter for the analysis of microbiota composition (n = 40). Feed residues from both feed bags remaining for 24 and 48 h in the fermenter were sampled on day 13 (n = 40). Additionally, the liquid after SAM treatment of both bags was sampled on day 13 (n = 20). Samples were stored at −20°C immediately after collection.
Target Amplicon Sequencing
DNA was extracted with the commercial DNA extraction kit FastDNA™ Spin Kit for soil (MP Biomedicals, LLC, Solon, OH, United States), following the manufacturer’s instructions. DNA was quantified with a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and stored at −20°C. Two sequencing libraries were prepared as previously described to assess bacterial (V1–V2 region; Kaewtapee et al., 2017) and archaeal communities (Arch349 and Arch806 primers; Lee et al., 2012; Deusch et al., 2017). In brief, targeted regions of the 16S rRNA gene were amplified by a first polymerase chain reaction (PCR) using forward primers with 6-nt barcodes and a 2-nt linker attached and both primers were complementary to the Illumina adapters. Then, 1 μl of the resulting product was used as a template for a second PCR, this time using a reverse primer containing a specific sequence to multiplex and index primers. Obtained amplicons were checked by agarose gel electrophoresis and normalized using the SequalPrep Normalization Kit (Invitrogen Inc., Carlsbad, CA, United States). Sequencing was performed on the 250 bp paired-end Illumina NovaSeq 6000 platform.
Calculations
The estimated microbial protein synthesis (EMPS) was calculated in accordance with Boguhn et al. (2006) and Hildebrand et al. (2011). In brief, the daily outflow of microbial N from a fermenter (Nm; mg/d) was calculated as the sum of N originating from the LAM in the effluent (NLAM–E; mg/d) and SAM (NSAM; mg/d) fractions:
Therefore, NLAM–E was calculated with the sum of the daily input of 15N via buffer solution and feed (15Nin; μg/d), the sum of the daily output of 15N via FR and NH3-N (15Nout; μg/d), and the analyzed abundance of 15N in NLAM–E (15NLAM–E; μg/mg N):
The NSAM was calculated using the amount of N in FR (NFR; mg/day), the analyzed abundances of 15N in the FR N (15NFR; μg/mg N), the feed N (15Nfeed; μg/mg N), SAM N (15NSAM, μg/mg N), and the assumed natural abundance of 15N in unlabeled NSAM (3.66303 μg/mg):
Microbial CP (CPm; mg/d) was calculated as Nm (Eq. 1) multiplied by 6.25:
In the next step, EMPS [g/kg degraded organic matter (OM)] was determined using the CPm (mg/d), the degraded organic matter (OMdeg; mg/d), and the amount of OM originating from SAM in the FR (OMSAM; mg/d):
Where OMSAM was calculated as described by Boguhn et al. (2006) with the amount of NSAM (Eq. 3), the analyzed N concentration in SAM (N%SAM; %), the concentration of ash in SAM [12%, Boguhn et al. (2006)], and the proportion of DM in the isolated SAM fraction [0.93, Boguhn et al. (2006)]:
Statistical Analyses and 16S Bioinformatics
Statistical analysis of gas data, nutrient degradation, NH3-N, VFA, and EMPS was done with a one-way ANOVA in SAS 9.4 using the mixed procedure. The treatment was the fixed effect, and the run, circulation thermostat, fermenter, and day were used as random effects. When treatment differences were identified, multiple t-tests (Fisher’s LSD test) were used to distinguish between treatments. The residuals were checked graphically for the normal distribution and homogeneity of variance.
Fastq files were demultiplexed with Sabre1 and processed by Qiime2 (v.2021.2/8; Bolyen et al., 2019). Primers were removed by the q2-cutadapt plugin (Martin, 2011). Reads were quality filtered, error corrected, dereplicated, and merged by the q2-dada2 plugin (Callahan et al., 2016). Taxonomy assignment of produced amplicon sequence variants (ASVs) was performed using VSEARCH-based consensus (Rognes et al., 2016) and pre-fitted sklearn-based classifiers (Pedregosa et al., 2011) against the Silva SSU-rRNA database (v.138.1, 16S 99%; Quast et al., 2013), formatted by RESCRIPt (Robeson et al., 2021). Reads from organelles and unassigned sequences were removed from the analysis. A phylogenetic tree was constructed with the q2-phylogeny plugin, implementing MAFFT 7.3 (Katoh and Standley, 2013) and FastTree 2.1 (Price et al., 2010). For diversity assessment, ASV tables were rarefied to 3580 (archaea) and 5856 (bacteria) sampling depths. Alpha diversity was estimated by Faith’s phylogenetic diversity (Faith, 1992) and Shannon’s entropy (Shannon, 1948) indices. For beta diversity, Bray–Curtis (Bray and Curtis, 1957) and Jaccard (Jaccard, 1912) distances were employed. Ordination of the beta-diversity distances was implemented with a principal-coordinate analysis (PCoA; Halko et al., 2011). Differences in diversity metrics were tested with the Kruskal–Wallis H-test (Kruskal and Wallis, 1952; alpha diversity and relative abundances), followed by Dunn’s pairwise test (Dunn, 1964) and a PERMANOVA test (Anderson, 2001; beta diversity) with 999 permutations. The obtained P-values were corrected for multiple comparisons using the Benjamini–Hochberg procedure (Benjamini and Hochberg, 1995).
For correlation analysis, ASV tables were filtered to remove genera with less than 1% relative abundances. Absolute abundances of the remaining genera were correlated using “Sparse Cooccurrence Network Investigation for Compositional Data” or the SCNIC tool (Shaffer et al., 2020). Genera were correlated using the SparCC method (Friedman and Alm, 2012) and with numerical metadata using the Spearman correlation (Dodge, 2010). All P-values were considered to be significant if ≤0.05 (except for q2-ALDEx2 output) and correlation coefficients (r) if their absolute values were ≥0.3.
The q2-ALDEx2 differential abundance plugin was used to test the effect of seaweed supplementation on ASVs raw counts (relative abundance ≥1%) with a significance threshold of Wilcoxon test 0.1 (Fernandes et al., 2013, 2014). Sequences are available at the European Nucleotide Archive (ENA) under accession number PRJEB50942 and bioinformatics2.
Results
Gas Production, Nutrient Degradation, and Microbial Protein Synthesis
The total gas production was decreased by both seaweeds compared to TMR alone, with a greater effect for the 5% seaweed inclusion level (P < 0.001, Table 2). No significant difference between AN and FV was detected for total gas production. For the CH4 concentration of total gas and the methane production per g of degraded OM, only the 5% inclusion level of the two seaweeds showed a difference compared to TMR alone (P < 0.001 and 0.017, respectively). The reduction was bigger for AN5 (8.9 and 16.9% reduction for CH4 concentration and CH4/g degraded OM) compared to FV5 (3.6 and 11.2% reduction for CH4 concentration and CH4/g degraded OM).
Table 2. Total gas production and methane concentration of the produced gas (days 7–13) and methane production per g of degraded organic matter (OM; days 7–12) of the five treatments.
Seaweed supplementation also decreased the degradation of all analyzed nutrients in the feed bags (P < 0.001, Table 3). Only for ADFom, AN2.5 showed no significant difference to TMR alone. In all other cases, nutrient degradation was significantly lower with seaweed supplementation. The 5% inclusion level caused at least a numerically greater reduction than the 2.5% inclusion. The highest decrease compared to TMR alone was observed for CP, with a minimum reduction of 8.9% percentage points. The two seaweed products only differed at the 5% inclusion level for DM and OM degradation, where FV showed a greater decrease than AN. The EMPS was decreased by both seaweeds, with a greater effect of the 5% inclusion level (P < 0.001, Table 3). No significant difference was detected between the inclusion of AN and FV.
Table 3. Nutrient degradation in the feed bags of the five treatments (days 7–12) and estimated microbial protein synthesis (days 7–13).
Fermentation Characteristics
The temperature (P = 0.576; 38.6°C) and redox potential measured in FL (−241 mV; P = 0.062) did not differ among the treatments. The pH in the fermenter was different between the TMR alone (pH = 6.82) and the seaweed supplemented treatments (all pH = 6.85; P = 0.004). At both supplementation levels, seaweed decreased analyzed NH3-N and VFA production compared to TMR alone, except for propionate production in FV2.5 (Table 4; P ≤ 0.001). For isovalerate and valerate, the 5% inclusion level of both seaweeds showed a less pronounced reduction than the 2.5% inclusion level. For butyrate, this was only the case for seaweed AN and acetate only for seaweed FV. For NH3-N and the other VFA, the decrease was equal to or greater for the 5% inclusion level than the 2.5%.
Table 4. NH3-N and VFA production analyzed in the effluent of the five treatments (days 7–13).
Alpha and Beta Diversity
After sequencing and demultiplexing, a total of 4,843,225 and 5,851,973 read pairs were obtained for archaeal and bacterial datasets, respectively. Quality filtering, denoising, merging of paired reads, and chimera removal resulted in 2,648,201 archaeal and 4,237,733 bacterial ASVs. Archaeal Faith’s PD and Shannon indices were approximately equal in RF and RSP at the start of the experiment (Figures 3A,C and Supplementary Table 2) and failed to reject the null hypothesis when tested against each other (Faith’s PD P = 0.881, Shannon P = 0.531). Archaeal phylogenetic diversity was greater in the Rusitec samples compared to the inoculum; meanwhile, Shannon entropy fluctuated in the system among sampling days. Compared to RF, Faith’s PD was greater in E and FL at day 7 (P = 0.001 and 0.044) and in E at day 13 (P = 0.026). E phylogenetic diversity at day 13 decreased compared to day 7 (P = 0.048). Shannon entropy at day 7 was greater in E than RF (P = 0.006). The entropy of E and FL sample types was reduced from days 7 to 13 (P = 0.001). At day 13, both metrics demonstrated no differences between FR and SAM samples (all P > 0.05) but were greater in FR samples from bags with a 48-h incubation period compared to 24 h (Faith’s PD P = 0.028, Shannon P = 0.001). In the bacterial dataset, the alpha diversity of samples at day 0 (RF and RSP) was also at the same level (all P > 0.05) (Figures 3B,D and Supplementary Table 2). Unlike archaea, bacterial communities demonstrated no changes in E and FL samples compared by sampling days and the initial RF at day 0 (all P > 0.05). Meanwhile, the Shannon index decreased in E and FL at day 7 (P = 0.029 and 0.017) and in E at day 13 (P = 0.019). Both tested alpha diversity metrics of FR and SAM revealed no differences between these sample types (all P > 0.05) but were lower for both of them than in all samples from days 0 to 7 and E and FL at day 13 (all P < 0.05). When testing FR samples based on incubation time, no differences were observed for phylogenetic diversity, but Shannon entropy was lower in 48-h samples (P = 0.026).
Figure 3. Alpha diversity indices by sample type and sampling day. Faith’s phylogenetic diversity (A,B) and Shannon entropy index (C,D) are plotted for archaea (A,C) and bacteria (B,D). Boxplots visualize the alpha diversity metric distribution across all samples inside a given group. Gray dashed lines subdivide boxes into groups, where pairwise comparisons were performed. Significant after Benjamini–Hochberg correction, P-values of Dunn’s pairwise test are denoted by black asterisks. For FR samples, additional (dots) were plotted, indicating samples taken from bags with 24 and 48 h incubation time. Gray asterisks are plotted above FR samples if the Kruskal–Wallis test indicated a significant difference between samples with different incubation times. Significant P-values denoted as: *P ≤ 0.05, **P ≤ 0.01, and ***P ≤ 0.001.
In archaea, seaweed additives affected the Shannon entropy of E at days 7 and 13 and FL at day 13 (Supplementary Figure 1). At both days 7 and 13, diets AN5, FV2.5, and FV5 resulted in lower entropy than the TMR alone (all P < 0.05). In FL, the Shannon index of FV2.5 and FV5 samples was also lower than that of TMR (all P < 0.05). Among bacteria, no significant changes between treatments were observed according to both alpha diversity metrics (all P > 0.05).
Concerning beta diversity, PCoA analysis (Figure 4) was performed on the ASVs, using Jaccard and Bray–Curtis distances. Archaeal and bacterial samples clustered into three big groups (RF and RSP from day 0 as one group, E and FL samples from days 7 to 13 as another one, and FR and SAM from day 13 as the third), with clearer clusters in the archaeal dataset. PERMANOVA analysis revealed that microbial composition was affected by the sample type (all P < 0.05), sampling day (all P < 0.05), and between FR samples incubated for 24 and 48 h (all P < 0.01), except for archaeal Jaccard distances.
Figure 4. Principal coordinates analysis (PCoA) plots of beta diversity metrics. Jaccard (A,C) and Bray–Curtis (B,D) distances are plotted for archaea (A,B) and bacteria (C,D). Points are sized by the sampling day (days 0–13), shaped by the sample type, and colored by treatments [as stated in the legend on subplot (D)]. Label “n/a” in the legend means no treatment was applied.
Seaweed inclusion affected archaeal beta diversity according to both Jaccard and Bray–Curtis distances of E samples at days 7 and 13 and FL at day 13 (all P < 0.05). Additionally, Bray–Curtis distances were also affected in FR and SAM samples (all P < 0.05). Pairwise tests with correction for multiple comparisons found differences in Bray–Curtis distances of FR samples between the FV5 diet and TMR alone and AN2.5 (all P = 0.025). No effect of seaweed on bacterial beta diversity was detected (all P > 0.05).
Archaeal and Bacterial Taxonomy Composition
Most of the archaeal sequences (≈ 95% of the reads) were assigned to Methanomicrobium, Methanobrevibacter, unclassified to genus level members of the Methanomethylophilaceae family, Candidatus Methanomethylophilus, and Methanimicrococcus (Figure 5A). Methanobrevibacter was the most abundant genus in RSP and RF samples, with a gradual shift toward Methanomicrobium dominance in E and FL. Considering all samples, regardless of the treatment by sample type and day (Supplementary Table 3), relative abundances of Methanomicrobium were barely detectable in the samples that served as inoculum RF and RSP but increased up to the dominance in samples from the Rusitec (all P < 0.05). On the contrary, Methanobrevibacter ratios mostly decreased in the Rusitec samples (all P < 0.05, except for E at day 7 compared to RF: P = 0.065). Moreover, Methanomicrobium relative abundances in E and FL were also greater at day 13 than at day 7 (P = 0.01 and 0.001); meanwhile, by the same comparison, the presence of Methanobrevibacter at day 13 was lower for E (P = 0.03). Except for SAM, the proportions of unclassified Methanomethylophilaceae in RF were lower than in Rusitec sample types (all P < 0.05). Methanomicrobium in FL at day 13 (P = 0.023), unidentified Methanomethylophilaceae in E at day 7 (P = 0.029), Methanimicrococcus in E at day 7 (P = 0.023), and FL at day 13 (P = 0.008) had higher relative abundances in FV5 samples than TMR alone (Supplementary Table 4). At day 13, Methanimicrococcus was also more abundant in FV2.5 samples from SAM (P = 0.043).
Figure 5. Stacked taxonomy bar plots. Relative abundances of (A) archaeal and (B) bacterial genera are plotted for each sample type and sampling day by treatments (only genera with relative abundance ≥3.5% among at least one column were included). If certain features were not annotated to genus level, the last available taxonomy unit was shown. Minor Silva annotation identifiers, such as clades, taxonomic groups, clones, and others not included in commonly used nomenclature, are indicated in the round brackets.
On average, among bacteria, Lactobacillus showed the highest abundance across all samples, followed by Pseudoscardovia, unclassified Bacteroidales, marked by Silva database as F082 (further F082), Limosilactobacillus, and Prevotella (Figure 5B). The high abundance of Lactobacillus was mostly observed in FR and SAM samples. In E and FL samples, F082 dominated the community. Prevotella accounted for the highest abundance in the inoculum (day 0) but was less abundant in E, FL, FR, and SAM (all P < 0.05) compared to RF, except for E at day 7. Compared with sample types and sampling days (Supplementary Table 3), relative abundances of Lactobacillus, Pseudoscardovia, and F082 were about the same in E and FL at both days 7 and 13 and increased their proportions compared to cow samples (all P < 0.05, except for Pseudoscardovia at day 13) after incubation in the Rusitec. Notably, Lactobacillus and Pseudoscardovia were poorly represented in RF and RSP. Treatment AN5 decreased relative abundances of F082 in E samples on day 7 (P = 0.025) and FR on day 13 (P = 0.005) compared to TMR alone. Meanwhile, Ruminobacter proportions in E (P = 0.044) and SAM samples (P = 0.05) increased at day 13. FV5 also increased relative abundances of Ruminobacter in FL samples at day 13 (P = 0.022) (Supplementary Table 4).
Effect of Seaweed Additives on the Microbial Amplicon Sequence Variant Abundances
In E and FL sample types, 10 archaeal and 10 bacterial ASVs were differentially abundant in seaweed-supplemented diets compared to TMR alone (Figure 6). In E, archaeal ASV with the first four characters in corresponding id digest 7f9f (further indicated in round brackets), assigned to Candidatus Methanomethylophilus, was more abundant in treatments AN5 (P = 0.054) and FV5 (P = 0.049). Methanobrevibacter wolinii ASVs (80e7, b43f) were less abundant in FV5-supplemented treatments (P = 0.013, 0.01); meanwhile, Methanomicrococcus (48e8) abundances increased (P = 0.002). Unclassified Methanomethylophilaceae ASVs (d2f4, 8b4b, 9af5, 898c, and 37c0) were mostly associated with seaweed-supplemented samples, with one ASV (c795) more abundant in TMR alone and decreased when the FV2.5 (Wilcoxon P = 0.008) was added. In FL, Methanomicrococcus (48e8) (P = 0.005) and Methanomethylophilaceae (9af5, d2f4) (P = 0.067, 0.069) were more abundant in FV5 diet. Methanomethylophilaceae (d2f4) also increased when FV2.5 seaweed was added (P = 0.052).
Figure 6. Relative abundances of differentially abundant features according to ALDEx2. Relative abundances are normalized by rows. Archaeal (A,B) and bacterial (C,D) features from effluent (A,C) and fermenter liquid (B,D) sample types. All seaweed additives were compared with TMR. ALDEx2 significant abundances are marked by asterisks. All supplements were compared to TMR alone and separated by black dashed lines.
Among bacteria of E and FL sample types, ASVs (7718, ea75), assigned to the Bacteroidales F082 group, Pseudoscardovia (9a4e), Atopobium (e1d3), Prevotellaceae (8969), and Succiniclasticum (1814) were mostly associated with TMR alone. However, Prevotellaceae (8969) (P = 0.083) and Pseudoscardovia (9a4e) (P = 0.085) were more abundant in FL samples of AN5. Prevotella (30d0) was more abundant in AN2.5, AN5, and FV5 samples (P = 0.002, 0.001, and 0.011 in FL and 0.002, 0.003, and 0.005 in E). Rikenellaceae RC9 gut group (81ff, 01fe) also increased in E and FL sample types when AN seaweed was added (P in the range of 0.002–0.003). Finally, in E unclassified Spirochaetaceae (8023) was more abundant in AN5 (P = 0.004) and FV5 (P = 0.04) supplemented diets.
Correlations
A correlation analysis was performed to investigate potential relationships between microbiota genera, total gas, CH4 concentration, and VFA production (Supplementary Tables 5–7). At the genus level, correlations of microbial abundances with total gas production, CH4 concentration, and VFA were performed for E, FL, and FR sample types (Supplementary Figure 2), with only significant correlations plotted (absolute value of Spearman r ≥ 0.3 and P ≤ 0.05). The largest number of genera that negatively correlated with the above-mentioned traits was detected in E samples; meanwhile, in FL and FR, positive correlations prevailed.
Discussion
Seaweed Effects
A significant reduction of CH4 concentration in the total gas was found to be caused by the supplementation of both seaweed species at the 5% inclusion level. However, the effect was not as high as previously reported for A. taxiformis (Kinley et al., 2016; Machado et al., 2018; Roque et al., 2019). Additionally, the lower CH4 concentration in total gas was accompanied by a reduced rate of overall fermentation. The microbiota produced not only less methane but also fewer metabolites relevant to animals. This was indicated by the decreased total gas production, VFA, NH3-N, and nutrient degradation. Reduced fermentation and nutrient degradation would generally mean a diminished amount of nutrients and energy available for the host animal. However, in the case of CP, a reduced degradation could benefit the animal because the ruminal undegradable protein could be used as a bypass protein like it was already shown for other seaweed species (Tayyab et al., 2016). The reduced CP degradation may have resulted from the formation of complexes between the phlorotannins contained in both seaweeds and other proteins (Belanche et al., 2016) and resulted in a lower EMPS, but more research is needed to confirm this.
The CH4 reduction effect of A. taxiformis is attributed to the presence of bromoform (Machado et al., 2016), which is known to be a carcinogen (National Toxicology Program, 1989). The suggested modes of action in the brown seaweeds AN and FV used in the present study were likely related to phlorotannins as the bromoform content [8.0 (AN) and <0.8 μg/kg DM (FV)] was negligibly lower compared to the one analyzed by Machado et al. (2016) for A. taxiformis (1,723,000 μg/kg DM). The TPC, an indicator of phlorotannin content, was very similar in both seaweeds (7.9 vs. 7.4 phloroglucinol equivalent g/100 g sample). The greater CH4 reduction with the supplementation of AN compared to FV could be explained by either the slightly greater TPC or the greater bromoform content in this sample, even if it is still very low. However, both species showed the same CH4 concentration of total gas at the low seaweed inclusion level. Both bromoform and terrestrial tannins influence methanogens (Cieslak et al., 2014; Machado et al., 2018). Additionally, it has been shown that the composition of archaeal methanogens is associated with CH4 production (Tapio et al., 2017). Consistent with this finding, our study showed that the decline in methane concentration of total gas in seaweed supplemented samples was accompanied by the separation of archaeal and bacterial ASVs in a dose-dependent manner, indirectly supporting the importance of methanogen composition. Interestingly, AN5 seaweed supplementation resulted in more differentiation of microbial ASVs abundances and a stronger reduction of CH4 concentration when compared to TMR than FV5.
The treatment showing the highest CH4 reduction in this study, AN5, also contained the highest iodine concentration (Supplementary Table 1). The iodine content of AN was more than 10 times greater than that of FV. Therefore, the application of practical feeding should be done with caution. Iodine toxicity in dairy cows is reported by the National Research Council (U.S.) (2001) for a concentration of 5 mg iodine/kg dietary DM. The European Food Safety Authority recommends a maximum iodine content of 2 mg/kg of complete feed (EFSA, 2013). In the current experiment, the iodine concentration was 35 mg/kg DM for treatment AN2.5 and 2.75 mg/kg DM for treatment FV2.5. Additionally, it was shown that the iodine concentration in milk follows a dose-response relationship with iodine intake (Antaya et al., 2015; Newton et al., 2021). Humans, especially children, are even more sensitive to iodine poisoning than ruminants (Zimmermann et al., 2005). Therefore, the seaweed inclusion levels used in the present study do not apply to practical ruminant feeding, except for FV2.5. Seaweed supplementation at an exceptional iodine inclusion level would therefore not significantly reduce methane compared to TMR.
Methanogenesis
In our study, Prevotella, unclassified members of Prevotellaceae and Clostridia negatively correlated with CH4 concentration (Supplementary Table 6) and archaeal methanogens (Supplementary Table 7). Prevotella spp. and other members of the Prevotellaceae family are known for their ability to consume H2 and produce propionate (Strobel, 1992; Mitsumori et al., 2012; Denman et al., 2015), which in turn is negatively associated with CH4 formation and able to act as an alternative hydrogen sink (Ungerfeld, 2015). Inverse associations of Prevotella spp. with CH4 concentration were previously reported (Mitsumori et al., 2012; Aguilar-Marin et al., 2020). Positive correlations with CH4 concentration were observed for Fibrobacter, Lactobacillus, Pseudoramibacter, Olsenella, Shuttleworthia, and unclassified Oscillospiraceae (NK4A214). Fibrobacter, a fibrolytic bacteria that produce formate, can also be utilized by methanogens as a substrate for CH4 production (Rychlik and May, 2000). Positive correlations of the Fibrobacter genus with CH4 formation were previously reported (Wang et al., 2016). The addition of some Lactobacillus and Limosilactobacillus spp. (L. mucosae) in an in vitro rumen fermentation technique led to the rise of CH4 production (Soriano et al., 2014). Pseudoramibacter can utilize carbohydrates as an energy source (Deusch et al., 2017), producing VFA during fermentation, including butyrate, acetate, formate, and hydrogen (Palakawong Na Ayudthaya et al., 2018), and increasing methane production by methanogens. Olsenella contains genes encoding choline trimethylamine lyase (Kelly et al., 2019) and is involved in CH4 production (Broad and Dawson, 1976; Neill et al., 1978). Shuttleworthia is positively correlated with CH4 formation in dairy heifers (Cunha et al., 2019) and produces acetate and butyrate through glucose fermentation (Downes et al., 2002). Oscillospiraceae members were linked with acetate production (Tanca et al., 2017) and, in our study, correlated with acetate, the A:P ratio, and CH4 concentration. In general, those correlations are consistent with previous studies. It was shown that the addition of AN to an in vitro fermentation system led to a decrease in the growth of Fibrobacter succinogenes (positively correlated with CH4 concentration) and an increase in Prevotella bryantii (negatively correlated with CH4 concentration) (Wang et al., 2009). As it was reviewed by Abbott et al. (2020), it is possible that seaweeds modulate the growth of cellulolytic rumen bacteria by altering polysaccharide availability (Lee et al., 2019).
Adaptation of Microbiota Composition
Although we allowed rumen microbiota to adapt for 7 days in the Rusitec prior to sampling, archaeal community distribution demonstrated a shift in the dominant genera at the end of the adaptation period and revealed a variation between the samplings at days 7 and 13. Partially, it can be explained by lower metabolite digestibility in the Rusitec than in vivo conditions and lack of or low protozoa counts (Martínez et al., 2010; Hristov et al., 2012). Relative abundances of Methanobrevibacter, which is the most abundant methanogen genus in cow samples, continuously decreased in the Rusitec. This tendency was somehow more profound in seaweed-treated samples, although no statistically significant differences with TMR were detected. In the Rusitec samples, dominance gradually shifted from Methanobrevibacter to Methanomicrobium by day 13. Although it has been shown that the Rusitec provides a stable condition for the methanogens (Lengowski et al., 2016), additional studies with longer in vitro incubation periods may be performed to investigate Methanobrevibacter/Methanomicrobium ratio dynamics in such systems since abundances of Methanobrevibacter relative to other archaea may affect CH4 production (Cunha et al., 2019). When the total gas production and methane concentration are additionally statistically analyzed not only by treatment but also by sampling day and the interaction between treatment and day, it becomes evident that both traits increased during the experimental period (P = 0.005 and <0.001, respectively; Supplementary Table 8). This coincides with changes in the proportion of the most abundant methanogens in the system (Methanobrevibacter abundances decreased, while Methanomicrobium increased) between days 7 and 13 and indicates that adaptation of the microbes to the system was not finalized after day 7. Previous studies have shown a decrease in the protozoa population, also known as a methane producer, in Rusitec fermenters over time (Ziemer et al., 2000; Martínez et al., 2010; Morgavi et al., 2012). Protozoa were not analyzed in the current study, but this potential shift could also have led to an altered archaeal community distribution.
Conclusion
The North Atlantic seaweed species used herein can modify the microbiota in the rumen toward reduced methane production. However, methane reduction was not great, and the seaweed inclusion in the ration caused reduced in vitro fermentation overall. Therefore, the tested seaweed species are not a viable option to be used as feed for ruminants to mitigate methane production.
Data Availability Statement
The data presented in the study are deposited in the European Nucleotide Archive (ENA) repository, accession number PRJEB50942 (https://www.ebi.ac.uk/ena/browser/view/PRJEB50942).
Ethics Statement
The animal study was reviewed and approved by the Regierungspräsidium Stuttgart, Germany.
Author Contributions
SK, KW, MR, AC-S, HG, and AP conceived and designed the experiment. SK, KW, and HP performed the experiments. SK, KW, HP, and MR performed the in vitro data analysis. TY and AC-S performed the microbial data analysis. AP and HG carried out chemical data analysis. SK and TY drafted the manuscript. All authors revised and approved the manuscript revisions.
Funding
This study was funded by the EIT Food, the innovation community on the food of the European Institute of Innovation and Technology (EIT), a body of the EU, under Horizon 2020, the EU Framework Programme for Research and Innovation via the activity SeaCH4NGE [seaweed supplementation to mitigate methane (CH4) emissions by cattle; No. 19084].
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We gratefully acknowledge the support of Christiane Horn, Helga Ott, Sibylle Rupp, Helga Terry, Brynja Einarsdóttir, Branka Borojevic, and Natasa Desnica in conducting the chemical analyses. We also acknowledge support by the High Performance and Cloud Computing Group at the Zentrum für Datenverarbeitung of the University of Tübingen, the state of Baden-Württemberg through bwHPC, and the German Research Foundation (DFG) through grant no. INST 37/935-1 FUGG.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.889618/full#supplementary-material
Supplementary Figure 1 | Alpha diversity indices by treatment, sample-type and day of sampling. Faith’s phylogenetic diversity and Shannon entropy index are plotted for archaea and bacteria. Each subplots represents one of alpha diversity metrics at a certain day of sampling. Boxplots visualize alpha diversity metrics distribution across all treatments inside of a certain group. Kruskal–Wallis H-test was performed to test the differences between treatments alpha diversity and corresponding P-values are plotted on the top of subplots. Significant P-values denoted as ‘*’ if P ≤ 0.05.
Supplementary Figure 2 | SCNIC based correlations of microbial relative abundances among effluent, fermenter liquid and feed residue sample types with total gas production, methane as a percentage of total gas [CH4(%)], ammonium nitrogen in effluent (NH3-N), VFA production and acetate to propionate ratio (A:P ratio). Only correlation coeffitients (r) with absolute values ≥0.3 and adjusted P-values ≤ 0.05 are plotted. Positive correlations are colored in read and negative in blue.
Abbreviations
ADFom, acid detergent fiber on ash-free basis; AN, Ascophyllum nodosum; aNDFom, Neutral detergent fiber on ash free-basis; A:P, acetate to propionate ratio; ASV, amplicon sequence variant; CH4, methane; CO2, carbon dioxide; CP, crude protein, DM, dry matter; E, effluent; EMPS, estimated microbial protein synthesis; FL, fermenter liquid; FR, feed residues; FV, Fucus vesiculosus; H2, hydrogen; LAM, liquid-associated microbes; SAM, solid-associated microbes; OTU, operational taxonomic units; PCoA, principle-coordinate analysis; PCR, polymerase chain reaction; r, correlation coefficient; RF, rumen fluid; RSP, rumen solid phase; TMR, total mixed ration; TPC, total phenolic content; VFA, volatile fatty acids.
Footnotes
References
Abbott, D. W., Aasen, I. M., Beauchemin, K. A., Grondahl, F., Gruninger, R., Hayes, M., et al. (2020). Seaweed and seaweed bioactives for mitigation of enteric methane: challenges and opportunities. Animals 10:2432. doi: 10.3390/ani10122432
Aguilar-Marin, S. B., Betancur-Murillo, C. L., Isaza, G. A., Mesa, H., and Jovel, J. (2020). Lower methane emissions were associated with higher abundance of ruminal Prevotella in a cohort of colombian buffalos. BMC Microbiol. 20:364. doi: 10.1186/s12866-020-02037-6
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral. Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Antaya, N. T., Soder, K. J., Kraft, J., Whitehouse, N. L., Guindon, N. E., Erickson, P. S., et al. (2015). Incremental amounts of Ascophyllum nodosum meal do not improve animal performance but do increase milk iodine output in early lactation dairy cows fed high-forage diets. J. Dairy. Sci. 98, 1991–2004. doi: 10.3168/jds.2014-8851
Belanche, A., Jones, E., Parveen, I., and Newbold, C. J. (2016). A metagenomics approach to evaluate the impact of dietary supplementation with Ascophyllum nodosum or Laminaria digitata on rumen function in rusitec fermenters. Front. Microbiol. 7:299. doi: 10.3389/fmicb.2016.00299
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Boguhn, J., Kluth, H., and Rodehutscord, M. (2006). Effect of total mixed ration composition on fermentation and efficiency of ruminal microbial crude protein synthesis in vitro. J. Dairy Sci. 89, 1580–1591. doi: 10.3168/jds.S0022-0302(06)72225-1
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of southern wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268
Britton, D., Schmid, M., Revill, A. T., Virtue, P., Nichols, P. D., Hurd, C. L., et al. (2021). Seasonal and site-specific variation in the nutritional quality of temperate seaweed assemblages: implications for grazing invertebrates and the commercial exploitation of seaweeds. J. Appl. Phycol. 33, 603–616. doi: 10.1007/s10811-020-02302-1
Broad, T. E., and Dawson, R. M. (1976). Role of choline in the nutrition of the rumen protozoon entodinium caudatum. J. Gen. Microbiol. 92, 391–397. doi: 10.1099/00221287-92-2-391
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods. 13, 581–583. doi: 10.1038/nmeth.3869
Cieslak, A., Zmora, P., Pers-Kamczyc, E., Stochmal, A., Sadowinska, A., Salem, A. Z., et al. (2014). Effects of two sources of tannins (Quercus L. and Vaccinium vitis idaea L.) on rumen microbial fermentation: an in vitro study. Ital. J. Anim. Sci. 13:3133. doi: 10.4081/ijas.2014.3133
Cunha, C. S., Marcondes, M. I., Veloso, C. M., Mantovani, H. C., Pereira, L. G. R., Tomich, T. R., et al. (2019). Compositional and structural dynamics of the ruminal microbiota in dairy heifers and its relationship to methane production. J. Sci. Food Agric. 99, 210–218. doi: 10.1002/jsfa.9162
Czerkawski, J. W., and Breckenridge, G. (1977). Design and development of a long-term rumen simulation technique (rusitec). Br. J. Nutr. 38, 371–384. doi: 10.1079/bjn19770102
de Menezes, A. B., Lewis, E., O’Donovan, M., O’Neill, B. F., Clipson, N., and Doyle, E. M. (2011). Microbiome analysis of dairy cows fed pasture or total mixed ration diets. FEMS Microbiol. Ecol. 78, 256–265. doi: 10.1111/j.1574-6941.2011.01151.x
Denman, S. E., Martinez Fernandez, G., Shinkai, T., Mitsumori, M., and McSweeney, C. S. (2015). Metagenomic analysis of the rumen microbial community following inhibition of methane formation by a halogenated methane analog. Front. Microbiol. 6:1087. doi: 10.3389/fmicb.2015.01087
Deusch, S., Camarinha-Silva, A., Conrad, J., Beifuss, U., Rodehutscord, M., and Seifert, J. (2017). A structural and functional elucidation of the rumen microbiome influenced by various diets and microenvironments. Front. Microbiol. 8:1605. doi: 10.3389/fmicb.2017.01605
Downes, J., Munson, M. A., Radford, D. R., Spratt, D. A., and Wade, W. G. (2002). Shuttleworthia satelles gen. nov., sp. nov., isolated from the human oral cavity. Int. J. Syst. Evol. 52, 1469–1475. doi: 10.1099/00207713-52-5-1469
Dunn, O. J. (1964). Multiple comparisons using rank sums. Technometrics 6, 241–252. doi: 10.2307/1266041
EFSA (2013). Scientific Opinion on the safety and efficacy of iodine compounds (E2) as feed additives for all species: calcium iodate anhydrous and potassium iodide, based on a dossier submitted by HELM AG. EFSA J. 11:3101. doi: 10.2903/j.efsa.2013.3101
EPA (1996). Method 8260B: Volatile Organic Compounds by Gas Chromatography/Mass Spectrometry (GC/MS) (Washington, DC: United States Environmental Protection Agency).
Evans, F. D., and Critchley, A. T. (2014). Seaweeds for animal production use. J. Appl. Phycol. 26, 891–899. doi: 10.1007/s10811-013-0162-9
Faith, D. P. (1992). Conservation evaluation and phylogenetic diversity. Biol. Conserv. 61, 1–10. doi: 10.1016/0006-3207(92)91201-3
Fernandes, A. D., Macklaim, J. M., Linn, T. G., Reid, G., and Gloor, G. B. (2013). ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq. PLoS One 8:e67019. doi: 10.1371/journal.pone.0067019
Fernandes, A. D., Reid, J. N., Macklaim, J. M., McMurrough, T. A., Edgell, D. R., and Gloor, G. B. (2014). Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2:15. doi: 10.1186/2049-2618-2-15
Friedman, J., and Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687. doi: 10.1371/journal.pcbi.1002687
Halko, N., Martinsson, P.-G., Shkolnisky, Y., and Tygert, M. (2011). An algorithm for the principal component analysis of large data sets. SIAM J. Sci. Comput. 33, 2580–2594. doi: 10.1137/100804139
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., and Janssen, P. H. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5:14567. doi: 10.1038/srep14567
Hildebrand, B., Boguhn, J., and Rodehutscord, M. (2011). Effect of maize silage to grass silage ratio and feed particle size on protein synthesis and amino acid profile in different microbial fractions in a semi-continuous rumen simulation. Animal 5, 537–546. doi: 10.1017/S1751731110002156
Hristov, A. N., Lee, C., Hristova, R., Huhtanen, P., and Firkins, J. L. (2012). A meta-analysis of variability in continuous-culture ruminal fermentation and digestibility data. J. Dairy Sci. 95, 5299–5307. doi: 10.3168/jds.2012-5533
Jaccard, P. (1912). The distribution of the flora in the alpine zone. New Phytol. 11, 37–50. doi: 10.1111/j.1469-8137.1912.tb05611.x
Kaewtapee, C., Burbach, K., Tomforde, G., Hartinger, T., Camarinha-Silva, A., Heinritz, S., et al. (2017). Effect of Bacillus subtilis and Bacillus licheniformis supplementation in diets with low- and high-protein content on ileal crude protein and amino acid digestibility and intestinal microbiota composition of growing pigs. J. Anim. Sci. Biotechnol. 8:37. doi: 10.1186/s40104-017-0168-2
Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kelly, W. J., Leahy, S. C., Kamke, J., Soni, P., Koike, S., Mackie, R., et al. (2019). Occurrence and expression of genes encoding methyl-compound production in rumen bacteria. Anim. Microbiome 1:15. doi: 10.1186/s42523-019-0016-0
Kinley, R. D., de Nys, R., Vucko, M. J., Machado, L., and Tomkins, N. W. (2016). The red macroalgae Asparagopsis taxiformis is a potent natural antimethanogenic that reduces methane production during in vitro fermentation with rumen fluid. Anim. Prod. Sci. 56, 282–289. doi: 10.1071/AN15576
Kruskal, W. H., and Wallis, W. A. (1952). Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 47, 583–621. doi: 10.2307/2280779
Lassey, K. R. (2007). Livestock methane emission: from the individual grazing animal through national inventories to the global methane cycle. Agric. For. Meteorol. 142, 120–132. doi: 10.1016/j.agrformet.2006.03.028
Leahy, S. C., Kelly, W. J., Altermann, E., Ronimus, R. S., Yeoman, C. J., Pacheco, D. M., et al. (2010). The genome sequence of the rumen methanogen methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions. PLoS One 5:e8926. doi: 10.1371/journal.pone.0008926
Lee, H. J., Jung, J. Y., Oh, Y. K., Lee, S.-S., Madsen, E. L., and Jeon, C. O. (2012). Comparative survey of rumen microbial communities and metabolites across one caprine and three bovine groups, using bar-coded pyrosequencing and 1H nuclear magnetic resonance spectroscopy. Appl. Environ. Microbiol. 78, 5983–5993. doi: 10.1128/AEM.00104-12
Lee, S. J., Jeong, J. S., Shin, N. H., Lee, S. K., Kim, H. S., Eom, J. S., et al. (2019). Impact of Ecklonia stolonifera extract on in vitro ruminal fermentation characteristics, methanogenesis, and microbial populations. Asian Australas. J. Anim. Sci. 92, 1864–1872. doi: 10.5713/ajas.19.0092
Lengowski, M. B., Zuber, K. H. R., Witzig, M., Möhring, J., Boguhn, J., and Rodehutscord, M. (2016). Changes in rumen microbial community composition during adaption to an in vitro system and the impact of different forages. PLoS One 11:e0150115. doi: 10.1371/journal.pone.0150115
Machado, L., Magnusson, M., Paul, N. A., Kinley, R., de Nys, R., and Tomkins, N. (2016). Identification of bioactives from the red seaweed Asparagopsis taxiformis that promote antimethanogenic activity in vitro. J. Appl. Phycol. 28, 3117–3126. doi: 10.1007/s10811-016-0830-7
Machado, L., Tomkins, N., Magnusson, M., Midgley, D. J., de Nys, R., and Rosewarne, C. P. (2018). In vitro response of rumen microbiota to the antimethanogenic red macroalga Asparagopsis taxiformis. Microb. Ecol. 75, 811–818. doi: 10.1007/s00248-017-1086-8
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200
Martínez, M. E., Ranilla, M. J., Tejido, M. L., Saro, C., and Carro, M. D. (2010). Comparison of fermentation of diets of variable composition and microbial populations in the rumen of sheep and rusitec fermenters. II. protozoa population and diversity of bacterial communities. J. Dairy Sci. 93, 3699–3712. doi: 10.3168/jds.2009-2934
McDougall, E. I. (1948). Studies on ruminant saliva. 1. the composition and output of sheep’s saliva. Biochem. J. 43, 99–109. doi: 10.1042/bj0430099
Mitsumori, M., Shinkai, T., Takenaka, A., Enishi, O., Higuchi, K., Kobayashi, Y., et al. (2012). Responses in digestion, rumen fermentation and microbial populations to inhibition of methane formation by a halogenated methane analogue. Br. J. Nutr. 108, 482–491. doi: 10.1017/S0007114511005794
Molina-Alcaide, E., Carro, M. D., Roleda, M. Y., Weisbjerg, M. R., Lind, V., and Novoa-Garrido, M. (2017). In vitro ruminal fermentation and methane production of different seaweed species. Anim. Feed Sci. Technol. 228, 1–12. doi: 10.1016/j.anifeedsci.2017.03.012
Morgavi, D. P., Martin, C., Jouany, J.-P., and Ranilla, M. J. (2012). Rumen protozoa and methanogenesis: not a simple cause-effect relationship. Br. J. Nutr. 107, 388–397. doi: 10.1017/S0007114511002935
National Research Council (U.S.) (2001). Nutrient Requirements of Dairy Cattle. Washington D.C: National Academy Press.
National Toxicology Program (1989). Toxicology and carcinogenesis studies of tribromomethane (bromoform) (CAS No. 75-25-2) in F344/N rats and B6C3F1 mice (gavage studies). Natl. Toxicol. Program. Tech. Rep. Ser. 350, 1–198.
Neill, A. R., Grime, D. W., and Dawson, R. M. (1978). Conversion of choline methyl groups through trimethylamine into methane in the rumen. Biochem. J. 170, 529–535. doi: 10.1042/bj1700529
Newton, E. E., Pétursdóttir, Á. H., Ríkharðsson, G., Beaumal, C., Desnica, N., Giannakopoulou, K., et al. (2021). Effect of dietary seaweed supplementation in cows on milk macrominerals, trace elements and heavy metal concentrations. Foods 10:1526. doi: 10.3390/foods10071526
NMKL (2007). Trace Elements - As, Cd, Hg, Pb and other Elements. Determination by ICP-MS After Pressure Digestion (NMKL 186) (Bergen: Nordic Committee on Food Analysis).
Paiva, L., Lima, E., Neto, A. I., and Baptista, J. (2018). Seasonal variability of the biochemical composition and antioxidant properties of Fucus spiralis at two Azorean islands. Mar. Drugs 16:248. doi: 10.3390/md16080248
Palakawong Na Ayudthaya, S., van de Weijer, A. H. P., van Gelder, A. H., Stams, A. J. M., Vos, W. M., and de, et al. (2018). Organic acid production from potato starch waste fermentation by rumen microbial communities from dutch and thai dairy cows. Biotechnol. Biofuels 11: 13. doi: 10.1186/s13068-018-1012-4
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011). Scikit-learn: machine learning in python. J. Mach. Learn. Res. 12, 2825–2830.
Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2 - approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ranilla, M. J., and Carro, M. D. (2003). Diet and procedures used to detach particle-associated microbes from ruminal digesta influence chemical composition of microbes and estimation of microbial growth in rusitec fermenters. J. Anim. Sci. 81, 537–544. doi: 10.2527/2003.812537x
Robeson, M. S., O’Rourke, D. R., Kaehler, B. D., Dillon, M. R., and Bokulich, N. A. (2021). RESCRIPt: reproducible sequence taxonomy reference database management. PLoS Comput. Biol. 17:e1009581. doi: 10.1371/journal.pcbi.1009581
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Roque, B. M., Brooke, C. G., Ladau, J., Polley, T., Marsh, L. J., Najafi, N., et al. (2019). Effect of the macroalgae Asparagopsis taxiformis on methane production and rumen microbiome assemblage. Anim. Microbiome 1:3. doi: 10.1186/s42523-019-0004-4
Rychlik, J. L., and May, T. (2000). The effect of a methanogen, methanobrevibacter smithii, on the growth rate, organic acid production, and specific ATP activity of three predominant ruminal cellulolytic bacteria. Curr. Microbiol. 40, 176–180. doi: 10.1007/s002849910035
Shaffer, M., Thurimella, K., and Lozupone, C. A. (2020). SCNIC: sparse correlation network investigation for compositional data. BioRxiv [Preprint]. doi: 10.1101/2020.11.13.380733
Shannon, C. E. (1948). A mathematical theory of communication. Bell Labs Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Singleton, V. L., and Rossi, J. A. (1965). Colorimetry of total phenolics with phosphomolybdic-phosphotungstic acid reagents. Am. J. Enol. Vitic. 16, 144–158.
Soriano, A. P., Mamuad, L. L., Kim, S.-H., Choi, Y. J., Jeong, C. D., Bae, G. S., et al. (2014). Effect of Lactobacillus mucosae on in vitro rumen fermentation characteristics of dried brewers grain, methane production and bacterial diversity. Asian Australas. J. Anim. Sci. 27, 1562–1570. doi: 10.5713/ajas.2014.14517
Strobel, H. J. (1992). Vitamin B12-dependent propionate production by the ruminal bacterium Prevotella ruminicola 23. Appl. Environ. Microbiol. 58, 2331–2333. doi: 10.1128/aem.58.7.2331-2333.1992
Tanca, A., Fraumene, C., Manghina, V., Palomba, A., Abbondio, M., Deligios, M., et al. (2017). Diversity and functions of the sheep faecal microbiota: a multi-omic characterization. Microb. Biotechnol. 10, 541–554. doi: 10.1111/1751-7915.12462
Tapio, I., Snelling, T. J., Strozzi, F., and Wallace, R. J. (2017). The ruminal microbiome associated with methane emissions from ruminant livestock. J. Anim. Sci. Biotechnol. 8:7. doi: 10.1186/s40104-017-0141-0
Tayyab, U., Novoa-Garrido, M., Roleda, M. Y., Lind, V., and Weisbjerg, M. R. (2016). Ruminal and intestinal protein degradability of various seaweed species measured in situ in dairy cows. Anim. Feed Sci. Technol. 213, 44–54. doi: 10.1016/j.anifeedsci.2016.01.003
Thauer, R. K. (2012). The wolfe cycle comes full circle. Proc. Natl. Acad. Sci. U.S.A. 109, 15084–15085. doi: 10.1073/pnas.1213193109
Ungerfeld, E. M. (2015). Shifts in metabolic hydrogen sinks in the methanogenesis-inhibited ruminal fermentation: a meta-analysis. Front. Microbiol. 6:37. doi: 10.3389/fmicb.2015.00037
Verband Deutscher Landwirtschaftlicher Untersuchungs- und Forschungsanstalten (1976). Handbuch Der Landwirtschaftlichen Versuchs- und Untersuchungsmethodik (VDLUFA-Methodenbuch), BD. III Die chemische Untersuchungs von Futtermitteln mit 1.-8. Ergänzuungslieferung (1983-2012). Darmstadt: VDLUFA-Verlag.
Wang, J., Liu, M., Wu, Y., Wang, L., Liu, J., Jiang, L., et al. (2016). Medicinal herbs as a potential strategy to decrease methane production by rumen microbiota: a systematic evaluation with a focus on perilla frutescens seed extract. Appl. Microbiol. Biotechnol. 100, 9757–9771. doi: 10.1007/s00253-016-7830-z
Wang, Y., Alexander, T. W., and McAllister, T. A. (2009). In vitro effects of phlorotannins from Ascophyllum nodosum (brown seaweed) on rumen bacterial populations and fermentation. J. Sci. Food Agric. 89, 2252–2260. doi: 10.1002/jsfa.3717
Wischer, G., Boguhn, J. Steingaß, H., Schollenberger, M., Hartung, K., and Rodehutscord, M. (2013). Effect of monensin on in vitro fermentation of silages and microbial protein synthesis.Arch. Anim. Nutr. 67, 219–234. doi: 10.1080/1745039X.2013.793050
Ziemer, C. J., Sharp, R., Stern, M. D., Cotta, M. A., Whitehead, T. R., and Stahl, D. A. (2000). Comparison of microbial populations in model and natural rumens using 16S ribosomal RNA-targeted probes. Environ. Microbiol. 2, 632–643. doi: 10.1046/j.1462-2920.2000.00146.x
Keywords: seaweed, macro algae, rumen, methane, Rusitec, microbiota, 16S rRNA gene
Citation: Künzel S, Yergaliyev T, Wild KJ, Philippi H, Petursdottir AH, Gunnlaugsdottir H, Reynolds CK, Humphries DJ, Camarinha-Silva A and Rodehutscord M (2022) Methane Reduction Potential of Brown Seaweeds and Their Influence on Nutrient Degradation and Microbiota Composition in a Rumen Simulation Technique. Front. Microbiol. 13:889618. doi: 10.3389/fmicb.2022.889618
Received: 04 March 2022; Accepted: 17 May 2022;
Published: 28 June 2022.
Edited by:
Diego P. Morgavi, INRAE Clermont-Auvergne-Rhône-Alpes, FranceReviewed by:
Gilbert Kamgan Nkuekam, University of Pretoria, South AfricaCátia Pacífico, Biome Diagnostics GmbH, Austria
Stephanie A. Terry, Lethbridge Research Center, Agriculture and Agri-Food Canada, Canada
Copyright © 2022 Künzel, Yergaliyev, Wild, Philippi, Petursdottir, Gunnlaugsdottir, Reynolds, Humphries, Camarinha-Silva and Rodehutscord. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Markus Rodehutscord, aW5zdDQ1MEB1bmktaG9oZW5oZWltLmRl
†These authors have contributed equally to this work and share first authorship