The greatest determinants of enteric CH4 production in ruminants are the quantity and digestibility of feed ingested. The relationship between feed intake and daily CH4 production rate (MPR) has been determined experimentally over many years down to the present(Reference Kriss1–Reference Warner, Bannink and Hatew6). Whilst it has been found that there is variability in MPR between animals fed the same ration, and variability between days for the same animal(Reference Blaxter and Clapperton7), it has also been demonstrated that there is a strong positive, linear relationship between MPR and the level of intake of a mainly forage diet that holds generally constant over a large range of intakes and types of forages(Reference Charmley, Williams and Moate8).
Sub-Saharan Africa has a ruminant livestock sector which is responsible for a disproportionally large level of anthropogenic greenhouse gas emissions in economies that are dominated by agriculture(Reference Tubiello, Salvatore and Cóndor Golec9), and there is a pressing need to comprehensively and accurately estimate these emissions. A key plank in the development of emission factors is the employment of accurate CH4 conversion factors (Y m: CH4 (MJ)/gross energy intake (MJ)). Equations developed by Charmley et al. (Reference Charmley, Williams and Moate8) have been derived from measurements of both Bos taurus and Bos indicus cattle consuming a wide range of tropical and temperate forages and thus may be considered more applicable to African cattle than equations derived solely from temperate regions. However, a significant limitation to the utility of these equations is that measurements have been conducted only on animals fed for production (i.e. at maintenance and above), whereas it has been recently demonstrated that cattle regularly experience episodes of significant seasonal weight loss in African smallholder systems(Reference Goopy, Onyango and Dickhoefer10,Reference Ndung’u, Bebe and Ondiek11) .
It has been established that higher CH4 yield (MY: g CH4/kg DM intake) is associated with longer mean rumen retention time (MRT) in sheep(Reference Pinares-Patiño, Ulyatt and Lassey12,Reference Goopy, Donaldson and Hegarty13) , while both feed digestibility and MRT increase with declining intake in limit-fed cattle(Reference Shaver, Nytes and Satter14). In contrast, the relationship between intake and MY in animals consuming feed at levels below maintenance requirements is not well understood and there is circumstantial evidence to suggest that the relationship may be altered when intake is very low. Several studies have demonstrated that MRT increases in ruminants fed food of a given quality well below maintenance, while apparent digestibility remains unchanged, or in some cases, actually declines(Reference Grimaud, Richard and Kanwé15–Reference Doreau, Michalet-Doreau and Grimaud18). Increasing the time digesta spends within the rumen may expose it to a kind of futile cycling, where material which may otherwise be enzymatically digested by the host is subjected to further microbial degradation in the rumen, resulting in higher production of volatile fatty acids, H+ and CH4, at the expense of microbial protein. Refermentation of digesta has been observed in sheep fed specialised diets(Reference Rowe, Loughnan and Nolan19) and has been posited as a mechanism where sheep fed low-quality diets produced more CH4 than predicted(Reference Goopy, Nolan and Pelaelo-Grand20).
At present, the Intergovernmental Panel on Climate Change guidelines recommend adopting a Y m of 6·3 % to estimate enteric emissions from cattle consuming a tropical forage-based diet, which has been substantially corroborated by the extensive study of Charmley et al. (Reference Charmley, Williams and Moate8). However, if below-maintenance intake results in a higher than expected MPR, this will have significant implications, both for inventory and for intervention options in systems where ruminants often can only access below maintenance rations for part of the year – as in Sub-Saharan Africa. Thus, this study investigated the relationship between feed intake, rumen kinetics, apparent digestibility and CH4 production in animals consuming diets at levels well below maintenance requirements. We hypothesised that MY and thus Y m would increase with increasing severity of sub-maintenance feeding.
Materials and methods
Animals and experimental design
All animal procedures were carried out adhering to the international standards for animal care and use for scientific purposes, reviewed by the Institutional Animal Use and Care Committee of the International Livestock Research Institute permit no: IACUC-RC2016-11.
Boran (B. indicus) yearling steers (n 12; live weight (LW): 162·3 (sem 3·77) kg) were sourced from a commercial ranch in Lakipia County (Northern Kenya). Before the commencement of the trial, animals were treated with an anthelmintic, an acaricide wash, vaccinated for foot and mouth disease and clostridial diseases, ear tagged and placed under quarantine for 21 d as part of the standard induction procedure of the institution. Steers were housed in open individual pens (1·90 × 2·87 m), covered with shade-cloth sails, with clean water supplied ad libitum from automatic waterers.
The experimental design was a 4 × 4 Latin square with four levels of feeding (maintenance energy requirement (MER): 1·0 MER, 0·8 MER, 0·6 MER, 0·4 MER). Steers were stratified by LW and randomly allocated to the four treatment groups (three animals per group, with the individual animal being the experimental unit). The groups were maintained throughout the experimental period and were randomly allocated to one of the four feeding levels in the first period (Table 1).
* Groups each with three Boran steers.
Diets and feeding
Experimental diets were based on an allocation of chaffed Rhodes grass late-cut hay (Chloris gayana cv. Boma DM: 875 g/kg; digestible energy: 8·4 MJ/kg DM; crude protein (CP): 73·1 g/kg DM), plus the addition of a small amount (equal to 10 % each of the calculated energy content of the ration) of cotton seed meal (DM: 947 g/kg; digestible energy: 12·7 MJ/kg DM; CP: 324·4 g/kg DM) and molasses (DM: 728 g/kg; digestible energy: 14·2 MJ/kg DM; CP: 46 g/kg DM) to the ration of animals being fed at 1·0 MER (in order to achieve required intake). Initially, MER was estimated using a regression equation developed from maintenance requirements for non-lactating dairy cattle(21), and rations were formulated based on the proximate analysis of feed samples. The rations were subsequently reduced by 20 % after observing the animal response during pre-trial feeding; thus, the rations were set based on the following equation:
Treatment periods were 35 d, including an adaption period of 14 d prior to each measurement period. During the final 14 d of each period, total urine and faecal collection were undertaken and three, 22·5 h measurements of enteric CH4 production for each animal in every period were conducted on alternate days in open-circuit respiratory chambers. Because most animals were expected to lose weight during the experimental periods, each measurement period was followed by a 2-week ‘feed-up’ period where animals received 1 kg each of cotton seed meal and molasses in days 1–3, then 1·5 kg of each for days 4–14, plus ad libitum hay. Rations for animals in each experimental period were set on their LW on the first day of each experimental period, based on the animals’ estimated MER (eq. 1) and proximate analysis of the feeds as outlined above.
Animals were fed twice daily to decrease loss of hay through spillage (09.30 hours, after removal and weighing of orts and 14.00 hours), with hay placed in the feeding trough, and the supplement (1·0 MER only) in a separate bucket.
Sample collection and analysis
Feed and refusals
After weighing, refusals were bulked daily by treatment, homogenised and a subsample (approximately 200 g) taken and stored in zipped polythene bags at −20°C. Samples of the basal diet were taken every 2 weeks and stored as for the refusals. At the end of each treatment period, the basal diet and refusals were bulked by treatment, homogenised and subsamples retained for subsequent processing and analysis.
Total faecal excretion was determined by daily collection over 6 d. Total faeces were weighed, and a sub-sample of approximately 500 g was transferred into labelled foil trays and then dried in a drying oven at 50°C until a constant weight was obtained for at least two consecutive days. Samples were cooled in desiccators, final weight determined, then stored in zipped polythene bags at room temperature until further analysis.
Dried faeces, feed and refusals were ground in a hammer mill through a 1 mm sieve and analysed as follows: true DM was determined at 105°C for 24 h; ash was determined by combustion in a muffle furnace at 550°C according to the methods of the AOAC (AOAC, 1990 methods no. 924.05). Feed and refusals were analysed for neutral-detergent fibre and acid-detergent fibre by the methods of Van Soest(Reference Van Soest, Robertson and Lewis22) using an Ankom200 fibre analyzer (ANKOM Technology), with alfa-amylase enzyme. Total N content in feed and faeces was determined by the micro Kjeldahl procedure of AOAC (AOAC, 1990, method no. 988.05) with Se catalyst tablets. Gross energy content of feed was determined by bomb caliometry (Parr 6300, Parr Instruments). The digestible energy content of original fed samples was estimated from in vitro organic matter digestibility from the equation from the National Research Council(24).
Marker application procedures, analysis and calculations
Passage of liquid and solid digesta through the gastrointestinal tract was determined using Co-EDTA(Reference Udén, Colucci and Van Soest25) and Yb-marked fibre particles(Reference Teeter, Owens and Mader26), respectively. At commencement of passage rate determination, individual animals were offered a quantity of Yb-marked fibre equivalent to 560 mg Yb/kg LW(Reference Richter and Schlecht27) mixed with 20 g molasses before morning feeding. After consumption of the marked fibre, each animal was drenched with Co-EDTA (23·56 mg Co/kg-LW). Time zero (t0) was taken as the completion of the Co-EDTA drench. To determine Yb and Co concentrations in faeces, about 60 g fresh faeces were sampled at 0 h (t0, then two hourly from 4 to 16 h after administration, four hourly from 16 to 40 h, six hourly from 40 to 88 h, then eight hourly from 88 to 136 h and finally at 148 h). Dried samples underwent sealed chamber digestion(Reference Anderson and Henderson28), and Yb and Co concentrations were determined by inductively coupled plasma optical emission spectroscopy (ICP-OES 5100 VDV, Agilent Technologies). The cumulative quantity of Yb and Co excreted during the total collection period (148 h) was determined as the concentration of the respective elements in individual faecal samples multiplied by total faecal mass at time ti (sampling time). PROC NLIN (method = dud) in SAS 9.1 (SAS Institute Inc.) was applied using the Type-N one-compartment Gamma-2 model of Richter & Schlecht(Reference Richter and Schlecht27) for parameters of liquid (Co) and solid (Yb) digesta passage as follows: first appearance of markers (TT; equivalent to post-ruminal laminar flow) and ruminal passage rate (λ ruminal retention time (2λ−1).
Enteric methane measurement and respiratory chamber operation
MPR was measured over three periods of approximately 22·5 h each, with each measurement period separated by 1 d; thus, six steers were assessed for MPR in each of weeks four and five in each experimental period (with total collection taking place in the alternate week) (see Fig. 1). Measurement was commenced at 09.30 hours with 1·5 h allocated for cleaning and unloading/loading animals. Thus, measurements were taken for all the steers over the same 22·5 h. On measurement days, immediately after refusals were collected, steers were taken from their pens to individual open-circuit respiratory chambers. The ration for the day was placed in the feed bin inside the chamber and the doors shut and sealed. Feed refusals were measured at the conclusion of each period, and steers were returned to their pens.
Respiration chambers, measuring equipment and operating conditions used in this study are presented here for the first time. Each custom-built chamber (3·08 m length × 1·50 m width × 2·00 m height; internal volume 8·90 m3) (No Pollution Industries) was constructed of injected polyurethane foam sandwich panels (60 mm thick), with external and internal lining of AISI 304 stainless steel. Two chambers shared a common wall, the third was independent. Large double-glazed, laminated glass panels were installed in side walls to provide visual communication between animals, with separate doors for entry and exit at opposite ends of each chamber also with windows fitted. Internally, chambers were equipped with supply and return ventilation grills, installed at ceiling height to mix the recirculation air and provide a slow movement of air across the chamber. The internal environment of each chamber was controlled by an air handling unit, treating a maximum volume of about 700 m3/h. Each air handling unit was equipped with a mixing, a filtering and a treatment section, composed of a cooling coil, a heating pack and a humidification circuit. A frequency-controlled main fan recirculates the air through the chamber and the air handling unit. Cold water from a chiller unit (MICS 0092 FF; Climaveneta Co.) passing through the cooling coil cooled and removed moisture from the air. The heating pack was electrically powered, and extraction of air from each chamber was achieved by a high-pressure blower (model: SCL K04-MS MOR, FPZ SpA), powered through a variable speed drive to allow control of the volume of air leaving each chamber, which was measured at the extraction point by a differential pressure system (model DP 2500-R8-AZ, Johnson controls Inc.), together with temperature and relative humidity. Overall control of the chambers, monitoring of conditions and events were carried out by custom designed software (Johnson controls metasys software, version 2.6, Johnson Controls Inc.), which integrates downloading capabilities for the main research parameters.
For the duration of the trial, the internal environment was set at 22°C, 50 % relative humidity and internal air circulation of 80 litres/s. Air was taken from the chambers under negative pressure by a pump at a flow of approximately 18 litres/s. MPR was calculated as chamber air flow multiplied by CH4 concentration in the chamber adjusted for CH4 concentration of the incoming air and temperature and atmospheric pressure in the chamber. The 22·5 h value was converted to MPR by multiplying by 24/22·5. Actual air flow through each chamber was measured using a venturi apparatus with differential pressure transducers (model DP 2500-R8-AZ, Johnson Controls Inc.).
Concentration of CH4 (parts per million per volume) was measured in the chamber incoming and exhaust air streams using a cavity ringdown laser absorption spectrometer (Picarro G2508 analyzer). The laser was validated using known CH4 concentrations. Sensors for relative humidity (model EE160 HCT01-00D, E+E Elektronic Ges.m.b.H, Elektronik Ges.m.b.H), CO2 (model EE85-10C35 E+E Elektronik Ges.m.b.H) and ambient temperature (model A99B-500C, Johnson Controls Inc.) monitored conditions in all chambers and recorded using the Sensing Science Laboratory software (Data Harvest Group) on a separate personal computer. Recoveries of CH4 were carried out for each chamber at the end of every experimental period by injecting known amounts of CH4 into the chamber using a gas-phase titration unit and measuring the total amount of CH4 found in the exhaust air.
Statistical analysis
Firstly, the effects of treatment (intake level) on intake, digestibility, CH4 production and rumen kinetics were analysed using R 3.0.3 (R development core team). Treatment and period effects were compared using ANOVA type 3 using a linear mixed model fitted by restricted maximum likelihood (REML) t tests with Satterthwaite approximations of df (LmerTest) were used. Differences between means were compared using Tukey’s honestly significant difference, and level of significance was determined at 0·05. We also explored relationships in the data using a mixed linear model (lme4 in R(Reference Bates, Mächler and Bolker29)) using either a linear (Y = a + bX) or two-factor polynomial (Y = a(X)2 + bX + c) model for treatments with period and animal ID included as random variables; both models used maximum likelihoods with Satterthwaite approximations to df. The ‘best’ model was selected after examination of the χ 2 values as well as the Akaike information criterion(Reference Burnham, Anderson and Huyvaert30), in order to assess if the extra complexity explained sufficient variability to justify its use. All linear and quadratic relationships are provided in online Supplementary Tables S1–S3.
Sample size justification
We hypothesised that digestibility would vary between the control (100 % MER) and the other treatment studied by 7–10 % of the overall mean. We then used a figure of 6·1 % in our calculations (from Doreau et al. (Reference Doreau, Michalet-Doreau and Béchet31), with a standard deviation of 4·4 % of the mean. Using a standard two-tail t test comparison, 80 % statistical power, and level of significance set at 5 %, a total of nine animals per treatment was arrived at based on the formula described by Charan & Kantharia(Reference Charan and Kantharia32). Because our experimental design was a 4 × 4 Latin square, we then settled on using three animals per treatment in each period which cumulatively gave us a total sample size of twelve at the end of the four periods. We used change in digestibility as our reference response variable instead of MPR because there were no experimental literatures exploring CH4 production at intakes below MER.
Results
The basal Rhodes grass diet varied in organic matter between periods, albeit over a narrow range (911–929 (sem 0·32) g/kg DM; P < 0·05), but not in CP (P = 0·347). However, refusals, were lower in CP compared with rations fed (Table 2), indicating that animals were selecting for higher quality parts of the hay, but animals receiving the MER 0·4 and MER 0·6 treatments consistently consumed all of their allocated ration and had no refusals.
Intakes differed between each treatment group as intended (Table 3), but intakes were lower in period 2 (P < 0·05) which was attributable to decreased intakes of the animals in 0·8 MER group, relative to the first period (P = 0·06). Animals in the 1·0 MER treatment group showed net weight gain, indicating that, even after adjusting the initial ration calculations, we overestimated MER, but the other three treatment groups lost weight in a linear manner (Table 3).
a,b,c,d Mean values within a row with unlike superscript letters were significantly different (P < 0·05).
* Calculated over the full 35-d feeding period.
Intakes of the 1·0 and 0·8 MER groups were less during the 6 d collection and during the MPR measurement periods (P = 0·05) than intakes during the 21 d measurement period, but these differences were small (less than 2 %) and were not further considered (Table 4). However, there were no differences in faecal composition or apparent digestibility between any of the treatment groups, although there was a trend for MER 0·4 to have lower apparent digestibility of organic matter (Table 4).
a,b,c,d Mean values within a row with unlike superscript letters were significantly different (P < 0·05).
* Displayed a non-significant trend.
Recoveries of CH4 for individual chambers ranged from a mean of 96·3 to 102·7 % across the four measurement periods. Daily MPR differed between all treatment groups and showed a clear linear relationship to intake as anticipated (MPR = 109·8 × MER + 25·6; P < 0·001). Interestingly, MY increased modestly, also in a linear fashion as intakes and MPR fell further below 1·0 MER (MY = −5·23 × MER + 34·3; P = 0·001), with an increase of about 8 % (P < 0·001) in the 0·4 MER group compared with 1·0 MER (Table 5). CH4 produced per kg of digested organic matter (−10·25 × MER + 58·7; P < 0·001) and Y m (Y m = −1·81 × MER + 11·03) followed a similar trend, with increases of approximately 10 % in the 0·4 MER group compared with 1·0 MER (Table 5). In contrast, rumen MRT was prolonged in both liquid and solid phases, also in a linear fashion (MRTliquid (h) = −11·6 × MER + 31·8; P < 0·001; MRTsolid (h) = −30 × MER + 92·0; P < 0·001) in all sub-maintenance groups by similar amounts of time, compared with the 1·0 MER treatment (Table 5).
GEI, gross energy intake.
a,b,c,d Mean values within a row with unlike superscript letters were significantly different (P < 0·05).
Discussion
The overarching aim of this study was to determine if cattle consuming rations at levels well below MER for maintenance produce more CH4 than would be expected from the amount of feed ingested. We initially estimated MER from the National Research Council(21) recommendation for growing B. taurus cattle, subsequently decreasing allocations by 20 %, based on our empirical observations. However, the (nominally) 1·0 MER group still gained weight during the trial, and from (linear regression) analysis of LW records (eq. 2), we determined that our revised ration still overestimated MER by about 10 %. This is consistent with findings of a recent meta-analysis(Reference Oliveira33) which suggests that the MER of B. indicus × B. taurus was up to 26 % less than that of pure B. taurus cattle. Ultimately, the miscalculation does not detract from the validity of the study, as it is clear from LW change that animals in all other treatment groups were in energy deficit and thus feed intake was below MER. Our thesis that both MY and Y m would increase when cattle consume feed well-below MER was based partly on the idea that MRT is inversely related to feed intake, although differences in feed intake featured in the literature are most often related to the neutral-detergent fibre content of forage consumed ad libitum (Reference Van Amburgh, Grant and Cotanch34). Prolonged residence time was observed in this study in both liquid and solid digesta fractions, but this was due to the restriction in intake rather than the limitation imposed by physical fill. It is perhaps important that the largest increase in MRT occurred between the 1·0 and 0·8 MER groups, with little prolongation seen past this. The larger difference in MRT between 1·0 MER and the other groups might be associated with 1·0 MER receiving cotton seed meal and molasses as part of their ration, but this is unproven, and even if it is the case, does not detract from the major findings of this study because the larger and significant changes to MY and Y m occurred at the 0·4 and 0·6 levels. Thus, both MY and Y m continued to increase as intake levels fell, but the changes were not accompanied by similar changes in MRT, nor by increases in apparent digestibility of the feed eaten. The observed results align with a study of Doreau et al. (Reference Doreau, Michalet-Doreau and Béchet31), which found digestibility actually decreased at extremely low levels of intake and tellingly was not improved by additional rumen-degradable protein, indicating that reduced digestibility was probably not attributable to reduced rumen microbial activity. Instead, the amount of CH4 produced per unit of digested feed increased as feeding levels fell, indicating that a greater proportion of feed actually assimilated was being ultimately diverted to the production of CH4 – as suggested by the changes to both MY and Y m. The shift towards increased CH4 is consistent with a shift in volatile fatty acid production (towards acetate) and/or reduced production or refermentation of microbial protein, but this study does not furnish direct evidence of this.
Our findings differ from Y m values currently in use, but other evidence exists to support their likely validity. We found that the Y m of cattle fed at approximately maintenance (9·1 %) is higher than both Intergovernmental Panel on Climate Change recommendations (6·5 %) and recently published estimates(Reference Charmley, Williams and Moate8) of 6·3 % and our estimates of Y m for cattle fed at 40 % of MER was greater still, at 10·1 %. However, another trial feeding tropical forages under ad libitum conditions reported a greater range of Y m of 5·0–7·2 %(Reference Kennedy and Charmley35), with the highest values occurring when cattle were fed very low-quality grass, at what would have been effectively sub-maintenance requirements. Kaewpila & Sommart(Reference Kaewpila and Sommart36) calculated (from earlier studies) a Y m of 8·2 % for B. indicus cattle fed low-quality forage – but a limitation of this study is that it does not seek to separate the possible influences of feeding levels from feed quality. While each of the studies detailed above, as well as the one reported here, has been conducted using cattle of predominantly B. indicus origins, it could be tempting to attribute the higher Y m values to breed effects. However, we assert that we have clearly demonstrated the contribution of feeding levels apart from any genetic effects that might exist. To clarify this, it would be desirable to conduct further measurements of MPR in cattle consuming tropical forages ad libitum to determine the effect on Y m values.
Conclusion
Feeding low-quality tropical forages at restricted intakes substantially increase both the MY and Y m of the consumed forages, whilst severe intake restriction exacerbates the case by a further 8–10 %. While meriting further investigation to confirm, we believe these changes observed in animals fed well below their voluntary intake are quantitatively important for both greenhouse gas inventory and animal husbandry in tropical rain-fed systems. Firstly for ruminants in areas where there are marked seasonal nutritional deficits such as in Sub-Saharan Africa (and northern Australia), incorporation of a higher Y m value (as suggested here) for part of the year would increase calculated emission factors by 10 % or more, which may significantly alter current estimates of CH4 emissions from livestock, for example, in Africa(Reference Tubiello, Salvatore and Cóndor Golec9). More importantly, perhaps, is that this work suggests that improving the feeding of ruminants under nutritional stress has not only the potential to reduce their contribution to greenhouse gas emissions but also direct economic benefits for farmers via improved nutrient partitioning away from the production of CH4 towards growth.
Acknowledgements
We gratefully acknowledge the assistance of Mr. S. Mwangi, Mr. A. Mbithi, Mr. S. Gwatibaya, Mr. R. Rogito, Mr. G. N. Wanyama and Mr. V. Omondi during the animal experiment at International Livestock Research Institute’s Mazingira Centre (mazingira.ilri.org).
This study was funded by the German Federal Ministry for Economic Cooperation and Development within the research project ‘In situ assessment of greenhouse gas emissions from livestock production systems in East Africa’ (grant no. 55 21914). We acknowledge the CGIAR Fund Council, Australia (ACIAR), Irish Aid, European Union, International Fund for Agricultural Development (IFAD), the Netherlands, New Zealand, UK, USAID and Thailand for funding to the CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS) and the CGIAR Research Program on Livestock.
No funder (as stipulated above) had any role in the design, analysis or writing of this article.
J. P. G.: original intellectual concept and study design; assembly and analysis of data; primary responsibility for writing of article. D. K.: leading role in conduct of trial, acquisition and assembly of data; ,ajor role in data analysis and production of results; involved in writing and revisions of article. D. E. P.: critical review of study design, detailed statistical analysis, critical review of data results and significant contribution to discussion. A. I. M. A.: involvement in conduct of trial; responsible for acquisition and calculation of data on digesta passage parameters. S. E. W.: involved in conduct of trial; responsible for acquisition of feed intake data and revisions of the article. E. S.: supporting role in acquisition and analysis of passage rate data, critical revisions of the article. U. D. (U. D. has not responded to a request to specify her contribution). L. M.: supporting role in conduct of trial and acquisition of CH4 data, critical review of results and involved in writing and revisions of the article. K. B.-B.: critical review of study design, leading role in experimental setup of respiration chambers and data analysis, critical review of data and results, significant contribution to discussion and writing.
The authors declare that they have no conflicts of interest.
Supplementary material
For supplementary material referred to in this article, please visit https://doi.org/10.1017/S0007114519003350