Next Article in Journal
The Relationship between Uterine, Fecal, Bedding, and Airborne Dust Microbiota from Dairy Cows and Their Environment: A Pilot Study
Next Article in Special Issue
Methane and Carbon Dioxide Emission of Beef Heifers in Relation with Growth and Feed Efficiency
Previous Article in Journal
Effect of Dietary Supplementation of Lactobacillus Casei YYL3 and L. Plantarum YYL5 on Growth, Immune Response and Intestinal Microbiota in Channel Catfish
Previous Article in Special Issue
The Role of Chitosan as a Possible Agent for Enteric Methane Mitigation in Ruminants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Volatile Fatty Acids in Ruminal Fluid Can Be Used to Predict Methane Yield of Dairy Cows

Agriculture Victoria Research, Ellinbank, VIC 3821, Australia
*
Author to whom correspondence should be addressed.
Animals 2019, 9(12), 1006; https://doi.org/10.3390/ani9121006
Submission received: 9 November 2019 / Accepted: 13 November 2019 / Published: 20 November 2019
(This article belongs to the Special Issue Reducing Enteric Methane Emissions from Ruminants)

Abstract

:

Simple Summary

Methane emissions from cattle are difficult to measure, and some proxies used to estimate them have required information that is not always available. An example is predicting methane from the fatty-acid profile of milk, but this strategy is not suitable for non-lactating animals. We propose equations to predict the methane emitted per unit of feed eaten (methane yield) based on the volatile fatty acids within the rumen of the animals. Three of the seven equations we investigated were equally good at predicting the methane yield of dairy cows. Validation of these equations using previously published results indicated that the equations should also work for beef cattle. Being able to predict the methane yield for all classes of cattle means that a single strategy can be used, eliminating differences because of the use of different methods for different animal classes. Further work is necessary, but our strategy should be able to be adapted for use in cattle production environments. Being able to predict the methane of production animals will enable accurate estimates of the methane emissions from those animals, and assessment of strategies to reduce those emissions.

Abstract

The dry matter intake (DMI) of forage-fed cattle can be used to predict their methane emissions. However, many cattle are fed concentrate-rich diets that decrease their methane yield. A range of equations predicting methane yield exist, but most use information that is generally unavailable when animals are fed in groups or grazing. The aim of this research was to develop equations based on proportions of ruminal volatile-fatty-acids to predict methane yield of dairy cows fed forage-dominant as well as concentrate-rich diets. Data were collated from seven experiments with a total of 24 treatments, from 215 cows. Forage in the diets ranged from 440 to 1000 g/kg. Methane was measured either by open-circuit respiration chambers or a sulfur hexafluoride (SF6) technique. In all experiments, ruminal fluid was collected via the mouth approximately four hours after the start of feeding. Seven prediction equations were tested. Methane yield (MY) was equally best predicted by the following equations: MY = 4.08 × (acetate/propionate) + 7.05; MY = 3.28 × (acetate + butyrate)/propionate + 7.6; MY = 316/propionate + 4.4. These equations were validated against independent published data from both dairy and beef cattle consuming a wide range of diets. A concordance of 0.62 suggests these equations may be applicable for predicting methane yield from all cattle and not just dairy cows, with root mean-square error of prediction of 3.0 g CH4/kg dry matter intake.

Graphical Abstract

1. Introduction

Mitigation of enteric methane from ruminants requires methods to estimate the extent of mitigation. Many methods rely on complex equipment and skilled people, and while they are suitable for research purposes, they are unsuitable for routine use on farms [1,2]. An alternative has been to use proxies of methane production (g/d) such as dry matter intake (DMI) and milk fatty acids [3].
The total DMI of cattle fed high-forage diets can be used to accurately predict their methane emissions [4]. However, many animals are fed enriched diets in which concentrate constitutes more than 300 g/kg of the diet [5,6] and these diets decrease the methane yield (g/kg DMI) of those animals [7]. A range of equations to predict methane production (g/day, 34 equations) and methane yield, (7 equations) were compared by Niu et al. [8]. These equations were based on various combinations of DMI and chemical descriptors of the consumed feed (e.g., concentrations of neutral detergent fiber and ether extract) as well as animal descriptors such as bodyweight, milk yield, and milk fat concentration. Niu et al. [8] concluded that the more complex equations were better than a simple equation based on DMI alone. The more complex approach may be suitable for situations where the feed intake of individual animals can be accurately measured and where the composition of the feed consumed by each individual animal is known. However, such an approach is unlikely to apply where DMI is unknown or difficult to measure, such as animals fed in groups or grazing.
Milk fatty acids are an alternate proxy to predict methane production and methane yield [9,10]. However, measurement of long-chain fatty acids in milk by gas chromatography is relatively expensive and the approach has been found to be inaccurate at predicting methane emissions and yield if cows are fed diverse diets [11].
Production of volatile fatty acids (VFA) in ruminal fluid, including acetate, propionate, and total butyrate (n-butyrate plus iso-butyrate) have been related to methane production using stoichiometric equations [12,13]. In addition, changes in the acetate to propionate ratio because of treatment effects have been associated with changes in methane production both in vitro [14] and in vivo [15,16], while changes in (acetate + butyrate)/propionate have also been associated with changes in methane production in vivo [13,17]. While no correlations were reported between ruminal VFA proportions and methane yield, the associations suggest that VFA proportions in ruminal fluid could be used as a simple proxy to predict methane yield of dairy cows.
In previous work, researchers have noted associations between VFA proportions in ruminal fluid and methane production within their single experiment [18,19], but none have presented equations to predict methane yield from the proportions of volatile fatty acids in ruminal fluid. Therefore, the aim of our research was to use a database of methane yield and proportions of specific VFA in ruminal fluid to develop and validate equations for their accuracy in predicting methane yield of dairy cows.
We hypothesized that the methane yield of dairy cows offered a wide range of diets could be predicted from the proportions of VFA (mol/100 mol total VFA) in ruminal fluid.

2. Materials and Methods

In all collated experiments, animals were cared for according to the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes [20,21]. Animal use was approved by the Agricultural Research & Extension Animal Ethics Committee of the Department of Jobs, Precincts and Regions–Victoria (Approvals 2009-31, 16 December 2009; 2010-16, 23 August 2010; 2011-24, 06 December 2011; 2012-15, 21 August 2012; 2012-27, 13 December 2012; 2013-22, 16 January 2014; 2014-12, 23 September 2014).

2.1. Model Development

Models were developed using theoretical stoichiometry and previously published empirical associations. The abbreviations used in the models are listed in Table 1 with a brief description and their units.
A stoichiometric equation from Moss et al. [13] served as a basis for the development of equations to predict the methane yield of dairy cows. Equation (1a) is as presented by Moss et al. [13], with methane production designated by CH4P, and the substrates acetate by C2, propionate by C3, and butyrate by C4:
CH4P = 0.50C2 − 0.25C3 + 0.50C4,
Note that in Equation (1a) CH4P, C2, C3, and C4 have units of moles. However, individual VFA are often reported as moles per 100 moles of total VFA (also designated as M%), and total VFA is often reported as a concentration (mol/L). Therefore, it is necessary to do a conversion to take the units for each parameter into account. If we divide both sides of Equation (1a) by 100 moles of total VFA, then we can rewrite Equation (1a) as:
CH4Y = 0.50A − 0.25P + 0.50B,
where CH4Y is the yield of methane per 100 moles of total VFA in ruminal fluid, and A refers to acetate, P to propionate, and B to total butyrate (i.e., n-butyrate plus iso-butyrate), and these are expressed as moles per 100 moles of total VFA. If it is assumed that each kg DM consumed is fermented in the rumen to produce c moles of total VFA, then it follows that Equation (1b) can be expressed as:
MY = 16 × (0.50A − 0.25P + 0.50 B) × c/100,
where MY is defined as methane yield (g CH4/kg DMI), 16 refers to the molecular mass of methane (16 g/mole), and c has units of moles of total VFA/kg DMI. An inherent assumption in Equation (1c) is that A, P, and B at the time of sampling reflect the daily production rates of these VFA. However, the concentrations of VFA in a sample of ruminal fluid reflect the net accumulation of fatty acids because of the production and removal of these fatty acids from the rumen by the bulk flow of fluid from the rumen and by absorption of VFA across the ruminal wall into the blood. Each VFA have different molecular weights and are absorbed across the ruminal wall at different rates [22]. Therefore, to take account of these issues, we have examined additional models that are based on various combinations of A, B, and P and in which the units for MY, A, B, and P are the same as those in Equation (1c). In particular, to account for possible differences in net accumulation of A, B, and P within the rumen, we propose Equation (2).
MY = dA − eP + fB,
Equations (1c) and (2) are based purely on stoichiometry, while the following equations follow a more empirical approach. Equation (3) is based on a commonly noted association between methane production and A/P [15,16].
MY = g(A/P) + h,
Equation (4) is based on a commonly noted association between methane production and (A + B)/P [13,17].
MY = i(A + B)/P + j,
Methane production has been reported to be positively associated with the molar concentration of propionate in ruminal fluid [23]. However, it is usually suggested that the association is negative [13] as the production of methane and propionate compete for metabolic hydrogen. Accordingly, Equation (5) examines the relationship between MY and P.
MY = kP + m,
Equation (6) is similar to Equation (5), but [P] refers to the concentration (mM/L) of propionate in the ruminal fluid.
MY = n[P] + q,
Equation (7) is a simplification of Equations (3) and (4).
MY = s/P + t,
For Equations (2)–(7), the coefficients d, e, f, k, and q have units of g CH4·kg DMI−1·moles−1·100 moles total VFA, while the coefficients g, and i have units of g CH4·kg DMI−1, and the coefficient m has units of g CH4·kg DMI−1·L·mmol−1. The constants h, j, l, n, and s have the units of g CH4·kg DMI−1. All of the above coefficients and constants are determined by regression as described in the statistical analysis section.

2.2. Data

Individual animal data used to develop the models were collated from seven previously conducted experiments (Table 2).
There were 215 records of methane yield with corresponding VFA proportions. Methane production was measured either by open-circuit respiration chambers over two days [32] or a sulfur hexafluoride (SF6) technique over three to five days [33]. Methane yields of individual animals ranged from 6.9 to 32.4 g/kg DMI (Table 3). Methane yield was positively linearly related to proportions of acetate (r = 0.557) and butyrate (r = 0.642), and positively linearly related to the reciprocal of the proportion of propionate (r = 0.745). In each experiment, ruminal fluid was collected via the mouth on one occasion at the conclusion of the methane measurement period, approximately four hours after the start of the morning feeding. Samples were collected using an oro-ruminal sampling probe, similar to the one described by Geishauser [34], and a vacuum pump [35]. Concentrations of total VFA (mol/L) and of the major individual VFA (acetic, propionic, n-butyric, and iso-butyric; mol/100 mol of total VFA) in ruminal fluid were measured by gas chromatography [35].
Additional individual-cow data for external validation were taken from an experiment of similar design involving 32 Holstein cows [6] (Table 4). Further data for validation came from 67 dietary treatment means published in the scientific literature (Table 4). These data from the scientific literature included data from both dairy cows and beef animals consuming a wide range of diets. Only articles with methane production and DMI or methane yield, and individual VFA proportions where methane was measured in calorimeters were selected. Dietary treatments containing nitrate, sulphate, 3-nitrooxypropanol, or halogenated methane analogues were not used in the validation data since they can cause changes in the methane yield with little effect on ruminal VFA.

2.3. Statistical Analysis and Model Evaluation

Data from individual experiments in Table 2 and their pooled data were plotted and visually inspected for linearity or curvature between methane yield and proportions of acetate, propionate, butyrate, reciprocal of propionate and acetate to propionate ratio.
Equations (1)–(7) were modelled statistically in Genstat 18 (VSN International Ltd., Hemel Hempstead, UK) using a leave-one-experiment-out cross-validation process [52,53]. All the individual cow data from one experiment were put aside as a validation set and the remaining data were used to construct a linear mixed-effects meta-analysis, with fixed-effects given by Equations (1)–(7) above. All meta-analysis models had random effects for experiment, treatment within experiment, and animal within treatment. The estimated fixed effects from each model were then used to predict methane yield data from the omitted experiment using its VFA proportions data. In this way, each experiment was used in turn as validation data, thus providing a complete set of predicted methane yields for the seven experiments for each model. These predicted methane yields were compared to the measured methane yields using the root-mean-square-error of prediction (RMSEP) calculated to assess the performance of each model:
RMSEP   =   1 n i = 1 n ( y i ŷ i ) 2
where yi are observed validation data and ŷi are model predictions calculated independently of the validation data.
The root-mean-square-error of prediction can be interpreted as an average standard deviation of prediction. Lin’s concordance correlation coefficient (CCC) was also used as a measure of agreement between predicted and measured methane yield for each model [54]. The models were ranked according to their cross-validation RMSEP and CCC.
Each model was also fitted to the pooled training data of Table 2 using the mixed-effects meta-analysis model described above, to obtain a final calibrated equation for use in prediction beyond these data. The best ranked equations were validated using external data (Table 4). The methane yields of the external data were predicted by applying the equations, calibrated on the pooled training data to the VFA proportions reported in the external data. These predictions were compared to their matching reported methane yields using the measures of agreement, RMSEP and CCC. The square of RMSEP (MSEP) was decomposed into three components, these being the error-of-central-tendency, error-due-to-regression, and error-of-disturbance [55] also known as the squared-bias, non-unity-of-slope, and lack-of-correlation [56], respectively.

3. Results

Visual inspection of the correlations between methane yield and proportions of acetate, propionate and butyrate in the pooled data showed that methane yield was positively linearly related to the proportions of acetate and butyrate, and the reciprocal of propionate (Figure 1). Within experiments, methane yield was always positively linearly related to proportions of acetate and butyrate, and was always related to the reciprocal of propionate proportions. In some experiments, and certainly for the pooled data, the relationship between methane yield and propionate proportion did not appear to be linear, whereas the relationship between methane yield and the reciprocal of propionate proportion did appear to be linear (data not shown).
Most models had similar performance on the training data in terms of their cross-validated prediction error and concordance. The exception was Model 6, based on the concentration of propionate, which had RMSEP much greater than all other models (Table 5). The best three models were Model 3, Model 4, and Model 7. These models all had the same RMSEP of 3.2 g/kg DMI. While Model 3 had a slightly lower concordance (0.69) compared to Models 4 and 7 (0.70).
The symbolic forms of all models are shown in Table 5 alongside the coefficient estimates for these models, fitted to the pooled training data of all seven experiments. The equation derived from stoichiometry (Equation 1c, Table 5) was the second poorest model, in terms of having the second greatest RMSEP. Ranking of the concordance between measured and predicted methane yield for the models tested reflected the ranking of the RMSEP. The equations for Models 3, 4, and 7 are shown against the combined training data in Figure 2. Components of variance for deviations were 1.7 for experiment, 1.5 for treatment within experiment, and 6.9 for cow, for Model 4 (similar for Models 3 and 7). Each component was statistically significant (p < 0.05).
The RMSEP between the observed methane yield of individual cows in unrelated datasets [6] and the methane yield predicted by Equation (4) was 3.62 g/kg DMI and concordance correlation coefficient was 0.66. The components of MSEP were 3.3 for error-of-central-tendency, 0.47 for error-due-to-regression, 9.3 for error-of-disturbance. When group means from published experiments (Table 4) were used, the RMSEP was 2.99 g/kg DM and concordance was 0.62 (Figure 3). The components of MSEP were 0.62 for error-of-central-tendency, 0.20 for error-due-to-regression, and 8.1 for error-of-disturbance. Similar results were obtained for Equations (3) and (7) (data not shown).

4. Discussion

Methane yield of dairy cows can be predicted from the proportions of VFA in their ruminal fluid, confirming our hypothesis, but the concordance between observed and predicted methane yield differed between the seven models tested. Equation (3) is consistent with previous research that reported a proportional relationship between methane yield and acetate to propionate ratio [15,18,57]. However, there are also reports of no relationship [37,45,51]. The reason for this is not clear but may relate to the different dietary treatments used, and this is discussed later as a limitation of our prediction models. Equation (4) reflects previous reports of methane yield being proportional to (A + B)/P [17,19]. The relationship between methane yield and (A + B)/P was consistent both within and between the 7 experiments in our study, notwithstanding significant variation associated with experiment, treatments, and animals around the average line. Summing the components of variance estimates for experiment and treatment imply a lack-of-fit standard deviation of 1.8 g CH4/kg DMI associated with feed-type and context, with the remaining error (SD =2.6 g CH4/kg DMI) being due to animal or measurement. Equation (7) supports the contention that the production of methane and propionate are in competition [58] but no reports linking methane yield solely with the reciprocal of propionate proportion were found. Propionate proportion alone (Equation 5) was not a good linear predictor of methane yield, which is to be expected given that the data suggested a reciprocal relationship between methane yield and propionate proportion. Despite this, in a previous study on 532 Angus bulls (1 to 2 years age) fed a diet containing lucerne and oaten hay chaff, a significant positive correlation between methane yield and propionate proportion was reported [59]. However, in the work of Herd et al. [59], feed was offered once per day and ruminal fluid was sampled 24 h after the start of feeding compared to the twice per day feeding and sampling 4 h after the morning feed used in our investigations. Subsequent work by Smith et al. [60], in which ruminal fluid was sampled 3 h after the start of feeding, reported a weak negative correlation between methane yield and individual VFA proportions.
Validation of our best equations showed there was a strong concordance between the observed and predicted methane yield. The experiment of Moate et al. [6] used the same sampling and measurement processes as those experiments used to develop the equations, so strong concordance was expected. Despite the literature data being from conditions different to ours (sampling, diets, and animals), the concordance between the mean methane yield reported in the scientific literature and the methane yield as predicted by Equations (3), (4) and (7) was greater than 0.60. The components of MSEP show that the lack of agreement between the independently observed and predicted methane yields was dominated by error of disturbance, that is, random error rather than either mean or slope bias. This suggests that our best equations are robust and using them to predict methane yield from ruminal VFA proportions may be applicable for predicting methane yield from all cattle and not just dairy cows. Thus, our equations offer a low-cost approach to screen and rank large numbers of animals for methane yield.
Equations derived from stoichiometry can be used to predict the quantity of methane produced from the quantity of individual VFA in the rumen, but caution is required when applying them. It must be remembered that Equation (1a) predicts the quantity of methane from a quantity of ruminal VFA and as such has applicability at the one instant when the ruminal VFA are sampled. Even after translating Equation (1a) into terms of methane yield (Equation 1c), its performance was inferior to five of the other six models. Other research in this area has focused on using feed intake and feed composition to predict ruminal VFA production and then using stoichiometry (Equation 1a) to predict methane production from predicted VFA production [61]. Despite the sophistication of the model by Alemu et al. [61], it had poor correlations between measured and predicted VFA concentrations and hence the model of Alemu et al. [61] had poor accuracy at predicting methane production. It is important to note that individual VFA may be absorbed from the rumen at different rates [22]. A consequence of this is that proportions of individual VFA in ruminal fluid may not reflect the proportions at which they were produced, which may explain why Equation (1c) had a large RMSEP. What remains unclear is why Equation (2) performed only slightly better than Equation (1c), given that adjusting the coefficients should have taken into account any difference between proportions of individual VFA produced and their proportions in ruminal fluid.
Using the proportions of total butyrate instead of n-butyrate in our prediction equations had no appreciable effect on the predicted methane yield. It was expected that n-butyrate might be a better parameter than total butyrate to predict methane yield. This is because iso-butyrate is not synthesized by ruminal bacteria, but is derived from deamination of branched chain amino acids and its formation does not provide protons for methane production [62]. Preliminary testing of our models using either proportions of total butyrate or n-butyrate showed there were only minor differences in the prediction equations and no significant differences in the predicted methane yields (data not shown). This outcome was probably due to n-butyrate being 93% of the total butyrate in our development data set. We also note that while the production of iso-butyrate may not directly result in the production of protons that can be utilized to produce methane, iso-butyrate has been reported to stimulate ruminal cellulose digestion [63] and we speculate that this indirect mechanism may lead to greater methane production. Given there was no appreciable difference in the accuracy of prediction whether we used total butyrate or n-butyrate and that many scientific publications have reported butyrate without qualification as to type, we opted to use total butyrate in our prediction equations so they could be used to predict methane yields on data from any experiment.
Our yield of total VFA per kilogram of DMI was similar to previous work. An early experiment used a combination of in vitro and in vivo methods to estimate the production rate of total VFA [64]. Using data reported by Stewart et al. [64] we calculate that their yield of total VFA was, depending on the diet, between 3.0 and 3.3 mol total VFA/kg DMI. In the study by Stewart et al. [64] there were several unverifiable assumptions that were made regarding the production rates of total VFA and they admit that their “production figures presented earlier should be considered as minimal.” More recently, Lee et al. [65] fermented dried grass in a RUSITEC system and reported a yield of 3.87 mol of total VFA/kg DMI, a value that is much closer to our estimate. Overall, we consider our estimate of 3.98 mol total VFA/kg DMI to be plausible.
The variability in the observed relationships between methane yield and ruminal VFA proportions may be due to activation of different fermentation pathways, measurement error, or the sampling protocols. In addition to the fermentation pathways proposed by Moss et al. [13], there are others as outlined by Bryant [66], and these may come into effect differently when cows consume a wide range of diets or have diverse ruminal microbial populations [8]. We speculate that these diverse fermentation pathways might alter the relationship between ruminal VFA proportions and methane production. Measurement error is present regardless of the technique used to estimate methane yield [67]. All equipment used by us was calibrated prior to use and validated after use, and our SF6 technique has been shown to be comparable to respiration chambers [33]. However, there is still variation because of the factors that cannot be accounted for, such as the extent of the effect that the previous day’s DMI has on today’s methane emission [68]. Sampling protocol, particularly the duration between feeding and sampling of ruminal fluid has been shown to affect the proportions of VFA in ruminal fluid and hence influence the strength of the correlation between methane yield and proportion of VFA [23,45]. While there are differences in sampling periods and in the intervals between feeding and sampling of ruminal fluid there will always be unaccounted for variation between methane yield and individual VFA proportions in ruminal fluid. Alternatively, if multiple samples of ruminal fluid were collected from individual cows over the same period as the methane measurement, then we speculate that mean proportions of specific VFA applied to the models tested in this research would result in a dataset and equations with greater correlations to methane yield and smaller errors of prediction than those which we have reported. However, this would require further research. If precision is not required, for example if animals are only being ranked, then our best equations may be suitable without further improvement.
Predictors of, or proxies for, methane yield have tended to focus on DMI and diet composition [69]. Although such models have had some success, it seems unlikely that they can take account of ruminal microbes adapting over time to a particular diet and this could cause considerable errors in prediction [70]. Furthermore, prediction equations based on DMI and diet composition may not be applicable to grazing situations where it is difficult to accurately estimate DMI or the chemical composition of what animals actually select from a mixed pasture [71]. Milk has also been investigated as a predictor of methane yield. Dijkstra et al. [10] developed an equation to predict methane yield from milk fatty acids (R2 = 0.73, RMSEP not reported, concordance not reported) using 50 observations covering 10 diets. Furthermore, Niu et al. [8] used a database of over 3000 records and reported several equations to predict methane yield from feed and milk characteristics. The best equation of Niu et al. [8] used concentrations of dietary fat and neutral detergent fiber, energy corrected milk yield, milk fat concentration, milk protein concentration, and body weight (RMSEP = 14.7% of mean, concordance 0.37). However, for grazing cows the detailed composition of feed eaten is seldom known. When Niu et al. [8] used a simpler model based on energy corrected milk yield and milk composition, it resulted in similar predictive power to their more complicated equation (RMSEP = 15.1% of mean, concordance = 0.30). In comparison, our Equation (4) had R2 = 0.60, RMSEP = 15.4% of mean observed methane yield and a concordance of 0.73 between predicted and observed methane yield. This error is similar to that of Niu et al. [8] but the concordance of our equation is much greater. In addition, our equation does not require any knowledge of diet and can be used in non-lactating animals. Sampling ease and analysis cost must also be taken into consideration. Milk samples are easy to obtain but currently more expensive to analyze for fatty acids compared to the cost of analyzing ruminal fluid for VFA. The collection of ruminal fluid is a common procedure in many nutrition related experiments on ruminants to enable the measurement of a range of parameters in ruminal fluid. Thus, the additional measurement of VFA in ruminal fluid is not necessarily onerous, but would require the use of a standard protocol and skilled technicians.
Using VFA proportions to predict methane production as distinct from methane yield, has not been very successful to date. Propionate concentration has been correlated with methane production of 12 Merino-cross ewes, but explained only 26% of the variation in methane production [23]. Robinson et al. [23] concluded that “none of the suite of VFA parameters assessed offers a useful tool to predict methane production in grazing sheep.” More recently, a review on this topic concluded the relationships between methane production and ruminal VFA were variable and not as straight forward as theory suggested [3]. These pessimistic conclusions regarding the prediction of methane production from ruminants are only partly warranted. The DMI of individual animals may vary greatly, and DMI is the principal determinant of methane production in ruminants [4,72]. If DMI of an individual animal can be measured or accurately estimated and since methane yield can be predicted from VFA proportions in ruminal fluid, then methane production can be estimated as the product of DMI and methane yield. Therefore, contrary to the conclusions of Robinson et al. [23] and Negussie et al. [3], we consider that VFA parameters in ruminal fluid, when combined with measured or estimated DMI, do offer a potential tool to predict methane production in ruminants.
There are limitations when using ruminal VFA to predict methane yield or methane production. In comparison to control diets, diets that contain either nitrate, sulphate, 3-nitrooxypropanol, or halogenated methane analogues can cause substantial reductions in methane yield and methane production with little or no effect on ruminal VFA proportions [47,73,74]. For example, while Equation (4) can predict the methane yield of control diets with errors less than 9%, the prediction errors associated with diets containing sulphate were 21%, nitrate 40%, nitrate plus sulphate 85%, and 3-nitrooxypropanol 67%. Similar prediction errors may also occur when using proxies based on milk fatty acid profiles since 3-nitrooxypropanol substantially reduces the methane yield but has little or no effect on the fatty acid profile of milk [75]. Another limitation of our equations is that the minimum methane yield it can predict is given by the value of the constant term. This is especially problematic if methane yield is to be predicted for diets that contain halogenated methane analogues (e.g., red seaweed) that may depress methane yield to a value approaching zero [73,76].
Further work is necessary to better understand the relationships between methane yield and ruminal VFA. More frequent VFA determination, or daily VFA proportions collected in breath, may be one way of improving the relationships. The breath collection apparatus as described by Deighton et al. [33] has been shown to be able to collect a sample that represents the entire day’s output from breathing and eructation. Using breath samples would require a reconstruction of the prediction equation since reported proportions of VFA in breath [77] are quite different to those we observed in ruminal fluid.
Potential applications of the equations that we have presented include: (1) They could be used to predict methane yields of individuals or groups of grazing cattle, especially under extensive grazing situations; (2) the equations could be used under confinement feeding situations such as in TMR dairies and in beef feedlots to estimate the methane yields of a representative sub-group of animals. In these circumstances, the total feed intake of a group of animals is usually known; therefore, when estimated methane yields are combined with group intake, the equations allow the estimation of group methane emissions. This approach is likely to be the cheapest and most practical way to estimate methane emissions in these situations. (3) Our equations could be used in research to estimate methane yield from VFA data previously published in many scientific articles, and then meta-analyses could be conducted to test many theories. (4) Lastly, our equations could be used in large-scale screening programs to identify low methane emitting animals.

5. Conclusions

VFA proportions in ruminal fluid can be used to predict methane yields of individual or groups of dairy or beef cattle.

Author Contributions

Conceptualization, methodology, investigation, formal analysis, data curation, S.R.O.W., M.C.H., and P.J.M.; resources, J.L.J. and W.J.W.; writing–original draft and editing, S.R.O.W., M.C.H., and P.J.M.; review J.L.J. and W.J.W.; project administration, funding acquisition, J.L.J. and W.J.W.

Funding

This investigation was funded by the Victorian Department of Jobs, Precincts and Regions; Meat and Livestock Australia; Dairy Australia; and the Australian Government Department of Agriculture.

Acknowledgments

This work would not have been possible without the contributions of the science and technical staff at Agriculture Victoria Research, Ellinbank Centre, Victoria, Australia.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hammond, K.J.; Humphries, D.J.; Crompton, L.A.; Green, C.; Reynolds, C.K. Methane emissions from cattle: Estimates from short-term measurements using a GreenFeed system compared with measurements obtained using respiration chambers or sulphur hexafluoride tracer. Anim. Feed Sci. Technol. 2015, 203, 41–52. [Google Scholar] [CrossRef]
  2. Hegarty, R.S. Applicability of short-term emission measurements for on-farm quantification of enteric methane. Animal 2013, 7, 401–408. [Google Scholar] [CrossRef] [PubMed]
  3. Negussie, E.; de Haas, Y.; Dehareng, F.; Dewhurst, R.J.; Dijkstra, J.; Gengler, N.; Morgavi, D.P.; Soyeurt, H.; van Gastelen, S.; Yan, T.; et al. Invited review: Large-scale indirect measurements for enteric methane emissions in dairy cattle: A review of proxies and their potential for use in management and breeding decisions. J. Dairy Sci. 2017, 100, 2433–2453. [Google Scholar] [CrossRef] [PubMed]
  4. Charmley, E.; Williams, S.R.O.; Moate, P.J.; Hegarty, R.S.; Herd, R.M.; Oddy, V.H.; Reyenga, P.; Staunton, K.M.; Anderson, A.; Hannah, M.C. A universal equation to predict methane production of forage-fed cattle in Australia. Anim. Prod. Sci. 2016, 56, 169–180. [Google Scholar] [CrossRef]
  5. Moate, P.J.; Deighton, M.H.; Williams, S.R.O.; Pryce, J.E.; Hayes, B.J.; Jacobs, J.L.; Eckard, R.J.; Hannah, M.C.; Wales, W.J. Reducing the carbon footprint of Australian milk production by mitigation of enteric methane emissions. Anim. Prod. Sci. 2016, 56, 1017–1034. [Google Scholar] [CrossRef]
  6. Moate, P.J.; Williams, S.R.O.; Jacobs, J.L.; Hannah, M.C.; Beauchemin, K.A.; Eckard, R.J.; Wales, W.J. Wheat is more potent than corn or barley for dietary mitigation of enteric methane emissions from dairy cows. J. Dairy Sci. 2017, 100, 7139–7153. [Google Scholar] [CrossRef]
  7. Sauvant, D.; Giger-Reverdin, S.; Serment, A.; Broudiscou, L. Influences des régimes et de leur fermentation dans le rumen sur la production de méthane par les ruminants. INRA Prod. Anim. 2011, 24, 433–446. [Google Scholar]
  8. Niu, M.; Kebreab, E.; Hristov, A.N.; Oh, J.; Arndt, C.; Bannink, A.; Bayat, A.R.; Brito, A.F.; Boland, T.; Casper, D.; et al. Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database. Glob. Chang. Biol. 2018, 24, 3368–3389. [Google Scholar] [CrossRef]
  9. Chilliard, Y.; Martin, C.; Rouel, J.; Doreau, M. Milk fatty acids in dairy cows fed whole crude linseed, extruded linseed, or linseed oil, and their relationship with methane output. J. Dairy Sci. 2009, 92, 5199–5211. [Google Scholar] [CrossRef]
  10. Dijkstra, J.; van Zijderveld, S.M.; Apajalahti, J.A.; Bannink, A.; Gerrits, W.J.J.; Newbold, J.R.; Perdok, H.B.; Berends, H. Relationships between methane production and milk fatty acid profiles in dairy cattle. Anim. Feed Sci. Technol. 2011, 166–167, 590–595. [Google Scholar] [CrossRef]
  11. Williams, S.R.O.; Moate, P.J.; Deighton, M.H.; Hannah, M.C.; Wales, W.J. Methane emissions of dairy cows cannot be predicted by the concentrations of C8:0 and total C18 fatty acids in milk. Anim. Prod. Sci. 2014, 54, 1757–1761. [Google Scholar] [CrossRef]
  12. Demeyer, D. Quantitative aspects of microbial metabolism in the rumen and hind gut. In Rumen Microbial Metabolism and Ruminant Digestion; Jouany, J.P., Ed.; INRA Editions: Paris, France, 1991; pp. 217–237. [Google Scholar]
  13. Moss, A.R.; Jouany, J.-P.; Newbold, J. Methane production by ruminants: Its contribution to global warming. Ann. Zoot. 2000, 49, 231–253. [Google Scholar] [CrossRef]
  14. Holtshausen, L.; Chaves, A.V.; Beauchemin, K.A.; McGinn, S.M.; McAllister, T.; Odongo, N. Feeding saponin-containing Yucca schidigera and Quillaja saponaria to decrease enteric methane production in dairy cows. J. Dairy Sci. 2009, 92, 2809–2821. [Google Scholar] [CrossRef] [PubMed]
  15. Mohammed, N.; Ajisaka, N.; Lila, Z.A.; Hara, K.; Mikuni, K.; Kanda, S.; Itabashi, H. Effect of Japanese horseradish oil on methane production and ruminal fermentation in vitro and in steers. J. Anim. Sci. 2004, 82, 1839–1846. [Google Scholar] [CrossRef] [PubMed]
  16. Eugène, M.; Martin, C.; Mialon, M.M.; Krauss, D.; Renand, G.; Doreau, M. Dietary linseed and starch supplementation decreases methane production of fattening bulls. Anim. Feed Sci. Technol. 2011, 166–167, 330–337. [Google Scholar] [CrossRef]
  17. Danielsson, R.; Schnürer, A.; Arthurson, V.; Bertilsson, J. Methanogenic population and CH4 production in Swedish dairy cows fed different levels of forage. Appl. Environ. Microbiol. 2012, 78, 6172–6179. [Google Scholar] [CrossRef] [PubMed]
  18. Foley, P.A.; Kenny, D.A.; Callan, J.J.; Boland, T.M.; O’Mara, F.P. Effect of DL-malic acid supplementation on feed intake, methane emission, and rumen fermentation in beef cattle. J. Anim. Sci. 2009, 87, 1048–1057. [Google Scholar] [CrossRef] [PubMed]
  19. Mohammed, R.; McGinn, S.M.; Beauchemin, K.A. Prediction of enteric methane output from milk fatty acid concentrations and rumen fermentation parameters in dairy cows fed sunflower, flax, or canola seeds. J. Dairy Sci. 2011, 94, 6057–6068. [Google Scholar] [CrossRef]
  20. NHMRC. Australian Code for the Care and Use of Animals for Scientific Purposes, 7th ed.; National Health and Medical Research Council: Canberra, Australia, 2004.
  21. NHMRC. Australian Code for the Care and Use of Animals for Scientific Purposes, 8th ed.; National Health and Medical Research Council: Canberra, Australia, 2013.
  22. Bannink, A.; France, J.; Lopez, S.; Gerrits, W.J.J.; Kebreab, E.; Tamminga, S.; Dijkstra, J. Modelling the implications of feeding strategy on rumen fermentation and functioning of the rumen wall. Anim. Feed Sci. Technol. 2008, 143, 3–26. [Google Scholar] [CrossRef]
  23. Robinson, D.L.; Goopy, J.; Hegarty, R.S. Can rumen methane production be predicted from volatile fatty acid concentrations? Anim. Prod. Sci. 2010, 50, 630–636. [Google Scholar] [CrossRef]
  24. Williams, S.R.O.; Chaves, A.V.; Deighton, M.H.; Jacobs, J.L.; Hannah, M.C.; Ribaux, B.E.; Morris, G.L.; Wales, W.J.; Moate, P.J. Influence of feeding supplements of almond hulls and ensiled citrus pulp on the milk production, milk composition, and methane emissions of dairy cows. J. Dairy Sci. 2018, 101, 2072–2083. [Google Scholar] [CrossRef] [PubMed]
  25. Moate, P.J.; Williams, S.R.O.; Hannah, M.C.; Eckard, R.J.; Auldist, M.J.; Ribaux, B.E.; Jacobs, J.L.; Wales, W.J. Effects of feeding algal meal high in docosahexaenoic acid on feed intake, milk production, and methane emissions in dairy cows. J. Dairy Sci. 2013, 96, 3177–3188. [Google Scholar] [CrossRef] [PubMed]
  26. Williams, S.R.O.; Hannah, M.C.; Eckard, R.J.; Wales, W.J.; Moate, P.J. Supplementing the diet of dairy cows with fat or tannin reduces methane yield, and additively when fed in combination. Animal. Under review.
  27. Williams, S.R.O.; Moate, P.J.; Deighton, M.H.; Hannah, M.C.; Wales, W.J.; Jacobs, J.L. Milk production and composition, and methane emissions from dairy cows fed lucerne hay with forage brassica or chicory. Anim. Prod. Sci. 2016, 56, 304–311. [Google Scholar] [CrossRef]
  28. Moate, P.J.; Williams, S.R.O.; Deighton, M.H.; Hannah, M.C.; Ribaux, B.E.; Morris, G.L.; Jacobs, J.L.; Hill, J.; Wales, W.J. Effects of feeding wheat or corn and of rumen fistulation on milk production and methane emissions of dairy cows. Anim. Prod. Sci. 2019, 59, 891–905. [Google Scholar] [CrossRef]
  29. Moate, P.J.; Deighton, M.H.; Jacobs, J.L.; Ribaux, B.E.; Morris, G.L.; Hannah, M.C.; Mapleson, D.; Islam, M.S.; Wales, W.J.; Williams, S.R.O. Influence of proportion of wheat in a pasture-based diet, on milk yield, methane emissions, methane yield and ruminal protozoa of dairy cows. J. Dairy Sci. Under review.
  30. Moate, P.J.; (Agriculture Victoria Research, Ellinbank, Victoria, Australia); Williams, S.R.O.; (Agriculture Victoria Research, Ellinbank, Victoria, Australia). Personal communication, 2018.
  31. Tyrrell, H.F.; Reid, J.T. Prediction of the Energy Value of Cow’s Milk. J. Dairy Sci. 1965, 48, 1215–1223. [Google Scholar] [CrossRef]
  32. Williams, S.R.O.; Clarke, T.; Hannah, M.C.; Marett, L.C.; Moate, P.J.; Auldist, M.J.; Wales, W.J. Energy partitioning in herbage-fed dairy cows offered supplementary grain during an extended lactation. J. Dairy Sci. 2013, 96, 484–494. [Google Scholar] [CrossRef]
  33. Deighton, M.H.; Williams, S.R.O.; Hannah, M.C.; Eckard, R.J.; Boland, T.M.; Wales, W.J.; Moate, P.J. A modified sulphur hexafluoride tracer technique enables accurate determination of enteric methane emissions from ruminants. Anim. Feed Sci. Technol. 2014, 197, 47–63. [Google Scholar] [CrossRef]
  34. Geishauser, T. An instrument for collection and transfer of ruminal fluid and for administration of water soluble drugs in adult cattle. Bov. Pract. 1993, 27, 38–42. [Google Scholar]
  35. Moate, P.J.; Williams, S.R.O.; Torok, V.; Hannah, M.C.; Ribaux, B.E.; Tavendale, M.; Eckard, R.J.; Jacobs, J.L.; Auldist, M.J.; Wales, W.J. Grape marc reduces methane emissions when fed to dairy cows. J. Dairy Sci. 2014, 97, 5073–5087. [Google Scholar] [CrossRef] [PubMed]
  36. Beauchemin, K.A.; McGinn, S.M. Methane emissions from feedlot cattle fed barley or corn diets. J. Anim. Sci. 2005, 83, 653–661. [Google Scholar] [CrossRef] [PubMed]
  37. Beauchemin, K.A.; McGinn, S.M. Methane emissions from beef cattle: Effects of fumaric acid, essential oil, and canola oil. J. Anim. Sci. 2006, 84, 1489–1496. [Google Scholar] [CrossRef] [PubMed]
  38. Beauchemin, K.A.; McGinn, S.M.; Martinez, T.F.; McAllister, T.A. Use of condensed tannin extract from quebracho trees to reduce methane emissions from cattle. J. Anim. Sci. 2007, 85, 1990–1996. [Google Scholar] [CrossRef] [PubMed]
  39. Beauchemin, K.A.; McGinn, S.M.; Benchaar, C.; Holtshausen, L. Crushed sunflower, flax, or canola seeds in lactating dairy cow diets: Effects on methane production, rumen fermentation, and milk production. J. Dairy Sci. 2009, 92, 2118–2127. [Google Scholar] [CrossRef]
  40. Benchaar, C.; Hassanat, F.; Gervais, R.; Chouinard, P.Y.; Julien, C.; Petit, H.V.; Massé, D.I. Effects of increasing amounts of corn dried distillers grains with solubles in dairy cow diets on methane production, ruminal fermentation, digestion, N balance, and milk production. J. Dairy Sci. 2013, 96, 2413–2427. [Google Scholar] [CrossRef]
  41. Bougouin, A.; Ferlay, A.; Doreau, M.; Martin, C. Effects of carbohydrate type or bicarbonate addition to grass silage-based diets on enteric methane emissions and milk fatty acid composition in dairy cows. J. Dairy Sci. 2018, 101, 6085–6097. [Google Scholar] [CrossRef] [Green Version]
  42. Bougouin, A.; Martin, C.; Doreau, M.; Ferlay, A. Effects of starch-rich or lipid-supplemented diets that induce milk fat depression on rumen biohydrogenation of fatty acids and methanogenesis in lactating dairy cows. Animal 2018, 13, 1421–1431. [Google Scholar] [CrossRef]
  43. Guyader, J.; Eugène, M.; Meunier, B.; Doreau, M.; Morgavi, D.P.; Silberberg, M.; Rochette, Y.; Gerard, C.; Loncke, C.; Martin, C. Additive methane-mitigating effect between linseed oil and nitrate fed to cattle. J. Anim. Sci. 2015, 93, 3564–3577. [Google Scholar] [CrossRef]
  44. Hassanat, F.; Gervais, R.; Massé, D.I.; Petit, H.V.; Benchaar, C. Methane production, nutrient digestion, ruminal fermentation, N balance, and milk production of cows fed timothy silage- or alfalfa silage-based diets. J. Dairy Sci. 2014, 97, 6463–6474. [Google Scholar] [CrossRef]
  45. Hatew, B.; Podesta, S.C.; Van Laar, H.; Pellikaan, W.F.; Ellis, J.L.; Dijkstra, J.; Bannink, A. Effects of dietary starch content and rate of fermentation on methane production in lactating dairy cows. J. Dairy Sci. 2015, 98, 486–499. [Google Scholar] [CrossRef]
  46. Klop, G.; Dijkstra, J.; Dieho, K.; Hendriks, W.H.; Bannink, A. Enteric methane production in lactating dairy cows with continuous feeding of essential oils or rotational feeding of essential oils and lauric acid. J. Dairy Sci. 2017, 100, 3563–3575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Martinez-Fernandez, G.; Duval, S.; Kindermann, M.; Schirra, H.J.; Denman, S.E.; McSweeney, C.S. 3-NOP vs. halogenated compound: Methane production, ruminal fermentation and microbial community response in forage fed cattle. Front. Microbiol. 2018, 9, 1582. [Google Scholar] [CrossRef] [PubMed]
  48. McGinn, S.M.; Beauchemin, K.A.; Coates, T.; Colombatto, D. Methane emissions from beef cattle: Effects of monensin, sunflower oil, enzymes, yeast, and fumaric acid. J. Anim. Sci. 2004, 82, 3346–3356. [Google Scholar] [CrossRef] [PubMed]
  49. Olijhoek, D.W.; Løvendahl, P.; Lassen, J.; Hellwing, A.L.F.; Höglund, J.K.; Weisbjerg, M.R.; Noel, S.J.; McLean, F.; Højberg, O.; Lund, P. Methane production, rumen fermentation, and diet digestibility of Holstein and Jersey dairy cows being divergent in residual feed intake and fed at 2 forage-to-concentrate ratios. J. Dairy Sci. 2018, 101, 9926–9940. [Google Scholar] [CrossRef] [Green Version]
  50. Van Gastelen, S.; Antunes-Fernandes, E.C.; Hettinga, K.A.; Klop, G.; Alferink, S.J.J.; Hendriks, W.H.; Dijkstra, J. Enteric methane production, rumen volatile fatty acid concentrations, and milk fatty acid composition in lactating Holstein-Friesian cows fed grass silage- or corn silage-based diets. J. Dairy Sci. 2015, 98, 1915–1927. [Google Scholar] [CrossRef] [Green Version]
  51. Van Zijderveld, S.M.; Fonken, B.; Dijkstra, J.; Gerrits, W.J.J.; Perdok, H.B.; Fokkink, W.; Newbold, J.R. Effects of a combination of feed additives on methane production, diet digestibility, and animal performance in lactating dairy cows. J. Dairy Sci. 2011, 94, 1445–1454. [Google Scholar] [CrossRef] [Green Version]
  52. Hastie, T.; Tibshirani, R.; Friedman, J. The Elements of Statistical Learning; Springer: Berlin, Germany, 2012. [Google Scholar]
  53. Stevens, A.; Miralles, I.; Van Wesemael, B. Soil organic carbon predictions by airborne imaging spectroscopy: Comparing cross-validation and validation. Soil Sci. Soc. Am. J. 2012, 76, 2174–2183. [Google Scholar] [CrossRef]
  54. Lin, L.K. A concordance correlation coefficient to evaluate reproducibility. Biometrics 1989, 45, 255–268. [Google Scholar] [CrossRef]
  55. Bibby, J.; Toutenburg, H. Prediction and Improved Estimation in Linear Models; J. Wiley & Sons: Chichester, NY, USA, 1977. [Google Scholar]
  56. Gauch, H.G., Jr.; Hwang, J.T.G.; Fick, G.W. Model Evaluation by Comparison of Model-Based Predictions and Measured Values. Agron. J. 2003, 95, 1442–1446. [Google Scholar] [CrossRef] [Green Version]
  57. Lopes, J.C.; de Matos, L.F.; Harper, M.T.; Giallongo, F.; Oh, J.; Gruen, D.; Ono, S.; Kindermann, M.; Duval, S.; Hristov, A.N. Effect of 3-nitrooxipropanol on methane and hydrogen emissions, methane isotopic signature, and ruminal fermentation in dairy cows. J. Dairy Sci. 2016, 99, 5335–5344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Van Soest, P.J. Nutritional Ecology of the Ruminant; Cornell University Press: Ithaca, NY, USA, 1982; p. 373. [Google Scholar]
  59. Herd, R.; Bird, S.; Donoghue, K.; Arthur, P.; Hegarty, R. Phenotypic associations between methane production traits, volatile fatty acids and animal breeding traits. In Proceedings of the 20th conference of the Association for the Advancement of Animal Breeding and Genetics, Napier, New Zealand, 20–23 October 2013; pp. 286–289. [Google Scholar]
  60. Smith, H.C.; Herd, R.M.; Donoghue, K.; Bird-Gardiner, T.; Arthur, P.F.; Hegarty, R.F. Associations of rumen volatile fatty acids with phenotypic and genetic variation in methane production traits in Angus cattle. Proc. Assoc. Adv. Anim. Breed. Genet. 2015, 21, 93–96. [Google Scholar]
  61. Alemu, A.W.; Dijkstra, J.; Bannink, A.; France, J.; Kebreab, E. Rumen stoichiometric models and their contribution and challenges in predicting enteric methane production. Anim. Feed Sci. Technol. 2011, 166–167, 761–778. [Google Scholar] [CrossRef]
  62. Hoffmann, E.M.; Selje-Assmann, N.; Becker, K. Dose studies on anti-proteolytic effects of a methanol extract from Knautia arvensis on in vitro ruminal fermentation. Anim. Feed Sci. Technol. 2008, 145, 285–301. [Google Scholar] [CrossRef]
  63. Bentley, O.G.; Johnson, R.R.; Hershberger, T.V.; Cline, J.H.; Moxon, A.L. Cellulolytic-factor activity of certain short-chain fatty acids for rumen microorganisms in vitro. J. Nutr. 1955, 57, 389–400. [Google Scholar] [CrossRef] [PubMed]
  64. Stewart, W.E.; Stewart, D.G.; Schultz, L.H. Rates of volatile fatty acid production in the bovine rumen. J. Anim. Sci. 1958, 17, 723–736. [Google Scholar] [CrossRef]
  65. Lee, M.R.F.; Merry, R.J.; Davies, D.R.; Moorby, J.M.; Humphreys, M.O.; Theodorou, M.K.; MacRae, J.C.; Scollan, N.D. Effect of increasing availability of water-soluble carbohydrates on in vitro rumen fermentation. Anim. Feed Sci. Technol. 2003, 104, 59–70. [Google Scholar] [CrossRef]
  66. Bryant, M.P. Microbial methane production—Theoretical aspects. J. Anim. Sci. 1979, 48, 193–201. [Google Scholar] [CrossRef]
  67. Hristov, A.N.; Kebreab, E.; Niu, M.; Oh, J.; Bannink, A.; Bayat, A.R.; Boland, T.B.; Brito, A.F.; Casper, D.P.; Crompton, L.A.; et al. Symposium review: Uncertainties in enteric methane inventories, measurement techniques, and prediction models. J. Dairy Sci. 2018, 101, 6655–6674. [Google Scholar] [CrossRef]
  68. Moate, P.J.; Deighton, M.H.; Hannah, M.C.; Williams, S.R.O. Intake effects on methane emissions from dairy cows. In Proceedings of the CCRSPI Conference, Melbourne, Australia, 27–29 November 2012; pp. 60–61. [Google Scholar]
  69. Ramin, M.; Huhtanen, P. Development of equations for predicting methane emissions from ruminants. J. Dairy Sci. 2013, 96, 2476–2493. [Google Scholar] [CrossRef]
  70. Moate, P.J.; Jacobs, J.L.; Hannah, M.C.; Morris, G.L.; Beauchemin, K.A.; Alvarez Hess, P.S.; Eckard, R.J.; Liu, Z.; Rochfort, S.; Wales, W.J.; et al. Adaptation responses in milk fat yield and methane emissions of dairy cows when wheat was included in their diet for 16 weeks. J. Dairy Sci. 2018, 101, 7117–7132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Reeves, M.; Fulkerson, W.J.; Kellaway, R.C.; Dove, H. A comparison of three techniques to determine the herbage intake of dairy cows grazing kikuyu (Pennisetum clandestinum) pasture. Aust. J. Exp. Agric. 1996, 36, 23–30. [Google Scholar] [CrossRef]
  72. Swainson, N.; Muetzel, S.; Clark, H. Updated predictions of enteric methane emissions from sheep suitable for use in the New Zealand national greenhouse gas inventory. Anim. Prod. Sci. 2018, 58, 973–979. [Google Scholar] [CrossRef]
  73. Li, X.; Norman, H.C.; Kinley, R.D.; Laurence, M.; Wilmot, M.; Bender, H.; De Nys, R.; Tomkins, N. Asparagopsis taxiformis decreases enteric methane production from sheep. Anim. Prod. Sci. 2018, 58, 681–688. [Google Scholar] [CrossRef]
  74. van Zijderveld, S.M.; Gerrits, W.J.J.; Apajalahti, J.A.; Newbold, J.R.; Dijkstra, J.; Leng, R.A.; Perdok, H.B. Nitrate and sulfate: Effective alternative hydrogen sinks for mitigation of ruminal methane production in sheep. J. Dairy Sci. 2010, 93, 5856–5866. [Google Scholar] [CrossRef] [Green Version]
  75. Reynolds, C.K.; Humphries, D.J.; Kirton, P.; Kindermann, M.; Duval, S.; Steinberg, W. Effects of 3-nitrooxypropanol on methane emission, digestion, and energy and nitrogen balance of lactating dairy cows. J. Dairy Sci. 2014, 97, 3777–3789. [Google Scholar] [CrossRef]
  76. Kinley, R.D.; De Nys, R.; Vucko, M.J.; MacHado, L.; Tomkins, N.W. The red macroalgae Asparagopsis taxiformis is a potent natural antimethanogenic that reduces methane production during in vitro fermentation with rumen fluid. Anim. Prod. Sci. 2016, 56, 282–289. [Google Scholar] [CrossRef]
  77. Oertel, P.; Küntzel, A.; Reinhold, P.; Köhler, H.; Schubert, J.K.; Kolb, J.; Miekisch, W. Continuous real-time breath analysis in ruminants: Effect of eructation on exhaled VOC profiles. J. Breath Res. 2018, 12. [Google Scholar] [CrossRef]
Figure 1. Methane yield (g/kg DMI) versus ruminal acetate, propionate and butyrate as mol per 100 mol of total volatile fatty acids from 215 cow records across 24 diets [24,25,26,27,28,29,30]. Solid lines are linear or reciprocal trend lines.
Figure 1. Methane yield (g/kg DMI) versus ruminal acetate, propionate and butyrate as mol per 100 mol of total volatile fatty acids from 215 cow records across 24 diets [24,25,26,27,28,29,30]. Solid lines are linear or reciprocal trend lines.
Animals 09 01006 g001
Figure 2. Combined methane yield (MY, g/kg DM) data from seven previously conducted experiments [24,25,26,27,28,29,30] plotted against ruminal fatty acid ratio A/P, (A + B)/P, and 1/P, with linear equations fitted by linear mixed effects meta-analysis (solid line).
Figure 2. Combined methane yield (MY, g/kg DM) data from seven previously conducted experiments [24,25,26,27,28,29,30] plotted against ruminal fatty acid ratio A/P, (A + B)/P, and 1/P, with linear equations fitted by linear mixed effects meta-analysis (solid line).
Animals 09 01006 g002
Figure 3. Independent validation of (a) Equation (3), (b) Equation (4), and (c) Equation (7) showing the observed methane yield against equation-predicted methane yield, for 32 individual [6], (open circles) and 67 treatment means from 16 previously published experiments [14,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51] (closed circles) with 1:1 line of agreement (dashed).
Figure 3. Independent validation of (a) Equation (3), (b) Equation (4), and (c) Equation (7) showing the observed methane yield against equation-predicted methane yield, for 32 individual [6], (open circles) and 67 treatment means from 16 previously published experiments [14,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51] (closed circles) with 1:1 line of agreement (dashed).
Animals 09 01006 g003aAnimals 09 01006 g003b
Table 1. Abbreviations, their descriptions and units for variables and parameters used in the equations.
Table 1. Abbreviations, their descriptions and units for variables and parameters used in the equations.
AbbreviationDescriptionUnits
Aacetate proportionmol/100 mol total VFA
Btotal butyrate proportionmol/100 mol total VFA
C2acetate quantitymol
C3propionate quantitymol
C4total butyrate quantitymol
CH4Pmethane production, stoichiometrymol
CH4Ymethane yield, stoichiometrymol/100 mol total VFA
DMIdry matter intakekg/d
MYmethane yield, predictedg methane/kg DMI
Ppropionate proportionmol/100 mol total VFA
[P]propionate concentrationmmol/L
VFAVolatile fatty acidsmol
Table 2. Details of the seven experiments and means for 24 treatments from which individual cow data were used to develop the prediction equations.
Table 2. Details of the seven experiments and means for 24 treatments from which individual cow data were used to develop the prediction equations.
ExptReferenceExpt Design 1Dietary Treatmentn CowsMethane MethodDIM Diet BaseSeasonMilk YieldECM 2Forage DMIConcentrate DMITotal DMI
(kg/Day)(kg/Day)(kg/Day)(kg/Day)(kg/Day)
1[24]Randomized blockControl (lucerne)12SF6179LucerneAutumn25.626.414.28.222.4
Almond hulls1017723.224.614.48.222.6
Citrus pulp1017625.925.413.27.821.0
2[25] Randomized blockDHA-0 g8Calorimeter 218LucerneAutumn22.225.418.25.924.1
DHA-25 g721526.024.018.06.124.1
DHA-50 g821623.121.916.56.222.7
DHA-75 g721522.321.615.36.221.5
3[26]Full crossoverControl (water)8Calorimeter90LucerneSpring32.332.319.15.824.9
Fat79534.532.017.85.923.7
Fat Plus Tannin89033.731.317.95.923.8
Tannin89031.130.018.75.924.6
4[27]Randomized blockControl (lucerne)12SF6149LucerneSummer23.922.715.45.420.8
Forage brassica1015027.525.415.25.420.6
Chicory1014520.419.312.35.417.7
5[28]Full crossoverCorn14Calorimeter191LucerneAutumn27.428.610.112.222.3
Wheat1419129.123.89.011.520.5
6[29]Randomized blockWheat-0 kg8SF657PastureSpring29.929.517.12.119.2
Wheat-3 kg85731.332.415.45.020.4
Wheat-6 kg85732.333.012.37.920.2
Wheat-9 kg85734.732.98.910.919.8
7[30]Randomized blockCorn 2 times/d8SF6206LucerneAutumn21.220.68.710.218.9
Wheat 2 times/d815721.322.08.810.219.0
Wheat 6 times/d816124.022.68.710.118.8
Wheat + buffer, 2 times/d819720.522.48.810.419.2
1 Expt = Experiment, DIM = days in milk, DMI = dry matter intake, DHA = docosahexanoic acid, control diet is unique within experiment. 2 ECM = energy corrected milk as per Tyrrell and Reid [31].
Table 3. Summary of data used to develop models between methane yield and individual volatile fatty acids.
Table 3. Summary of data used to develop models between methane yield and individual volatile fatty acids.
ParameterMethane Yield (g CH4/kg DMI)Acetate (M%)Propionate (M%)Butyrate (M%)
Characteristics
Mean (±s.d.)20.7 (4.81)63.3 (4.53)20.2 (5.44)9.9 (4.02)
Minimum6.951.312.90.7
Maximum32.470.740.215.6
Table 4. Summary of the data used to externally validate the prediction equations.
Table 4. Summary of the data used to externally validate the prediction equations.
ExperimentExperiment DesignN 1AnimalsMethane MethodDays in Milk Diet BaseMethane Yield AcetatePropionateButyrate
(g/kg DM)(M%)(M%)(M%)
[6]Randomized block32Dairy cowsSF671Lucerne hay10.1–26.951.5–67.416.3–38.56.5–15.3
[14]Full crossover3Dairy cowsCalorimeter155TMR, barley silage18.5–19.261.5–61.722.8–23.111.5–11.5
[36]Randomized block4Beef heifersCalorimeter-TMR, barley9.2–24.842.6–64.020.5–45.77.1–12.6
[37]Latin square4Beef heifersCalorimeter-Barley silage19.9–21.659.0–65.020.6–25.99.2–11.4
[38]Latin square3Beef cattleCalorimeter-Barley silage18.5–18.864.6–65.819.0–19.611.5–11.6
[39]Crossover4Dairy cowsCalorimeter96Barley silage and grain13.4–16.359.9–61.725.1–26.78.7–10.3
[40]Latin square4Dairy cowsCalorimeter99Alfalfa and corn silage18.9–20.660.1–63.421.8–23.112.3–14.4
[41]Latin square4Dairy cowsCalorimeter61Grass silage17.3–21.369.5–70.415.1–18.68.1–11.8
[42]Latin square4Dairy cowsCalorimeter61Corn silage18.7–20.264.7–66.118.0–19.810.4–12.8
[43]Randomized block2Dairy cowsCalorimeter-Grass hay20.7–25.067.0–70.2158–19.410.1–10.4
[44]Latin square3Dairy cowsCalorimeter92Silage19.7–20.166.6–67.819.2–19.510.2–11.0
[45]Randomized block4Dairy cowsCalorimeter215TMR, grass silage21.5–22.467.9–69.115.6–16.811.1–11.8
[46]Randomized block8Dairy cowsCalorimeter139TMR, corn19.7–23.164.2–66.718.0–20.811.5–12.4
[47]Randomized block2Beef steersCalorimeter-Grass hay22.9–23.270.0–73.516.1–16.67.6–10.2
[48]Latin square4Beef steersCalorimeter-TMR, barley silage18.8–22.664.0–68.018.3–21.68.6–10.4
[49]Crossover8Dairy cowsCalorimeter187Grass clover silage15.4–23.460.2–69.715.4–26.59.3–13.5
[50]Randomized block4Dairy cowsCalorimeter192Grass or corn silage22.0–25.063.6–66.017.1–18.912.8–16.3
[51]Randomized block2Dairy cowsCalorimeter176TMR, grass silage20.5–22.158.6–60.319.8–23.013.6–15.3
1 n is number of animals for [6] and number of treatments for all other sources.
Table 5. Seven models, used to predict methane yield (MY, g methane/kg dry matter intake) from ruminal acetate (A), propionate (P), and butyrate (B) expressed as mol/100 mol total volatile fatty acid, and the concentration of propionate in ruminal fluid ([P], mmol/L). Their coefficients 1 and constants 2 were estimated by fitting the models to training data, pooled from seven experiments by meta-analysis. Root-mean-square error of prediction (RMSEP, g CH4/kg dry matter intake) and Lin’s concordance (CCC) were calculated by a leave-one-experiment-out cross-validation procedure. These were used to calculate rank (1 = minimum RMSEP, maximum CCC) of equations tested.
Table 5. Seven models, used to predict methane yield (MY, g methane/kg dry matter intake) from ruminal acetate (A), propionate (P), and butyrate (B) expressed as mol/100 mol total volatile fatty acid, and the concentration of propionate in ruminal fluid ([P], mmol/L). Their coefficients 1 and constants 2 were estimated by fitting the models to training data, pooled from seven experiments by meta-analysis. Root-mean-square error of prediction (RMSEP, g CH4/kg dry matter intake) and Lin’s concordance (CCC) were calculated by a leave-one-experiment-out cross-validation procedure. These were used to calculate rank (1 = minimum RMSEP, maximum CCC) of equations tested.
Equation NumberModelEstimates ± S.E.RMSEPCCCRank
(1c)MY = 16 × (0.50A − 0.25P + 0.50B) × c/100c = 3.98 ± 0.153.70.524
(2)MY = dA − eP + fBd = 0.30 ± 0.042, e = 0.22 ± 0.055, f = 0.48 ± 0.163.60.634
(3)MY = g(A/P) + hg = 4.08 ± 0.36, h = 7.05 ± 1.403.20.691
(4)MY = i(A+B)/P + ji = 3.28 ± 0.29. j = 7.60 ± 1.283.20.701
(5)MY = kP + mk =−0.57 ± 0.057, m = 32.3 ± 1.43.50.634
(6)MY = n[P] + qn =−0.24 ± 0.035, q = 26.1 ± 1.34.20.367
(7)MY = s/P + ts = 316 ± 28, t = 4.4 ± 1.53.20.701
1 Coefficient c has units of moles of total VFA/kg DMI, coefficients d, e, f, k, and q have units of g CH4·kg DMI−1·moles−1·100 moles total VFA, the coefficients g, and i have units of g CH4·kg DMI−1 and the coefficient m has units of g CH4·kg DMI−1·L·mmol−1. 2 Constants h, j, l, n, and s have the units of g CH4·kg DMI−1.

Share and Cite

MDPI and ACS Style

Williams, S.R.O.; Hannah, M.C.; Jacobs, J.L.; Wales, W.J.; Moate, P.J. Volatile Fatty Acids in Ruminal Fluid Can Be Used to Predict Methane Yield of Dairy Cows. Animals 2019, 9, 1006. https://doi.org/10.3390/ani9121006

AMA Style

Williams SRO, Hannah MC, Jacobs JL, Wales WJ, Moate PJ. Volatile Fatty Acids in Ruminal Fluid Can Be Used to Predict Methane Yield of Dairy Cows. Animals. 2019; 9(12):1006. https://doi.org/10.3390/ani9121006

Chicago/Turabian Style

Williams, S. Richard O., Murray. C. Hannah, Joe L. Jacobs, William J. Wales, and Peter J. Moate. 2019. "Volatile Fatty Acids in Ruminal Fluid Can Be Used to Predict Methane Yield of Dairy Cows" Animals 9, no. 12: 1006. https://doi.org/10.3390/ani9121006

APA Style

Williams, S. R. O., Hannah, M. C., Jacobs, J. L., Wales, W. J., & Moate, P. J. (2019). Volatile Fatty Acids in Ruminal Fluid Can Be Used to Predict Methane Yield of Dairy Cows. Animals, 9(12), 1006. https://doi.org/10.3390/ani9121006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop