Research | Open | Published:
Dynamics of methanogenesis, ruminal fermentation and fiber digestibility in ruminants following elimination of protozoa: a meta-analysis
Journal of Animal Science and Biotechnologyvolume 9, Article number: 89 (2018)
Ruminal microbes are vital to the conversion of lignocellulose-rich plant materials into nutrients for ruminants. Although protozoa play a key role in linking ruminal microbial networks, the contribution of protozoa to rumen fermentation remains controversial; therefore, this meta-analysis was conducted to quantitatively summarize the temporal dynamics of methanogenesis, ruminal volatile fatty acid (VFA) profiles and dietary fiber digestibility in ruminants following the elimination of protozoa (also termed defaunation). A total of 49 studies from 22 publications were evaluated.
The results revealed that defaunation reduced methane production and shifted ruminal VFA profiles to consist of more propionate and less acetate and butyrate, but with a reduced total VFA concentration and decreased dietary fiber digestibility. However, these effects were diminished linearly, at different rates, with time during the first few weeks after defaunation, and eventually reached relative stability. The acetate to propionate ratio and methane production were increased at 7 and 11 wk after defaunation, respectively.
Elimination of protozoa initially shifted the rumen fermentation toward the production of more propionate and less methane, but eventually toward the production of less propionate and more methane over time.
The rumen provides an ideal habitat for protozoa, whose concentration can reach 105–106 cells/mL. In return, protozoa serve important functions in the rumen microbial ecosystem, such as predation, competition for nutrients, and involvement in symbiotic relationships with other microorganisms [1, 2]. Protozoa prey on bacteria and fungal spores, but are preferentially retained in the rumen, thus reducing the postruminal microbial protein supply . Protozoa compete with amylolytic bacteria for dietary starch, which is mostly fermented into acetate by protozoa  while mostly into propionate by amylolytic bacteria [3, 4]. For host animals, the energy recovery efficiency is reduced by 38% when the substrate (glucose) is fermented into acetate but increased by 9% when fermented into propionate [3, 4]. Protozoa are important ruminal hydrogen (H2) producers, and the produced H2 is mostly converted into methane (CH4) by methanogens situated inside protozoa or on their external surface [5,6,7]. The CH4 emissions from ruminants represent 2–12% dietary energy loss . Therefore, the presence of protozoa seems to adversely affect animals’ energy efficiency.
Complete removal of ruminal protozoa, termed defaunation, has been suggested as an efficient method for reducing CH4 emissions and enhancing propionate fermentation [9, 10], but these effects have not been consistently observed in studies investigating this 112method [11,12,13]. Hegarty et al.  and Morgavi et al.  suggested that the duration of defaunation might be responsible for this inconsistency, but the temporal dynamics of methanogenesis and ruminal volatile fatty acid (VFA) profiles after defaunation are difficult to determine experimentally because of the difficulties in raising defaunated animals over a long-term . Meta-analysis is a statistical method that combines the results from multiple studies to achieve a more precise estimate of treatment effects and to explore the potential sources of between-study heterogeneity [15, 16]. Two prior meta-analyses [14, 17] have summarized the combined responses of rumen fermentation to defaunation; however, the combined effects on ruminal VFA profiles were also inconsistent between them, and neither of them explored the time-dependent effects. Therefore, the current meta-analysis was conducted to quantitatively summarize the temporal dynamics of methanogenesis, ruminal VFA profiles and dietary fiber digestibility in ruminants after defaunation, and to explain the contribution of the defaunated duration to the between-study heterogeneity.
Literature search, screening and data extraction
A flowchart detailing the process of literature search, screening and data extraction is shown in Fig. 1. The inclusion criteria for the studies were as follows: (1) peer-reviewed and published in the English language; (2) complete defaunation in vivo; (3) inclusion of relevant variables for extraction. The relevant variables for this meta-analysis included the daily CH4 production, ruminal VFA profiles and dietary fiber total-tract digestibility.
The faunation states of control animals included the natural ruminal ecosystem without any treatment or reintroduction of protozoa after partial defaunation, which appeared to be restored quickly after withdrawal of the protozoa-inhibiting treatment [18, 19]. The control animals into which protozoa were reintroduced after complete defaunation were excluded from the analysis, because preliminary analysis showed that high heterogeneity existed between the faunation and refaunation subgroups (see Additional file 1). Protozoa-free animals were obtained through either the absence of protozoa from birth (BF) or artificial removal of protozoa from the natural ruminal ecosystem (AF). Artificial defaunation was conducted using chemical agents, such as alkanes and sodium lauryl sulfate, or by applying a rumen washing technique. The defaunation duration was calculated based on the schedule of experimental activities. The defaunation duration of BF animals was calculated as their age shortened by 4 wk, because ruminal protozoa did not appear when newborn calves were fed with milk for 30 d , and the concentration of protozoa quickly increased after 5 wk of age .
The final database included 22 publications, with 49 in vivo studies that satisfied the inclusion criteria for the meta-analysis. Summary descriptions of the selected studies are provided in Table 1. Briefly, among the 49 selected studies, 38 were conducted in sheep, and 11 were conducted in cattle; 18 measured CH4 production, 29 measured ruminal VFAs, and 15 measured total-tract fiber digestibility. CH4 emissions were measured using the respiration chamber technique or sulfur hexafluoride tracer technique. Methane production was presented in liters per day in most of the studies; thus, values presented in grams or kilojoules (kJ) were converted to liters per day, based on the assumption that one mole of CH4 weighs 16 g or contains 890 kJ of energy and occupies a volume of 24.5 L (under conditions of 25 °C and 1 atmospheric pressure).
The meta-analysis was performed using Stata 14.1 (Stata Corp., Texas, USA).
Assessment of heterogeneity and effect size
Between-study variability was quantified via the I2 statistic, which measures the percentage of variation due to heterogeneity [15, 21]. When the I2 value was less than 50%, indicating low heterogeneity, studies were combined using a fixed effects model, which was based on the assumption that the expected effect from each study was homogeneous. When the value of I2 was over 50%, indicating high heterogeneity, studies were combined using a random effects model, based on the assumption that the expected effect from each study was heterogeneous.
The differences in animal species or ages, daily sampling times or dietary forage percentages across the studies caused that the data of certain relevant variables to vary greatly across the studies (Table 2). To reduce these potential interferences, the effect size in this analysis was estimated via the standardized mean difference (SMD), which was calculated as the raw mean difference between the treatment and control groups divided by their pooled standard deviations . For example, although cattle CH4 production in the study by Schönhusen et al.  was higher than in other subgroup studies involving sheep, it was homogeneous with most of them (Additional file 1: Figures S1 and S2). The studies were weighted using the inverse of the variance of the differences in means. Details of the calculations used in the meta-analysis are provided by Lean et al. .
The meta-regression analysis was performed using the Knapp-Hartung restricted maximum likelihood method , with the SMD of the individual studies used as the response variable and the corresponding standard error of the SMD used as the variance. The percentage of between-study heterogeneity explained by the covariate (defaunation duration) was quantified via the adjusted R2 value.
Preliminary analysis showed that the temporal SMD dynamics after defaunation included a linear phase followed by a plateau phase. The durations of the linear phase and plateau phase for each outcome were dependent on the highest adjusted R2 of the meta-regression analysis and the minimum I2 of the heterogeneity analysis, respectively; therefore, the two phases might overlap over a short duration. The probability levels were set at P < 0.05 for significance and 0.05 ≤ P < 0.10 for a trend.
Effect size and heterogeneity across all the studies
The meta-analysis based on all the selected studies showed that elimination of rumen protozoa reduced (P < 0.05) CH4 production, ruminal VFA concentration, the proportion of butyrate and dietary fiber digestibility, and tended to increase (P = 0.083) the proportion of acetate (Table 2). However, the heterogeneity across the studies was considerable (I2 > 50%) for most of the responses to defaunation, except for the ruminal VFA concentration (I2 = 44.2%).
Methanogenesis dynamics during adaptation to defaunation
Compared with that of faunation, the effect size of defaunation on CH4 emissions presented a linear relationship over time during the first 12 wk (linear phase) after defaunation (Fig. 2 and Table 3): CH4 production was reduced by defaunation (intercept = − 5.484, P = 0.003), and the reduction decreased weekly by 0.486 (P = 0.003) until 12 wk. The defaunation duration explained 76.8% of the between-study heterogeneity (I2 = 76.2%) during the linear phase, and no between-study heterogeneity (I2 = 0.0%) was observed after 11 wk (plateau phase), suggesting that little fluctuation occurred during the plateau phase. Interestingly, the defaunated animals during the plateau phase presented higher CH4 productions (SMD = 0.313, P = 0.039) than the control animals.
Ruminal VFA profiles and total-tract fiber digestibility dynamics during adaptation to defaunation
Consistent with the temporal dynamics of CH4 production after defaunation, the ruminal VFA profiles dynamics also included a linear phase (≤ 11 wk) and a plateau phase (≥ 7 wk) (Table 3). After defaunation, decreases in the acetate proportion, butyrate proportion and A:P (intercept = − 4.086, − 7.059 and − 6.737, respectively), and an increase in the propionate proportion (intercept = 7.306) were estimated (P < 0.01). These alterations decreased linearly (P < 0.01) at different rates over the first 11 wk of defaunation; instead, the ruminal acetate proportion (SMD = 0.748, P < 0.001) and A:P (SMD = 0.915, P = 0.016) were higher, and the propionate proportion (SMD = − 0.366, P = 0.033) was lower in defaunated animals than faunated animals during the plateau phase.
Compared with faunated animals, defaunated animals exhibited a reduced total VFA concentration (intercept = − 1.883 and P = 0.008), and the reduction decreased weekly by 0.132 (P = 0.048) until 12 wk (Table 3). The duration of defaunation could explain 99.9% of the between-study heterogeneity during the linear phase, although the heterogeneity (35.5%) was low. The decrease in total VFA concentration was still observed (SMD = − 0.424, P < 0.001) in defaunated animals during the plateau phase (≥ 11 wk).
Compared with that of faunation, the effect sizes of defaunation on total-tract fiber digestibility were linearly related to defaunation duration (Table 3): total-tract NDF and ADF digestibility (intercept = − 4.458 and − 6.276, respectively) were reduced (P < 0.01) after defaunation, and the reductions decreased weekly (P < 0.05) by 0.153 and 0.213, respectively, until 23 wk, which was the longest studied duration in the available data.
Effect size and heterogeneity across all the studies
Based on all the selected studies, this meta-analysis showed that complete elimination of rumen protozoa generated adverse effects on the ruminal VFA concentration, butyrate proportion and dietary fiber digestibility; these findings were consistent with the results of previous meta-analyses [14, 18]. However, the heterogeneity across the studies was considerable for most of the responses to defaunation. Excess between-study variance increases the risk of incorrect average effect sizes when combining studies . For example, the present meta-analysis based on all the studies showed that defaunation tended to increase the proportion of ruminal acetate but had no effect on the proportion of propionate. These findings were consistent with a recent meta-analysis by Newbold et al.  but inconsistent with that of Eugène et al. , who reported that defaunation induced a reduction in the ruminal acetate proportion and an increase in the propionate proportion. Therefore, the potential source of heterogeneity among the studies needs to be explored to better understand the responses to treatment, and this additional exploration is also one of most important tasks of meta-analysis .
The role of protozoa in rumen carbohydrate metabolism
Despite the fact that protozoa make up a large portion of the rumen biomass, their role in ruminal fermentation and their contribution to the metabolism and nutrition of the host are still topics of substantial controversy [2, 14], due to the difficulty of pure cultivation of protozoa in vitro. Rumen protozoa are not essential to the animal for survival, and defaunation has therefore been used to estimate the role of ciliate protozoa in rumen function. However, the adapted alteration of other microbe after defaunation may interfere with such estimations. Hence, the estimated intercept from the meta-regression analysis more accurately reflected the role that protozoa played in rumen fermentation. Reductions in the ruminal acetate proportion, butyrate proportion and CH4 production after complete removal of ruminal protozoa would be expected because protozoa ferment carbohydrates into acetate, butyrate and H2 , and the H2 produced is mostly converted to CH4 by methanogens situated inside protozoa or on their external surface [5,6,7]. A strong correlation between CH4 emissions and protozoa concentration has been reported , and protozoa-associated methanogens have been estimated to be responsible for 37% of CH4 production by ruminants . Additionally, ruminal protozoa possess a full complement of hydrolytic enzymes for fermentation of the major components of feedstuffs, and certain ciliates present a wide range of fibrolytic enzyme genes, ingest small plant particles and use cell wall carbohydrates [18, 24]. Moreover, protozoa can indirectly contribute to ruminal degradation kinetics by maintaining a suitable rumen fermentation environment, for example, by scavenging oxygen to maintain anaerobiosis and slowing the rate of starch fermentation to maintain a proper ruminal pH , which favors the development and activity of bacteria and fungi [25, 26]. Therefore, the reductions in dietary fiber digestibility and ruminal total VFA concentration observed in this study would be expected after the complete removal of ruminal protozoa.
Temporal dynamics during adaptation to defaunation
Although protozoa are important ruminal H2 producers and exhibit interspecies H2 transfer with methanogens, we found that the effects of short- and long-term defaunation on CH4 production were opposite. This finding supports the conclusion of Morgavi et al. , who showed that there was not a simple cause-effect relationship between rumen protozoa and methanogenesis. Bird et al.  and Hegarty et al.  observed higher CH4 production (although not significantly so) in long-term defaunated ewes (11 and 26 wk) and lambs (12 to 33 wk) than in faunated animals. The significant increasing effect of long-term defaunation on CH4 production detected in this meta-analysis can be attributed to the pooled analysis, in which the number of replicate animals was increased by combining the results of relevant individual studies [15, 16]. The CH4 emissions from ruminants contribute to global greenhouse gas emissions and represent energy loss for the animals [8, 27, 28]. Therefore, the potential environmental protection and energy-saving values following defaunation were gradually lost and eventually became negative.
Acetate production during rumen fermentation is accompanied by reducing equivalents ([H]) production, whereas propionate production is accompanied by [H] consumption ; the excess [H] is converted to H2. The shift in the VFA profiles from propionate to acetate following defaunation increased the H2 available for methanogenesis, at least partially explaining the time-related changes in CH4 production observed in this study. When glucose is metabolized into acetate, propionate or butyrate, the energy efficiency relative to glucose for animal is 62%, 109% and 78%, respectively [3, 4]. Propionate fermentation is most energy efficient, due to assimilating energy from H2 and being the main precursor of gluconeogenesis in animals [3, 30]. In ruminants, the VFAs produced in the rumen satisfy up to 70% of energy requirements . Shabat et al.  and Weimer et al.  observed that the ruminal total VFA concentration and propionate proportion were higher in highly efficient cows than in cows with low efficiency. Therefore, the decreases in the ruminal total VFA concentration and propionate proportion during the plateau phase also suggested that the elimination of rumen protozoa adversely affected the energy supply of animals in a long run.
The time-related variations in CH4, VFA profiles and dietary fiber digestibility implied a series of complex changes in the ruminal ecosystem over the course of defaunation. Nguyen et al.  reported that rumen microbes had likely not stabilized after 12 wk of defaunation, which agrees with our results showing that the linear phase for ruminal VFA profiles and CH4 emissions lasted 11 wk and 12 wk, respectively. When sudden major changes are made in the diet, it takes approximately 2 wk for the new microbial population balance to become established ; the much longer linear phase associated with defaunation suggests that protozoa play an important role in the ruminal ecosystem. Ruminal methanogens appear to develop more slowly than bacteria following defaunation . Hristov et al.  noted that reductions in the population of protozoa-associated methanogens might be compensated by an increase in the population of bacteria- or rumen fluid-associated methanogens, and Mosoni et al.  found that long-term defaunation (2 yr) increased the abundance of methanogens. In addition, ruminal protozoa elimination results in increased bacterial abundance and changes in bacterial communities [35, 36]; defaunation has been shown to increase the anaerobic fungal population by two fold  and the Ruminococcaceae population by six fold . Frumholtz  found that long-term defaunation (6 mo) increased the abundance of cellulolytic bacteria. Similar to protozoa, fungi and cellulolytic bacteria are also the main ruminal cellulolytic and H2-producing microbes that generate acetate, butyrate and/or H2 as primary end products [29, 39]. Therefore, it can be concluded that the increase in the populations of methanogens, fungi and cellulolytic bacteria following defaunation gradually counteracts the defaunation-induced reductions in dietary fiber digestibility, ruminal A:P and CH4 production, which may confirm an earlier theory of Weimer  indicating that the multiple microbial taxa in the ruminal community show functional redundancy (overlap of physiological function) and may therefore be substitutable with little impact on ecosystem processes [41, 42]. As noted by Taxis et al.  regarding the relationship between ruminal ecosystems and function: the players may change but the game remains. These observations also suggest that defaunation is not a good model for estimating the role of protozoa in rumen function due to the compensation effects of fungi and bacteria. Further animal experiments are required to fully understand the succession of rumen bacterial and archaeal community structure and function following defaunation, and the metabolic characteristics of rumen protozoa need be revealed using their genome and transcriptome data.
The present meta-analysis summarized the temporal dynamics of methanogenesis, ruminal fermentation and dietary fiber digestibility in ruminants after defaunation, and the results showed that defaunation adversely affected dietary fiber digestibility and the ruminal VFAs available to the host animals, although the effects were lessened over time. Furthermore, the energy advantages of defaunation gained by reducing CH4 production and shifting ruminal VFA profiles to more propionate were gradually lost over time, and the effects eventually became disadvantageous. Therefore, elimination of rumen protozoa adversely affects the energy supply of animals over the long-term.
Acid detergent fiber
Acetate to propionate ratio
Neutral detergent fiber
Standardized mean difference
Volatile fatty acids
Ushida K, Jouany JP, Demeyer DI. Effects of presence or absence of rumen protozoa on the efficiency of utilization of concentrate and fibrous feeds. In: Tsuda T, Sasaki Y, Kawashima, editors. Physiological aspects of digestion and metabolism in ruminants. San Diego: Academic Press. 1991:625–54.
Williams AG, Coleman GS. The rumen protozoa. 1st ed. New York: Springer-Verlag; 1992.
Millen DD, Arrigoni MDB, Pacheco RDL. Rumenology. 1st ed. Springer: International Publishing; 2016.
Ryle M, Ørskov ER. Energy nutrition in ruminants. 1st ed. Netherlands: Springer; 1990.
Finlay BJ, Esteban G, Clarke KJ, Williams AG, Embley TM, Hirt RP. Some rumen ciliates have endosymbiotic methanogens. FEMS Microbiol Lett. 1994;117:157–61.
Morgavi DP, Forano E, Martin C, Newbold CJ. Microbial ecosystem and methanogenesis in ruminants. Animal. 2010;4:1024–36.
Belanche A, de la Fuente G, Newbold CJ. Study of methanogen communities associated with different rumen protozoal populations. FEMS Microbiol Ecol. 2014;90:663–77.
Johnson KA, Johnson DE. Methane emissions from cattle. J Anim Sci. 1995;73:2483–92.
Whitelaw FG, Eadie JM, Bruce LA, Shand WJ. Methane formation in faunated and ciliate-free cattle and its relationship with rumen volatile fatty acid proportions. Br J Nutr. 1984;52:261–75.
Faichney GJ, Graham NM, Walker DM. Rumen characteristics, methane emissions, and digestion in weaned lambs reared in isolation. Aust J Agric Res. 1999;50:1083–90.
Bird SH, Hegarty RS, Woodgate R. Persistence of defaunation effects on digestion and methane production in ewes. Aust J Exp Agric. 2008;48:152–5.
Hegarty RS, Bird SH, Vanselow BA, Woodgate R. Effects of the absence of protozoa from birth or from weaning on the growth and methane production of lambs. Br J Nutr. 2008;100:1220–7.
Morgavi DP, Martin C, Jouany JP, Ranilla MJ. Rumen protozoa and methanogenesis: not a simple cause–effect relationship. Br J Nutr. 2012;107:388–97.
Eugène M, Archimède H, Sauvant D. Quantitative meta-analysis on the effects of defaunation of the rumen on growth, intake and digestion in ruminants. Livest Prod Sci. 2004;85:81–97.
Lean IJ, Rabiee AR, Duffield TF, Dohoo IR. Invited review: use of meta-analysis in animal health and reproduction: methods and applications. J Dairy Sci. 2009;92:3545–65.
Viechtbauer W. Learning from the past: refining the way we study treatments. J Clin Epidemiol. 2010;63:980–2.
Newbold CJ, de la Fuente G, Belanche A, Ramos-Morales E, McEwan NR. The role of ciliate protozoa in the rumen. Front Microbiol. 2015;6:1313.
Sauer FD, Fellner V, Kinsman R, Kramer JK, Jackson HA, Lee AJ, et al. Methane output and lactation response in Holstein cattle with monensin or unsaturated fat added to the diet. J Anim Sci. 1998;76:906–14.
Guan H, Wittenberg KM, Ominski KH, Krause DO. Efficacy of ionophores in cattle diets for mitigation of enteric methane. J Anim Sci. 2006;84:1896–906.
Schönhusen U, Zitnan R, Kuhla S, Jentsch W, Derno M, Voigt J. Effects of protozoa on methane production in rumen and hindgut of calves around time of weaning. Arch Anim Nutr. 2003;57:279–95.
Higgins JPT, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.
Knapp G, Hartung J. Improved tests for a random effects meta-regression with a single covariate. Stat Med. 2003;22:2693–710.
Guyader J, Eugène M, Nozière P, Morgavi DP, Doreau M, Martin C. Influence of rumen protozoa on methane emission in ruminants: a meta-analysis approach. Animal. 2014;8:1816–25.
Nagaraja TG. Microbiology of the rumen. In: Millen DD, Arrigoni MDB, Pacheco RDL, editors. Rumenology. Cham: Springer; 2016. p. 39–61.
Jouany JP, Demeyer DI, Grain J. Effect of defaunating the rumen. Anim Feed Sci Technol. 1988;21:229–65.
Mosoni P, Martin C, Forano E, Morgavi DP. Long-term defaunation increases the abundance of cellulolytic ruminococci and methanogens but does not affect the bacterial and methanogen diversity in the rumen of sheep. J Anim Sci. 2011;89:783–91.
Karakurt I, Aydin G, Aydiner K. Sources and mitigation of methane emissions by sectors: a critical review. Renew Energy. 2012;39:40–8.
Hristov AN, Oh J, Lee C, Meinen R. Mitigation of greenhouse gas emissions in livestock production: a review of technical options for non-CO2 emissions. Food and Agriculture Organization of the United Nations: Rome; 2013.
Janssen PH. Influence of hydrogen on rumen methane formation and fermentation balances through microbial growth kinetics and fermentation thermodynamics. Anim Feed Sci Technol. 2010;160:1–22.
Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol Rev. 1990;70:567–90.
Shabat SKB, Sasson G, Doron-Faigenboim A, Durman T, Yaacoby S, Berg Miller ME, et al. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 2016. https://doi.org/10.1038/ismej.2016.62.
Weimer PJ, Cox MS, de Paula TV, Lin M, Hall MB, Suen G. Transient changes in milk production efficiency and bacterial community composition resulting from near-total exchange of ruminal contents between high- and low-efficiency Holstein cows. J Dairy Sci. 2017;100:7165–82.
Nguyen SH, Bremner G, Cameron M, Hegarty RS. Methane emissions, ruminal characteristics and nitrogen utilisation changes after refaunation of protozoa-free sheep. Small Rumin Res. 2016;144:48–55.
Reece WO. Dukes’ physiology of domestic animals. 12th ed. Comstock Pub. Associates: Ithaca; 2004.
Ozutsumi Y, Tajima K, Takenaka A, Itabashi H. The effect of protozoa on the composition of rumen bacteria in cattle using 16S rRNA gene clone libraries. Biosci Biotechnol Biochem. 2005;69:499–506.
Morgavi DP, Rathahao-Paris E, Popova M, Boccard J, Nielsen KF, Boudra H. Rumen microbial communities influence metabolic phenotypes in lambs. Front Microbiol. 2015;6:1060.
Newbold CJ, Hillman K. The effect of ciliate protozoa on the turnover of bacterial and fungal protien in the rumen of sheep. Lett Appl Microbiol. 1990;11:100–2.
Frumholtz PP. Manipulation of the rumen fermentation and its effects on digestive physiology: University of Aberdeen; 1991.
Morvan B, Rieu-Lesme F, Fonty G, Gouet P. In vitro interactions between rumen H2-producing cellulolytic microorganisms and H2-utilizing acetogenic and sulfate-reducing bacteria. Anaerobe. 1996;2:175–80.
Weimer PJ. Redundancy, resilience, and host specificity of the ruminal microbiota: implications for engineering improved ruminal fermentations. Front Microbiol. 2015;6:296.
Lawton JH, Brown VK. Redundancy in ecosystems. In: Schulze ED, Mooney HA, editors. Biodiversity and ecosystem function. Berlin: Springer; 1994. p. 255–70.
Rosenfeld JS. Functional redundancy in ecology and conservation. Oikos. 2002;98:156–62.
Taxis TM, Wolff S, Gregg SJ, Minton NO, Zhang C, Dai J, et al. The players may change but the game remains: network analyses of ruminal microbiomes suggest taxonomic differences mask functional similarity. Nucleic Acids Res. 2015. https://doi.org/10.1093/nar/gkv973.
Belanche A, Abecia L, Holtrop G, Guada JA, Castrillo C, de la Fuente G, et al. Study of the effect of presence or absence of protozoa on rumen fermentation and microbial protein contribution to the chyme. J Anim Sci. 2011;89:4163–74.
Chandramoni JSB, Tiwari CM, Haque N, Murarilal KMY. Energy metabolism and methane production in faunated and defaunated sheep fed two diets with same concentrate to roughage ratio (70:30) but varying in composition. Asian-Australas J Anim Sci. 2001;14:1238–44.
Chaudhary LC, Srivastava A. Performance of growing Murrah buffalo calves as affected by treatment with Manoxol and the presence of ciliate protozoa in the rumen. Anim Feed Sci Technol. 1995;51:281–6.
Eadie JM, Gill JC. The effect of the absence of rumen ciliate protozoa on growing lambs fed on a roughage–concentrate diet. Br J Nutr. 1971;26:155–67.
Eugène M, Sauvant D, Weisbecker JL, Archimède H. Effects of defaunation on digestion of fresh Digitaria decumbens grass and growth of lambs. Animal. 2010;4:439–45.
Kasuya N, Wada I, Shimada M, Kawai H, Itabashi H. Effect of presence of rumen protozoa on degradation of cell wall constituents in gastrointestinal tract of cattle. Anim Sci J. 2007;78:275–80.
Kreuzer M, Kirchgessner M, Müller HL. Effect of defaunation on the loss of energy in wethers fed different quantities of cellulose and normal or steamflaked maize starch. Anim Feed Sci Technol. 1986;16:233–41.
Nagaraja TG, Godfrey SI, Winslow SW, Rowe JB, Kemp KE. Effect of virginiamycin on ruminal fermentation in faunated or ciliate-free sheep overfed with barley grain. Small Rumin Res. 1995;17:1–8.
Santra A, Karim SA. Growth performance of faunated and defaunated Malpura weaner lambs. Anim Feed Sci Technol. 2000;86:251–60.
Santra A, Karim SA. Nutrient utilization and growth performance of defaunated and faunated lambs maintained on complete diets containing varying proportion of roughage and concentrate. Anim Feed Sci Technol. 2002;101:87–99.
Santra A, Karim SA, Chaturvedi OH. Rumen enzyme profile and fermentation characteristics in sheep as affected by treatment with sodium lauryl sulfate as defaunating agent and presence of ciliate protozoa. Small Rumin Res. 2007;67:126–37.
Sultana H, Miyazawa K, Kanda S, Itabashi H. Fatty acid composition of ruminal bacteria and protozoa, and effect of defaunation on fatty acid profile in the rumen with special reference to conjugated linoleic acid in cattle. Anim Sci J. 2011;82:434–40.
Williams PP, Dinusson WE. Ruminal volatile fatty acid concentrations and weight gains of calves reared with and without ruminal ciliated protozoa. J Anim Sci. 1973;36:588–91.
Yáñez-Ruiz DR, Williams S, Newbold CJ. The effect of absence of protozoa on rumen biohydrogenation and the fatty acid composition of lamb muscle. Br J Nutr. 2007;97:938–48.
Zhou YY, Mao HL, Jiang F, Wang JK, Liu JX, McSweeney CS. Inhibition of rumen methanogenesis by tea saponins with reference to fermentation pattern and microbial communities in Hu sheep. Anim Feed Sci Technol. 2011;166–167:93–100.
The authors thank Prof. Chang Xu (Chinese Evidence-Based Medicine Center and Chinese Cochrane Center, West China Hospital, Sichuan University) for help with data analysis.
This work was supported by the National Key Research and Development Program of China (grant number: 2017YFD0500500).
Availability of data and materials
All the datasets were presented in the main manuscript (reference list in Table 1) and available to readers.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Figure S1. Forest plot showing the results of the subgroup meta-analysis of the anti-methanogenic effect size of defaunation, grouped by faunation state and duration of defaunation (11 wk). BF = born and reared protozoa free; AF = artificial defaunation; SMD = standardized mean difference; 95% CI = 95% confidence interval. * I-squared = percentage of heterogeneity across studies; P-value of SMD = 0. Figure S2. Funnel plot for the effect size of defaunation on CH4 production in (A) all studies, (B) short-term defaunation, (C) long-term defaunation, and (D) refaunation. The P-value of publication bias is presented. SMD = standardized mean difference, se = standard error. (DOCX 1373 kb)