- 1Sydney School of Veterinary Science, Faculty of Science, The University of Sydney, Camperdown, NSW, Australia
- 2Department of Agricultural, Forest and Food Sciences, University of Turin, Turin, Italy
- 3Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
- 4Agriculture Victoria Research, Ellinbank, VIC, Australia
- 5Centre for Agricultural Innovation, School of Agriculture and Food, Faculty of Veterinary and Agricultural Sciences, The University of Melbourne, Melbourne, VIC, Australia
- 6School of Life and Environmental Sciences, Faculty of Science, The University of Sydney, Camperdown, NSW, Australia
Ozone (O3) is volatile, highly oxidative, and has theoretical potential to reduce ruminant enteric methanogenesis by interactions between archaea and bacteria, and substrate and oxygen. The effects of O3 on the rumen microbiota, fermentation parameters, and CH4 emissions were studied through in vitro fermentation using a RUSITEC apparatus with O3 dissolved in the salivary buffer. The substrate consisted of maize silage or grain concentrates, and the treatments were (1) control (no O3) and (2) O3 at 0.07 ± 0.022 mg/L in the buffer. A 4-day adaptation period followed by a 6-day experimental period was used for measuring gas production and composition, as well as fermentation characteristics, which included ruminal volatile fatty acids (VFA) and liquid- and solid-associated microbial communities. Ozone treatment decreased total gas production by 15.4%, most notably CH4 production by 20.4%, and CH4 gas concentration by 5.8%, without compromising dry matter digestibility (DMD) of either maize silage or grain concentrates. There were no significant effects of O3 treatment on VFA production or pH. Ozone treatment reduced the relative abundance of methanogens, particularly Methanomicrobium. This study demonstrates the potential use of O3 as a method to reduce ruminant enteric methanogenesis.
Introduction
Ruminant CH4 production is notoriously one of the largest contributing factors toward greenhouse gas (GHG) emissions, accounting for 50% of emissions from the agricultural sector in Australia, and 7% to total Australian GHG emissions (Hook et al., 2010; Duarte et al., 2017b). Research toward methods to reduce CH4 emissions is increasingly valuable with increasing worldwide demand for meat and dairy products; however, minimal progress has actually been achieved (Martin et al., 2010).
Studies have shown that manipulation of feed directly or through additives has been an effective strategy (Beauchemin et al., 2008; Hristov et al., 2013; Teoh et al., 2019) in reducing CH4 production. The rumen microbiome produces volatile fatty acids (VFA) from otherwise non-digestible plant polysaccharides and these VFA are the main source of energy to ruminants (Dijkstra, 1994; Henderson et al., 2015). Ruminal archaea, namely, the methanogens, are sensitive to changes in their environment, particularly oxygen availability, and produce CH4 as a metabolic by-product under hypoxic conditions (Issazadeh et al., 2013; Patra et al., 2017). Indeed, Moate et al. (2019) have recently shown that in comparison to non-fistulated cows, rumen-fistulated cows have a lower CH4 yield, presumably due to inhibition of the methanogenic populations by oxygen.
One common method of disinfecting water is ozonation (also known as O3 disinfection), using O3 gas to destroy bacteria and viruses (Remondino and Valdenassi, 2018). Ozone gas (O3) is formed upon collision of oxygen gas (O2) with free oxygen atoms. Currently, O3 is widely used in the food industry (Brodowska et al., 2018); however, it has not yet been considered as a dietary additive to mitigate CH4 production in ruminants. Methanogens should be inhibited more by O3 compared to O2, because they are extremely sensitive to oxygen and O3 is a more powerful oxidant than O2 (Wang et al., 2004). The majority of the rumen microbiome is composed of bacteria and protozoa by mass, with archaea accounting for less than 4% of the microbial community (Matthews et al., 2019). When compounds such as hemiacetal of chloral and starch and halogenated CH4 analogs such bromochloro-methane were included in the diets of ruminants, rumen methanogens and CH4 production were observed to be severely reduced without negative effects on dry matter intake, dry matter digestibility (DMD) or liveweight gain (Trei et al., 1972; Tomkins et al., 2009; Martinez-Fernandez et al., 2018).
Thus, we hypothesized that (1) if O3 was included in an in vitro ruminal fermentation process, it would decrease CH4 production and (2) would have a CH4-inhibiting effect without affecting feed digestibility and other ruminal fermentation parameters. As such, the objective of this study was to evaluate the effect of including O3 in an in vitro ruminal fermentation system on CH4 production, DMD, the rumen microbiota, and other fermentation parameters.
Materials and Methods
This study was performed in accordance with The University of Sydney Animal Ethics Committee (Approved Protocol Number 2015/835). The cow used for ruminal contents was housed at Corstorphine Dairy, The University of Sydney (Camden Farm Dairy, Cobbitty, NSW, Australia, 34° 04′S; 150°81 69′E).
Experimental Design and Treatments
A feed consisting of maize silage and concentrates was fed through a RUSITEC apparatus mimicking in vitro ruminal fermentation, to investigate the effects of O3 supplementation on fermentation parameters, the rumen microbiota and methanogen activity ultimately leading to CH4 production. The experiment consisted of a completely randomized design with two treatments: a control and an O3-buffered treatment, and three replicates of each over the course of each 10-day experiment (Klevenhusen et al., 2014; Ertl et al., 2015; Sarnataro et al., 2019). The O3 treatment contained O3 at 0.07 mg/L that was pumped into McDougall’s buffer and then synthesized into the feed, whereas the control consisted of maize silage and concentrates only.
The study was comprised of 2 × 10-day incubation runs, whereby fermentation parameters (including pH, oxidation–reduction potential, total gas production, CH4 gas production, and effluent volume production) were measured daily for 6 days following a 4-day adaptation period. These are similar procedures used by Soliva and Hess (2007), Klevenhusen et al. (2014), Ertl et al. (2015), Soliva et al. (2015), and Sarnataro et al. (2019) in Rusitec studies. Effluent samples were collected from days 5–10 for analysis. Solid- and liquid-associated microbes were sampled at days 0, 5, and 10 in each incubation run. DMD was measured on days 6–9 using retained nylon bags. The VFA production was calculated using results from daily (days 0–10) sampling of effluent contents and effluent volume production.
Preparation of Substrate, O3, and Rumen Inoculum
The basis for each feed (control and treatment) consisted of maize silage (50% w/w, 5.5 g DM-basis) and concentrate (“grain”; Table 1) (5.5 g DM-basis for each substrate), weighed into individual porous nylon bags (11 cm × 6.5 cm, pore size 50 μm)—the daily total mass being 11 g (DM basis). Substrate components were derived from Corstorphine farm, The University of Sydney. Details of climate, harvest, processing, and substrate feed analysis are detailed by Teoh et al. (2019).
The ruminal fluid donor, a non-lactating Holstein dairy cow, was fed pasture ad libitum (including oaten hay at 6 kg DM/d). The donor was ruminally fistulated and some of the rumen inoculum was collected on the first day of each experimental run, 3 h post-morning feeding. Using cheesecloth for separation, the solid and liquid portions of the ruminal contents were separated and stored in preheated thermos containers, to be immediately transported to the laboratory (Teoh et al., 2019).
The O3 was produced by a commercial ozonator (model O3 Series BO3, CSN Global Pty Ltd., Buderim, QLD 4556, Australia) and pumped into the buffer (McDougall’s buffer) using a fish tank aerator (Aqua One Precision 750—3.5 W; 19 kPa). The daily average concentration of O3 across the experiment was 0.07 mg/L, measured twice daily for each 10-day run with a dissolved O3 meter (SKU-I-2019 model; Ozone Solutions, Hull, IA, United States).
RUSITEC Fermentation
For each run, the RUSITEC apparatus was set up with six separate fermentation vessels (800 mL each) submerged in a water bath at 39°C. Three of these vessels were intended for the control, and three for the O3 treatment. On day 1, each vessel contained three smaller, porous nylon bags in a smaller vessel submerged in approximately 780 mL (until the entire vessel was filled) of the collected liquid ruminal fluid. Of the nylon bags, one contained approximately 70 g wet weight of the collected solid ruminal contents used unprocessed to inoculate the system and another two substrate bags as described above. Each fermentation vessel was equipped with an input port for O3-buffered infusion, and an outlet port for the effluent. The fermentation vessels were continuously infused with McDougall’s buffer at a dilution rate of 33 mL/h (control), while the three treatment vessels were continuously infused with the same buffer plus O3. To mimic rumination, electric motors moved the contents of the fermentation vessels up and down so as to mix solid and liquid particles and their associated microbes (Teoh et al., 2019).
The nylon bag containing ruminal solids was allowed to incubate for 24 h, before being removed and replaced with two bags: one containing maize silage and one containing concentrates. On day 1, each vessel contained three nylon bags, and on each subsequent day there were four at any given time—and each of these substrate bags was allowed to remain in the fermentation vessel for 48 h (Teoh et al., 2019). The last two bags that were placed in the vessel at day 9 were not used in DM digestibility calculations on day 10 as they had only undergone 24 h of incubation.
Sample Collection and Analysis
Total Gas and CH4 Production
Air-tight bags (Plastigas, Linde AG, Munchen, Germany) were connected to the effluent output ports of each fermentation vessel. Each day, the total volume of gas produced by each vessel was recorded using a drum-type gas meter (Terry et al., 2018). Each effluent gas bag was manually pressurized to produce a reading on the meter.
Methane (CH4) production was determined daily from day 5 until the end of each 10-day run. To measure daily gas production, 34 mL of gas was collected from each gas bag. Some gases from the previous day’s fermentation may potentially be trapped in the tube connecting the fermentation vessel to the gas sampling bag. To prevent these stale gases from contaminating the gas samples, the first 17 mL from each sample was evacuated, and the remaining 17 mL used for analysis. This retained 17 mL gas sample was transferred into a 10 mL evacuated exetainer (Labco Ltd., High Wycombe, United Kingdom). A 3 mL subsample was taken from the exetainer and CH4 percentage was determined by gas chromatography (GC; Agilent model 7890A). The GC conditions used were as described by Teoh et al. (2019). Methane production was calculated by multiplying the total gas volume by the percentage of CH4, with correction for temperature and pressure (Duarte et al., 2017b). Results were expressed as mg CH4/g DM.
Dry Matter Digestibility (DMD), pH, Oxidation–Reduction Potential
Dry matter digestibility was calculated using the contents of the nylon bags (maize silage and concentrates) that had been fermenting for 48 h. The collection of these bags took place on days 6–9, when DNA extraction was not required. Each of the bags was washed in a washing machine set to cold-delicate for 30 min, before being dried at 60°C until a constant weight was obtained (approximately 48 h). The residue weight was recorded and used in the calculation of DMD which was expressed as a percentage (%) of the original DM substrate.
The pH of each vessel was measured daily (during nylon bag exchange) using a pH meter (TPS pH-mV-Temp Meter, Model WP-80) calibrated at 39°C. Redox potential was also measured during this time using an oxidation–reduction potential meter calibrated at 39°C (MW 500 ORP meter range of ±999 mV, accuracy of ±5 mV, Milwaukee Instruments, Inc., Rocky Mount, NC, United States; Ramos et al., 2018).
Effluent Volume and Volatile Fatty Acids (VFA) Production
Effluent volume production was measured daily using 2 L glass flasks, connected to the effluent outport of each fermentation vessel and submerged in ice to halt fermentation and microbial growth. Daily effluent was measured using a measuring cylinder and expressed as mL/day. From days 5 to 10, a 2 mL sample of effluent contents was collected for determination of concentrations of VFA. This sample was centrifuged (13,500 × g for 2 min at 5°C) and a 1.5 mL subsample of the supernatant was transferred into a 2 mL centrifuge tube and acidified with 0.3 mL metaphosphoric acid (0.20; wt/v). The subsample was then frozen at −20°C until it was analyzed for VFA concentrations using GC (Agilent model 7820A). Crotonic acid was used as internal standard. The settings used were FID set up at 250°C, air flow 300 mL per min, makeup flow (N2) ran at 30 mL per min with a capillary column (DB-FFAP, 30 m × 0.32 mm ID × 1 μm). Helium was used as a carrier gas with a flow rate of 30 mL per min. The split inlet was heated to 225°C, and set at 9.526 PSI, with He constant flow 1.5 mL per min, and split ratio 50:1. The oven temperature was programmed to 150°C (held at that temperature for 1 min) and then ramped at 5°C per min to 195°C (held 3 min). Total and individual VFA were estimated daily by multiplying VFA concentration by the volume of effluent.
Sequencing and Analysis of the Archaeal and Bacterial 16S rRNA Gene
The V4 hypervariable region of archaea and bacteria was sequenced as previously described (Terry et al., 2018; Teoh et al., 2019) using an Illumina MiSeq instrument and the MiSeq Reagent Kit with 500 cycles (Illumina, Inc., San Diego, CA, United States) according to manufacturer’s instructions. DADA2 v. 1.12.1 (Callahan et al., 2016) in R v. 3.6.1 was used to process and quality-filter all sequences. Briefly, the forward and reverse reads were trimmed to 220 and 200 bp, respectively, merged with a minimum overlap of 75 bp, and chimeras removed. Taxonomy was assigned to the remaining sequences, referred to here as operational taxonomic units (OTUs) at 100% similarity, using the RDP naïve Bayesian classifier and the SILVA SSU database release 132 (Quast et al., 2012). The OTUs classified as chloroplasts, mitochondria, or eukaryotic in origin were removed prior to analyses.
The number of OTUs per sample, Shannon diversity, and inverse Simpson’s diversity indices were calculated in R using Phyloseq 1.28.0 (McMurdie and Holmes, 2013). The bacterial structure was assessed using Bray–Curtis dissimilarities which were calculated with vegan 2.5.6 and the effect of O3 treatment was determined using PERMANOVA (adonis2 function) in R with 10,000 permutations. The LAM and SAM samples were assessed separately as they are quite different from each other in terms of microbial diversity and structure (Duarte et al., 2017b).
All 16S rRNA gene sequences were submitted to the Sequence Read Archive under BioProject number PRJNA590488.
Statistical Analysis
Fermentation data were analyzed using the MIXED procedure of SAS (SAS, Inc., 2020; SAS Online Doc 9.1.4). The model included the fixed effects of O3, day, and O3 × day interaction. Therefore, the individual fermenter was used as the experimental unit for statistical analysis. The minimum values of Akaike’s information criterion were used to select the covariance structure among compound symmetry, heterogeneous compound symmetry, autoregressive, heterogeneous autoregressive, Toeplitz, unstructured, and banded for each parameter. Treatment means were compared using the least squares mean linear hypothesis test (LSMEANS/DIFF). Significance was declared at P ≤ 0.05 and a trend was discussed when 0.05 < P ≤ 0.10.
Associations between the relative abundance of the 10 most relatively abundant genera, as well as three methanogenic taxa, and fermentation parameters were assessed using Spearman’s rank correlation coefficient (ρ). P-values were adjusted for multiple comparisons the Benjamini–Hochberg method. Fermentation parameters were fit to the non-metric multidimensional scaling (NMDS) ordinations of the Bray–Curtis dissimilarities using the envfit function in vegan with 10,000 permutations. Differentially abundant OTUs between the O3 and control treatments were determined using DESeq2 v. 1.24.0 (Love et al., 2014) in R. To account for unequal sequencing depth, the 16S rRNA gene samples were randomly subsampled to 23,900 sequences per sample prior to analysis of the diversity measures and Bray–Curtis dissimilarities. DESeq2 uses a negative binomial generalized linear model and the Wald test to identify differentially abundant OTUs. For the DESeq2 analysis, samples were not randomly subsampled and only OTUs present in ≥25% of the samples were included.
Results
Effect of O3 on Fermentation Parameters
In comparison to the control, O3 treatment approximately doubled the measured oxidation–reduction potential (P ≤ 0.01) in the buffer solutions in the AM and PM, but there was no difference in the oxidation–reduction potential in the incubation fluid of the control and O3 treatments (Table 2). The addition of O3 caused a 15.4% decrease in total gas production (P ≤ 0.01), and a 20.4% decrease in CH4 production (P ≤ 0.01), without compromising DMD of either the maize silage (P = 0.81) or concentrates (P = 0.44), nor total VFA production (P = 0.49) (Table 2). The concentration of CH4 (%) was decreased by 5.8% in the O3-buffered treatments (P = 0.01). Additionally, O3 supplementation had no effect on the pH (P = 0.87). There were no effects of O3 treatment on acetate, propionate, branched-chain VFA, valerate, or caproate production (P ≥ 0.23); however, butyrate production (mmol/d) was 7.14% greater (P = 0.04) with O3 supplementation compared to the control treatment.
Table 2. Effect of ozone (O3) on oxidation–reduction potential (Redox), dry matter disappearance (DMD), gas, methane, and volatile fatty acids (VFA) production.
Effect of O3 on the Rumen Microbiota
A total of 3491 archaeal and bacterial OTUs were identified among all samples with an average sequencing depth of 53,411 ± 1454 SEM sequences per sample. The structure of the rumen solid-associated microbiota (SAM) was strongly affected by substrate (PERMANOVA: R2 = 0.63; P < 0.01; Supplementary Figure S1) and less so by sampling time (R2 = 0.048; P < 0.01; Figures 1A,B). However, the solid-associated rumen samples were not affected by O3 treatment within each substrate type (P > 0.05). The microbial community structure of the liquid-associated rumen samples (LAM) was also not affected by O3 treatment (P > 0.05) but there was a time effect (R2 = 0.25; P < 0.01; Figure 1C). Additionally, no differentially abundant OTUs were identified between the control and O3 treatments for either the LAM or SAM samples. Among the fermentation parameters, butyrate (R2 = 0.36; P < 0.01) and total VFA (R2 = 0.28; P < 0.05) were significantly associated with the NMDS ordination of the Bray–Curtis dissimilarities for the SAM samples from the grain concentrates substrate. Methane (R2 = 0.31; P < 0.05) and propionate (R2 = 0.33; P < 0.05) concentrations as well as pH (R2 = 0.27; P < 0.05) were significantly correlated with NMDS ordination for the maize silage SAM samples. There were no differences on the alpha diversity indices identified between the O3 and control for the LAM (Supplementary Table S1). Further analysis of the SAM microbiota (Table 3) demonstrated that on day 10 the supplementation with O3 decreased the number of OTUs (richness) and the Shannon diversity index (P = 0.02) compared to the control. On day 5, none of the alpha-diversity metrics were affected by O3 treatment (P ≥ 0.63).
Figure 1. Non-metric multidimensional scaling (NMDS) of the Bray–Curtis dissimilarities for archaeal and bacterial solid-associated rumen samples by ozone treatment and sampling day for the (A) grain concentrates and (B) maize silage substrates. (C) NMDS of the liquid-associated rumen samples by sampling day and ozone treatment. In (A,B), vectors having a statistically significant association (P < 0.05) with the ordinations are included. Vector length is proportional to the degree of correlation between the fermentation parameters and the ordination.
Table 3. Effect of ozone (O3) on richness (number of OTUs) and diversity indexes of the solid-associated microbe (SAM) samples for each sampling days 5 and 10.
Five unique archaeal genera along with unclassified Methanomethylophilaceae members belonging to the Methanobacteria and Thermoplasmata classes were identified in the rumen samples. All of these taxa are methanogens. The large majority of archaeal 16S rRNA gene sequences were identified as unclassified members of the Methanomethylophilaceae family (data not shown). The overall relative abundance of these methanogenic taxa was 5.15 ± 0.28% (SEM) and 0.43 ± 0.04% in the LAM and SAM samples, respectively. Ozone treatment decreased the relative abundance of the total methanogen population in the SAM microbiota on day 10 (P = 0.02; Table 4), as well as Methanomicrobium (P = 0.05) and unclassified Methanomethylophilaceae (P = 0.01). There was no effect of O3 treatment (P ≥ 0.14) on any of the methanogenic genera in the LAM microbiota (Table 5).
Table 4. Effects of ozone (O3) treatment on the percent relative abundance of methanogenic genera on days 5 and 10 in solid-associated microbes (SAM) samples.
Table 5. Effects of ozone (O3) treatment on the percent relative abundance of the methanogenic genera on days 5 and 10 in liquid-associated microbe (LAM) samples.
The effect of O3 treatment on the 10 most relatively abundant bacterial genera on days 5 and 10 in the SAM samples is displayed in Table 6 and Figure 2. Nearly all of these genera differed in relative abundance between the grain concentrates and silage diets; however, only the Schwartzia genus was significantly affected by the O3 treatment, and only on day 5 in the grain diet. None of the other relatively abundant genera were affected by O3 supplementation. There were no significant differences in the relative abundance of the 10 most relatively abundant bacterial genera in the LAM samples by O3 treatment or between the grain concentrates and maize silage diets (P > 0.05) (Figure 3). Prevotella and Bifidobacterium were the most relatively abundant genera in the LAM and SAM samples, respectively. Associations between the 10 most relatively abundant genera as well as the three most relatively abundant methanogenic taxa and fermentation parameters were determined using Spearman’s rank correlation coefficient (Figure 4). Only the relative abundance of Methanomicrobium and Ruminobacter was significantly associated with any of the parameters (P < 0.05) and both were negatively correlated with pH. Some of the positive associations of note (Spearman’s rank correlation coefficient ρ > 0.50; P < 0.10) were Prevotellaceae YAB2003 with CH4 and propionate, Megasphaera with BCVFA and caproate, and Methanobrevibacter with BCVFA.
Table 6. Effects of ozone (O3) treatment on the 10 most relatively abundant bacterial genera on days 5 and 10 in solid-associated microbe (SAM) samples.
Figure 2. Box and whisker plots of the 10 most relatively abundant bacterial genera in the solid-associated microbe (SAM) rumen samples by ozone treatment and day (D5 or D10) for the grain (G) and maize silage (MS) substrates. The box in the box plots indicates the interquartile range (IQR) (middle 50% of the data), the middle line represents the median value, and the whiskers represent 1.5 times the IQR.
Figure 3. Box and whisker plots of the 10 most relatively abundant bacterial genera in the liquid-associated microbe (LAM) rumen samples ozone treatment. The box in the box plots indicates the interquartile range (IQR) (middle 50% of the data), the middle line represents the median value, and the whiskers represent 1.5 times the IQR.
Figure 4. Heatmap of the association between the 10 most relatively abundant bacterial genera and three methanogenic taxa and fermentation parameters based on Spearman’s rank correlation coefficient. Correlations that are statistically significant (P < 0.05) are indicated with “**” and those with P-values less than 0.10 are denoted with “*”.
Discussion
It is known that manipulation of diet through changes in composition or additives is an effective strategy to reduce enteric methanogenesis in the ruminant (Beauchemin et al., 2008; Hristov et al., 2013; Teoh et al., 2019). It has been suggested that certain bacterial and archaeal species in the rumen are more susceptible than others to unfavorable conditions (Wintsche et al., 2018). In the present study, inclusion of O3 in the buffer during in vitro fermentation led to a 15.4% decrease in total gas production, and a 20.4% decrease in CH4 production, confirming our first hypothesis. As far as we can ascertain, this is the first study to demonstrate in vitro, that ruminal methanogenesis can be inhibited by O3.
The O3 treatment reduced enteric CH4 production without compromising DMD. Although treatment with O3 did not significantly alter the bacterial and archaeal microbiota in the LAM samples and had no effect on the SAM bacterial populations, it was associated with a reduction in relative abundance of methanogens in the SAM samples and an increase in the production of butyric acid. In particular, significant effects were observed on both the SAM methanogenic microbiota and bacterial diversity and richness on day 10. Feng et al. (2013) demonstrated similar findings with a decreased methanogen diversity in rice-paddy soil exposed to elevated ground-level O3. While the most relatively abundant genera differ between the present study and that of Feng et al. (2013), similar patterns of increasing richness (number of OTUs) over time were observed.
The relative abundance of the Methanomicrobium genus, unclassified Methanomethylophilaceae, and total methanogenic 16S rRNA gene sequences was significantly reduced on day 10 in the O3-treated SAM samples (P ≤ 0.05). Similarly, Zhang et al. (2016) indicated a decrease in the abundance of methanogens in both O3-sensitive and O3-tolerant rice cultivars subjected to elevated ground-level O3. In the present study, none of the methanogenic genera in either the rumen solid or liquid samples were affected by the O3 treatment on day 5. However, by day 10, the relative abundance of methanogens in the solid samples was reduced by approximately 31%, indicating that the effect of O3 on methanogens increased with time. Methanomicrobium and other hydrogenotrophs including Methanobrevibacter utilize recycled hydrogen and carbon dioxide for methanogenesis. Ozone molecules are highly reactive and when present in the rumen will preferentially utilize H2 for conversion into H2O, reducing hydrogen ion availability for methanogens to synthesize CH4 (Zabranska and Pokorna, 2018). Competition for hydrogen thus has the ability to reduce methanogen activity, as seen in Zhao et al. (2018). A similar reduction in methanogenesis by using this hydrogen-competition mechanism has been achieved with other feed supplements such as nitrate (Zhao et al., 2018). However, this hydrogen competition mechanism is unlikely to be the only mechanism whereby ozone inhibits methanogenesis. As referred to in the introduction, ozone has a general biocidal effect against many types of microorganisms, and given the sensitivity of methanogens to oxygen, it seems that ozone may also have a specific biocidal effect against methanogens, or at least a direct inhibitory effect on methanogenesis (Brodowska et al., 2018; Remondino and Valdenassi, 2018).
Methanomethylophilaceae is a newly proposed family within the Methanomassiliicoccales order (Gaci et al., 2014). The characterized taxa in this order are hydrogen-dependent methylotrophs that produce CH4 through the reduction of methanol or methylamine (Lang et al., 2015) and are relatively abundant in the bovine rumen (Henderson et al., 2015; Holman and Gzyl, 2019). Interestingly, we recently reported a decrease in Methanomethylophilaceae members in response to hardwood biochar supplementation in a similar RUSITEC experiment (Teoh et al., 2019). We did not find any strong associations between any of the three most relatively abundant methanogenic taxa and CH4 production although the PCR primers used likely do not amplify all 16S rRNA genes found in methanogens. Methanobrevibacter was positively associated with BCVFA concentration while Methanomicrobium was negatively correlated. This may be related to the fact that Methanobrevibacter spp. and Methanomicrobium spp. have opposite growth requirements for BCVFA (Tanner and Wolfe, 1988; Miller and Lin, 2002). Schwartzia was the only bacterial genus among the relatively abundant genera that was significantly affected by O3 treatment and only in the grain substrate diet at day 5. The only species in this genus is Schwartzia succinivorans which was originally isolated from the bovine rumen and utilizes succinate as a sole energy source producing propionate (van Gylswyk et al., 1997). Schwartzia has also been previously reported to be negatively associated with CH4 emissions in cattle (Cunha et al., 2017).
Ozone is a highly reactive free radical with strong positive redox potential enabling it to attack glycoproteins and glycolipids in the bacterial cell membrane resulting in rupture and destruction of the cell (Greene et al., 2012; Megahed et al., 2018). In the present research, O3 treatment almost doubled the redox potential in the buffer solution added to the “O3” incubations. Methanogenic archaea have unique cell walls composed of pseudopeptidoglycan, or pseudomurein, components that make their cell wall similar to those of Gram-positive bacteria and this may contribute to their sensitivity to antimicrobials (Varnava et al., 2017; Doyle et al., 2019) and destructive agents such as O3 that target glycosidic bonds. Friedman et al. (2016) suggested that the redox potential effect is caused by a reduction in the number of anti-reactive oxygen species proteins in the genome of the most abundant methanogens, particularly Methanomicrobium spp. Additionally, O3 is able to cause a complete loss of function by disruption of essential enzymatic activity, as the disturbance of membrane-bound enzymes, proteins, and DNA leads to cell wall lysis (Komanapalli and Lau, 1996) and ultimately a reduction in methanogenesis and total gas production.
Differences in individual and herd CH4 emissions can be attributed to intrinsic factors such as the rumen microbiota and variances in particle retention time in the rumen (Martin et al., 2010), as well as extrinsic factors such as manipulation of feed components including silage derivatives, concentrates, and additives. Methanogenesis is largely dependent on the metabolic properties and function of certain members of the rumen bacterial microbiota, which are in turn dependent on substrate type and availability (Evans et al., 2019). Although forage substrate preparation has been shown to have little effect on fermentation parameters (Duarte et al., 2017a), forage quality is known to impact the rumen microbiota, with cattle having a higher feed conversion efficiency producing approximately 20% less CH4 (Nkrumah et al., 2006; Hegarty et al., 2007). Feed used in the present study is based on a combination of maize silage and grain concentrates at 50% w/w on a 5.5 g DM basis, consisting of 11.95% CP and 29.56% NDF. Similar studies have seen a correlation between the use of relatively high-quality forage and a reduction in ruminant CH4 production (Ramos et al., 2018); though further investigation is required to demonstrate a causal relationship as other experimental conditions are highly variable (Benchaar et al., 2001).
High concentrate diets with a greater presence of readily fermentable substrates (e.g., starch) compared to high forage diets have resulted in a decrease in CH4 emissions (Popova et al., 2011). Starch in grain silages favors propionate production rather than acetate as propionate acts as a sink for hydrogen ions and reduces the availability of H2 for methanogens (Beauchemin et al., 2008). Thus, an increase in the abundance of propionate-producing bacteria should reduce methanogenesis through the diversion of hydrogen (Myer et al., 2015). Indeed, a decreased acetate: propionate (A:P) ratio has been observed in numerous feed manipulation studies that have resulted in decreased CH4 production (Ramos et al., 2018; Williams et al., 2019). It should also be noted that the A:P ratio tended to be greater (P = 0.07) in the O3 treatment; this despite a high relative abundance of Prevotella spp., some of which are propionate producers (McCabe et al., 2015; Emerson and Weimer, 2017). This also corresponds with the positive association observed between the related Prevotellaceae YAB2003 group and propionate concentration. Taken together, our findings that ozone had no effect on the A/P ratio, no effects on the relative abundance of ruminal bacteria, but reduced the total number of methanogens and reduced production of methane indicate that ozone treatment specifically inhibits methanogenesis.
While the present study has shown O3 treatment can reduce methanogenesis in bovine ruminal fluid, there are many questions to be answered before this could be implemented as a strategy to reduce CH4 emissions from the livestock industries. For example, what would be the best and cheapest way to administer O3 to ruminants? Ozonation of drinking water offers a potential solution since substantial production benefits have been realized when using clean drinking water in cattle (Schütz, 2012). It is also unknown whether O3 administration reduces enteric CH4 production in vivo. Results from this study indicate the importance and promise of in vivo testing, although certain factors including rumen transit time, individual rumen microbiota composition, and fermentation characteristics will differ between individual cattle. Additionally, does O3 administration to ruminants have adverse impacts on their health and production? It is hoped that this experiment may stimulate further research designed to answer these questions and that O3 may eventually play a role in reducing enteric CH4 emissions from livestock industries.
Conclusion
The inclusion of O3 into a fermenter vessel via buffer decreased total gas production, notably CH4 production by 20.4%. It is important to note that neither DMD nor total VFA production was affected. Ozone, upon interaction with pre-existing oxygen in the buffer, doubled the buffer redox potential. There was no effect by O3 on other fermentation parameters, with the exception of a 7% increase in butyrate production. Ozone treatment decreased both the richness (number of OTUs) and diversity (Shannon diversity) in the solid rumen samples on day 10 (P = 0.02). Also, on day 10, O3 treatment decreased the relative abundance of the total methanogen population in the SAM samples, particularly Methanomicrobium spp. The decrease in CH4 gas production may be explained by disruption of the members of the methanogenic genera; however, further research is required in vivo to ensure the safety and efficiency of O3 for commercial agriculture.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by The University of Sydney Animal Ethics Committee (Approved Protocol Number 2015/835).
Author Contributions
AC designed the studies. LZ, EC, DH, KG, PM, and AC acquired data, read and critically revised drafts for intellectual contents, and approved the final manuscript. LZ, EC, DH, PM, and AC conducted laboratory analysis. DH and KG conducted bioinformatics. AC and DH ran statistical analysis. LZ, EC, DH, PM, and AC wrote the manuscript. All authors contributed to the article and approved the submitted version.
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.
Acknowledgments
We acknowledge Dana Tribbia (UNSW) and Alyssa Martino (UWA) for helping collecting data and running the RUSITEC experiment during their summer holidays. The authors would also like to thank Denison and Science Research Experience Scholarships (Faculty of Science, The University of Sydney, Australia) for their financial support of D. Tribbia and A. Martino.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.571537/full#supplementary-material
Supplementary Figure 1 | Non-metric multidimensional scaling (NMDS) of the Bray–Curtis dissimilarities for the archaeal and bacterial solid-associated rumen samples by substrate and treatment.
References
Beauchemin, K., Kreuzer, M., O’Mara, F., and McAllister, T. A. (2008). Nutritional management for enteric methane abatement: a review. Anim. Prod. Sci. 48, 21–27. doi: 10.1071/EA07199
Benchaar, C., Pomar, C., and Chiquette, J. (2001). Evaluation of dietary strategies to reduce methane production in ruminants: a modelling approach. Can. J. Anim. Sci. 81, 563–574. doi: 10.4141/A00-119
Brodowska, A. J., Nowak, A., and Śmigielski, K. (2018). Ozone in the food industry: principles of ozone treatment, mechanisms of action, and applications: an overview. Crit. Rev. Food Sci. Nutr. 58, 2176–2201. doi: 10.1080/10408398.2017.1308313
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. doi: 10.1038/nmeth.3869
Cunha, C. S., Veloso, C. M., Marcondes, M. I., Mantovani, H. C., Tomich, T. R., Pereira, L., et al. (2017). Assessing the impact of rumen microbial communities on methane emissions and production traits in Holstein cows in a tropical climate. Syst. Appl. Microbiol. 40, 492–499. doi: 10.1016/j.syapm.2017.07.008
Dijkstra, J. (1994). Production and absorption of volatile fatty acids in the rumen. Livest. Prod. Sci. 39, 61–69. doi: 10.1016/0301-6226(94)90154-6
Doyle, N., Mbandlwa, P., Kelly, W. J., Attwood, G., Li, Y., Ross, R. P., et al. (2019). Use of lactic acid bacteria to reduce methane production in ruminants, a critical review. Front. Microbiol. 10:2207. doi: 10.3389/fmicb.2019.02207
Duarte, A. C., Holman, D. B., Alexander, T. W., Durmic, Z., Vercoe, P. E., and Chaves, A. V. (2017a). The type of forage substrate preparation included as substrate in a RUSITEC system affects the ruminal microbiota and fermentation characteristics. Front. Microbiol. 8:704. doi: 10.3389/fmicb.2017.00704
Duarte, A. C., Holman, D. B., Alexander, T. W., Kiri, K., Breves, G., and Chaves, A. V. (2017b). Incubation temperature, but not pequi oil supplementation, affects methane production, and the ruminal microbiota in a rumen simulation technique (Rusitec) system. Front. Microbiol. 8:1076. doi: 10.3389/fmicb.2017.01076
Emerson, E. L., and Weimer, P. J. (2017). Fermentation of model hemicelluloses by Prevotella strains and Butyrivibrio fibrisolvens in pure culture and in ruminal enrichment cultures. Appl. Microbiol. Biotechnol. 101, 4269–4278. doi: 10.1007/s00253-017-8150-7
Ertl, P., Knaus, W., Metzler-Zebeli, B. U., Klevenhusen, F., Khiaosa-Ard, R., and Zebeli, Q. (2015). Substitution of common concentrates with by-products modulated ruminal fermentation, nutrient degradation, and microbial community composition in vitro. J. Dairy Sci. 98, 4762–4771. doi: 10.3168/jds.2014-9063
Evans, P. N., Boyd, J. A., Leu, A. O., Woodcroft, B. J., Parks, D. H., Hugenholtz, P., et al. (2019). An evolving view of methane metabolism in the Archaea. Nat. Rev. Microbiol. 17, 219–232. doi: 10.1038/s41579-018-0136-7
Feng, Y., Lin, X., Yu, Y., Zhang, H., Chu, H., and Zhu, J. (2013). Elevated ground-level O3 negatively influences paddy methanogenic archaeal community. Sci. Rep. 3:3193. doi: 10.1038/srep03193
Friedman, N., Shriker, E., Gold, B., Durman, T., Zarecki, R., Ruppin, E., et al. (2016). Diet-induced changes of redox potential underlie compositional shifts in the rumen archaeal community. Environ. Microbiol. 19, 174–184. doi: 10.1111/1462-2920.13551
Gaci, N., Borrel, G., Tottey, W., O’Toole, P. W., and Brugère, J. F. (2014). Archaea and the human gut: new beginning of an old story. World J. Gastroenterol. 20, 16062–16078. doi: 10.3748/wjg.v20.i43.16062
Greene, A. K., Guzel-Seydim, Z. B., and Seydim, A. C. (2012). “Chemical and physical properties of ozone,” in Ozone in food processing, eds C. O’Donnell, B. K. Tiwari, P. J. Cullen, and R. G. Rice (Oxford: Blackwell Publishing Ltd.), 19–32. doi: 10.1002/9781118307472.ch3
Hegarty, R. S., Goopy, J. P., Herd, R. M., and McCorkell, B. (2007). Cattle selected for lower residual feed intake have reduced daily methane production. J. Anim. Sci. 85, 1479–1486. doi: 10.2527/jas.2006-236
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., Collaborators, G. R. C., et al. (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
Holman, D. B., and Gzyl, K. E. (2019). A meta-analysis of the bovine gastrointestinal tract microbiota. FEMS Microbiol. Ecol. 95:fiz072. doi: 10.1093/femsec/fiz072
Hook, S. E., Wright, A. D., and McBride, B. W. (2010). Methanogens: methane producers of the rumen and mitigation strategies. Archaea 2010:945785. doi: 10.1155/2010/945785
Hristov, A. N., Oh, J., Lee, C., Meinen, R., Montes, F., Ott, T., et al. (2013). “Mitigation of greenhouse gas emissions in livestock production – a review of technical options for non-CO2 emissions,” in FAO Animal Production and Health Paper No. 177, eds P. J. Gerber, B. Henderson, and H. P. S. Makkar (Rome: FAO), 1–149.
Issazadeh, K., Nejati, P., Zare, F., and Laczai, O. (2013). Diversity of methanogenic bacteria in ecological niches. Ann. Biol. Res. 4, 36–42.
Klevenhusen, F., Deckardt, K., Sizmaz, Ö, Wimmer, S., Muro-Reyes, A., Khiaosa-ard, R., et al. (2014). Effects of black seed oil and Ferula elaeochytris supplementation on ruminal fermentation as tested in vitro with the rumen simulation technique (Rusitec). Anim. Prod. Sci. 55, 736–744. doi: 10.1071/AN13332
Komanapalli, I. R., and Lau, B. H. S. (1996). Ozone-induced damage of Escherichia coli K-12. Appl. Microbiol. Biol. 46, 610–614. doi: 10.1007/s002530050869
Lang, K., Schuldes, J., Klingl, A., Poehlein, A., Daniel, R., and Brunea, A. (2015). New mode of energy metabolism in the seventh order of methanogens as revealed by comparative genome analysis of Candidatus Methanoplasma termitum. Appl. Environ. Microbiol. 81, 1338–1352. doi: 10.1128/AEM.03389-14
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Martin, C., Morgavi, D. P., and Doreau, M. (2010). Methane mitigation in ruminants: from microbe to the farm scale. Animal 4, 351–365. doi: 10.1017/S1751731109990620
Martinez-Fernandez, G., Duval, S., Kindermann, M., Schirra, H. J., Denman, S. E., and McSweeney, C. S. (2018). 3-NOP vs. halogenated compound: methane production, ruminal fermentation and microbial community response in forage fed cattle. Front. Microbiol. 9:1582. doi: 10.3389/fmicb.2018.01582
Matthews, C., Crispie, F., Lewis, E., Reid, M., O’Toole, P. W., and Cotter, P. D. (2019). The rumen microbiome: a crucial consideration when optimising milk and meat production and nitrogen utilisation efficiency. Gut Microbes 10, 115–132. doi: 10.1080/19490976.2018.1505176
McCabe, M. S., Cormican, P., Keogh, K., O’Connor, A., O’Hara, E., Palladino, R. A., et al. (2015). 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 10:e0133234. doi: 10.1371/journal.pone.0133234
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Megahed, A., Aldridge, B., and Lowe, J. (2018). The microbial killing capacity of aqueous and gaseous ozone on different surfaces contaminated with dairy cattle manure. PLoS One 13:e0196555. doi: 10.1371/journal.pone.0196555
Miller, T. L., and Lin, C. (2002). Description of Methanobrevibacter gottschalkii sp. nov., Methanobrevibacter thaueri sp. nov., Methanobrevibacter woesei sp. nov. and Methanobrevibacter wolinii sp. nov. Int. J. Syst. Evol. Microbiol. 52, 819–822. doi: 10.1099/00207713-52-3-819
Moate, P. J., Williams, S. R. O., Deighton, M. H., Hannah, M. C., Ribaux, B. E., Morris, G. L., et al. (2019). Effects of feeding wheat or corn and rumen fistulation on milk production and methane emissions of dairy cows. Anim. Prod. Sci. 59, 891–905. doi: 10.1071/AN17433
Myer, P. R., Smith, T. P., Wells, J. E., Kuehn, L. A., and Freetly, H. C. (2015). Rumen microbiome from steers differing in feed efficiency. PLoS One 10:e0129174. doi: 10.1371/journal.pone.0129174
Nkrumah, J. D., Okine, E. K., Mathison, J. W., Schmid, K., Li, C., Basarab, J. A., et al. (2006). Relationships of feedlot feed efficiency, performance, and feeding behavior with metabolic rate, methane production, and energy partitioning in beef cattle. J. Anim. Sci. 84, 145–153. doi: 10.2527/2006.841145x
Patra, A., Park, T., Kim, M., and Zhongtang, Y. (2017). Rumen methanogens and mitigation of methane emission by anti-methanogenic compounds and substances. J. Anim. Sci. Biotechnol. 8:13. doi: 10.1186/s40104-017-0145-9
Popova, M., Martin, C., Eugène, M., Mialon, M. M., Doreau, M., and Morgavi, D. P. (2011). Effect of fibre- and starch-rich finishing diets on methanogenic Archaea diversity and activity in the rumen of feedlot bulls. Anim. Feed Sci. Technol. 16, 113–121. doi: 10.1016/j.anifeedsci.2011.04.060
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). 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
Ramos, A. F. O., Terry, S. A., Holman, D. B., Breves, G., Pereira, L. G. R., Silva, A. G. M., et al. (2018). Tucumã oil shifted ruminal fermentation, reducing methane production and altering the microbiome but decreased substrate digestibility within a RUSITEC fed a mixed hay – concentrate diet. Front. Microbiol. 9:1647. doi: 10.3389/fmicb.2018.01647
Remondino, M., and Valdenassi, L. (2018). Different uses of ozone: environmental and corporate sustainability. Literature review and case study. Sustainability 10:4783. doi: 10.3390/su10124783
Sarnataro, C., Petri, R. M., Spanghero, M., Zebeli, Q., and Klevenhusen, F. (2019). A nutritional and rumen ecological evaluation of the biorefinery by-product alfalfa silage cake supplemented with Scrophularia striata extract using the rumen simulation technique. J. Sci. Food Agric. 99, 4414–4422. doi: 10.1002/jsfa.9676
Schütz, K. (2012). Effects of Providing Clean Water on the Health and Productivity of Cattle. Available at: https://envirolink.govt.nz/assets/Envirolink/1051-NLRC142-Effects-of-providing-clean-water-on-the-health-and-productivity-of-cattle.pdf (accessed on 11 June 2020)
Soliva, C., and Hess, D. (2007). ““Measuring methane emission of ruminants by in vitro and in vivo techniques,”,” in Measuring Methane Production from Ruminants, eds H. P. S. Makkar and P. E. Vercoe (Dordrecht: Springer), 15–135. doi: 10.1007/978-1-4020-6133-2_2
Soliva, C. R., Amelchanka, S. L., and Kreuzer, M. (2015). The requirements for rumen-degradable protein per unit of fermentable organic matter differ between fibrous feed sources. Front. Microbiol. 6:715. doi: 10.3389/fmicb.2015.00715
Tanner, R. S., and Wolfe, R. S. (1988). Nutritional requirements of Methanomicrobium mobile. Appl. Environ. Microbiol. 54, 625–628. doi: 10.1128/aem.54.3.625-628.1988
Teoh, R., Caro, E., Holman, D. B., Joseph, S., Meale, S. J., and Chaves, A. V. (2019). Effects of hardwood biochar on methane production, fermentation characteristics, and the rumen microbiota using rumen simulation. Front. Microbiol. 10:1534. doi: 10.3389/fmicb.2019.01534
Terry, S. A., Ramos, A. F. O., Holman, D. B., McAllister, T. A., Breves, G., and Chaves, A. V. (2018). Humic substances alter ammonia production and the microbial populations within a RUSITEC fed a mixed hay–concentrate diet. Front. Microbiol. 9:1410. doi: 10.3389/fmicb.2018.01410
Tomkins, N. W., Colegate, S. M., and Hunter, R. A. (2009). A bromochloromethane formulation reduces enteric methanogenesis in cattle fed grain-based diets. Anim. Prod. Sci. 49, 1053–1058. doi: 10.1071/EA08223
Trei, J. E., Scott, G. C., and Parish, R. C. (1972). Influence of methane inhibition on energetic efficiency of lambs. J. Anim. Sci. 34, 510–515. doi: 10.2527/jas1972.343510x
van Gylswyk, N. O., Hippe, H., and Rainey, F. A. (1997). Schwartzia succinivorans gen. nov., sp. nov., another ruminal bacterium utilizing succinate as the sole energy source. Int. J. Syst. Bacteriol. 47, 155–159. doi: 10.1099/00207713-47-1-155
Varnava, K. G., Ronimus, R. S., and Sarojini, V. (2017). A review on comparative mechanistic studies of antimicrobial peptides against archaea. Biotechnol. Bioeng. 114, 2457–2473. doi: 10.1002/bit.26387
Wang, G., Bates-Kenney, S. R., Tao, J. Q., Phelps, D. S., and Floros, J. (2004). Differences in biochemical properties and in biological function between human SP-A1 and SP-A2 variants, and the impact of ozone-induced oxidation. Biochemistry 43, 4227–4239. doi: 10.1021/bi036023i
Williams, S. R. O., Hannah, M. C., Jacobs, J. L., Wales, W. J., and Moate, P. J. (2019). Volatile fatty acids in ruminal fluid can be used to predict methane yield of dairy cows. Animals 9:1006. doi: 10.3390/ani9121006
Wintsche, B., Jehmlich, N., Popp, D., Harms, H., and Kleinsteuber, S. (2018). Metabolic adaptation of methanogens in anaerobic digesters upon trace element limitation. Front. Microbiol. 9:405. doi: 10.3389/fmicb.2018.00405
Zabranska, J., and Pokorna, D. (2018). Bioconversion of carbon dioxide to methane using hydrogen and hydrogenotrophic methanogens. Biotechnol. Adv. 36, 707–720. doi: 10.1016/j.biotechadv.2017.12.003
Zhang, J., Tang, H., Zhu, J., Lin, X., and Feng, Y. (2016). Divergent responses of methanogenic archaeal communities in two rice cultivars to elevated ground-level O3. Environ. Pollut. 213, 127–134. doi: 10.1016/j.envpol.2016.01.062
Keywords: enteric, methanogenesis, rumen, cattle, O3, CH4 (methane)
Citation: Zhao L, Caro E, Holman DB, Gzyl KE, Moate PJ and Chaves AV (2020) Ozone Decreased Enteric Methane Production by 20% in an in vitro Rumen Fermentation System. Front. Microbiol. 11:571537. doi: 10.3389/fmicb.2020.571537
Received: 11 June 2020; Accepted: 30 September 2020;
Published: 02 November 2020.
Edited by:
George Tsiamis, University of Patras, GreeceReviewed by:
Birbal Singh, Indian Veterinary Research Institute (IVRI), IndiaMin Wang, Nanjing Agricultural University, China
Zhi Yuan Ma, Chinese Academy of Sciences (CAS), China
Copyright © 2020 Zhao, Caro, Holman, Gzyl, Moate and Chaves. 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: Alex V. Chaves, alex.chaves@sydney.edu.au
†These authors have contributed equally to this work