Skip to main content

High throughput accurate method for estimating in vitro dry matter digestibility of maize silage



Since the introduction of studies on maize silage digestibility at the end of the nineteenth century, protocols to estimate dry matter digestibility have not stopped evolving. Since the early 1980s, the protocol developed by Aufrère became a benchmark in many laboratories to estimate in vitro dry matter digestibility. In order to increase its throughput, to facilitate its execution and to decipher the impact of the different parameters of the protocol we decided to test the combination of 7 parameters in 21 different protocols.


We thus tested the impact of (1) the presence or absence of pepsin in HCl solution, (2) the temperature of incubation during enzymatic hydrolysis, (3) the presence or absence of a gelatinization step, (4) washing/rinsing versus neutralization step, (5) the presence or absence of α-amyloglucosidase in enzymatic solution, (6) the duration of cellulase incubation, and (7) the concentration of the cellulase solution. The major result of our work highlighted that it was essential to carry out a gelatinization step to correctly estimate the in vitro dry matter digestibility of maize silage.


The proposed protocol in this paper is innovative, reliable, highthroughput and easy to implement in many laboratories to accurately quantity in vitro dry matter digestibility.


First introduction of study on maize silage digestibility appeared in France with the work of Auguste Goffart “Manuel de la culture et de ensilage des maïs” [1] for cow winter nutrition and milk production. During the middle of the twentieth century, the concept of digestibility was first introduced by comparing in vivo silage digestibility and in vitro dry matter digestibility (IVDMD) using rumen fluid [2, 3] (Table 1). Tilley and Terry [2] were the first to propose IVDMD quantification following a two step protocol. The first step consists in the action of rumen fluid directly added on the dry matter (DM) and the second step consists in the action of acidic solution and enzymes [2] (Table 1). However, the use of rumen fluid in protocols caused supply difficulties and its utilization was restricted to few teams [4]. To encompass the use of rumen fluid, IVDMD protocols using fungal enzymes appeared in North America during the 70’s (Table 1). Indeed, in 1975, Jones and Hayward [5] proposed a two step digestibility protocol with fungal enzymes (Table 1). The first step of this protocol consists in an enzymatic pretreatment with pepsin (2 g L−1) in acidic solution (HCl 0.1 N) during 24 h at 40 °C and the second step in an acidic digestion with Cellulase (Trichoderma viride) during 48 h at 40 °C (Table 1). To established this protocol Jones and Hayward [5] studied the impact of different parameters such as the influence of different pretreatments, the order of the two step, and the impact of different cellulase origins. Their results on grasses demonstrated that a protocol with two step [step 1: Pepsin; step 2: cellulase solution (T. viride)] was more accurate (r = 0.96) to account for in vivo digestibility values and variation than a protocol without pepsin pretreatment step (r = 0.91). Moreover, they also showed that applying pepsin solution before cellulase cocktail reflected better in vivo digestibility rather than after cellulase cocktail [5].

Table 1 Reported correlations between in vivo and in vitro digestibilities since 1955.

Thereafter, even if the different protocols proposed and used by the scientific community have their own specificities, they all retained two step procedure as demonstrated by the protocols reported in Table 1 (adapted from [4, 6]). Aufrère [6, 7] tested several combinations to adapt the digestibility protocol of Jones and Hayward [5] on temperate grasses. During that experimentation, the impact of the 5 following parameters on IVDMD was studied: the acid concentration during the pretreatment step as well as the duration and temperature of this pretreatment, the cellulase concentration in enzymatic treatment step and the addition of a starch gelatinization step to solubilize the high amount of starch provided by the cob in the case of maize dry matter (Table 1) [6, 8, 9]. The starch gelatinization allowed the breaking down of intermolecular bonds of starch granules and more precisely of amylopectin the more crystalline part of starch.

This has led to the establishment of the “Aufrère protocol” [7] to quantify maize IVDMD at silage stage: this protocol first proceeds to a pretreatment step with a 0.1 N HCl solution with 2% pepsin at 40 °C for 24 h. Afterwards a gelatinization step (80 °C for 30 min) is applied. Finally, after filtration and washing, an enzymatic hydrolysis step is carried out at 40 °C with an Onozuka R10 cellulase solution (1 mg mL−1) for 24 h. The obtained IVDMD is expressed as the mass percentage of matter lost during these successive steps. Resulting IVDMD quantification obtained with this protocol accounts in an acceptable way for more than 60% of in vivo digestibility observed variation [10].

Nowadays, the protocol developed by Aufrère is a benchmark in many laboratories to estimate IVDMD [8, 9, 11, 12] of silage maize. In its original version, this protocol allows to estimate digestibility but it remains restrictive: first of all because 500 mg of DM are digested and this requires high amount of solutions and space to manipulate and secondly because the gelatinization step is followed by a tedious rinsing step before the addition of enzymatic solution, reducing the throughput of the protocol.

In our group we developed a few years ago a highthroughput protocol to characterize starch-free maize samples digestibility. This protocol tardily published [13] was an adaptation of the Aufrère protocol developed in 1982 for starch-free. This protocol first proceeds to a pretreatment step of 30 mg of dry matter with a 0.1 N HCl solution at 40 °C for 24 h. Afterwards a neutralization step was performed with NaOH 2 N before an enzymatic hydrolysis step carried out at 50 °C with an Onozuka R10 cellulase solution (8 mg mL−1) for 72 h.

Herein, we developed a protocol dedicated to starch containing maize samples and capable of enabling highthroughput analyzes. With the aim to facilitate protocol execution, we propose to reduce the dry matter sampling to 30 mg and, instead of rinsing, to neutralize the solution after the gelatinization step as proposed in Virlouvet et al. [13]. We also tested the influence of different parameters which differ significantly between Aufrère and Doreau [7] and Virlouvet et al. [13] (Fig. 1a) forthe estimation of IVDMD of maize silage. We also tested the impact of α-amyloglucosidase in enzymatic solution as proposed by Limagrain [8] (Fig. 1a). For that we compared 21 protocols to test the impact of 7 parameters: (1) presence or absence of pepsin in HCl solution, (2) temperature of incubation during enzymatic hydrolysis, (3) presence (G) or absence (WG) of gelatinization step, (4) washing/rinsing (R) or neutralization (N) step, (5) presence or absence of α-amyloglucosidase in enzymatic solution, (6) duration (24 or 72 h) of cellulase incubation, and (7) concentration (1 or 8 mg mL−1) of the cellulase solution.

Fig. 1
figure 1

A Methodologies for the 3 enzymatic methods selected to test the effect of 7 key parameters (adapted from Andrieu et al. [8]). B Protocol 2 adopted following the results presented in this paper as a mixture between the Aufrère et al. (2007) protocol (pink) and Virlouvet et al. [13] protocol (green)

This study leads to the establishment of an accurate, highthroughput and easy to handle protocol to quantify IVDMD of maize samples.


Reference values of IVDMD

IVDMD reference values were determined by LANO laboratory on a selection of 6 maize dry matter samples that cover a large range of variation for dry matter digestibility. These IVDMD reference values varied from 65 to 76%DM (Fig. 2). Protocol 1 (Table 2) is the one which most closely resembles that proposed by Aufrère [6], taken up by the LANO laboratory.

Fig. 2
figure 2

Plot correlation between value of digestibility protocol of two labs. RV is digestibility reference value of LANO labs Aufrere [7] C1 is Aufrere [7] digestibility protocol adapted at IJBP, black line y = x, red dotted lines is the correlation line, R2 is the determination coefficient, RMSE is root-mean-square error

Table 2 Protocols analyzed

We applied protocol 1 to the 6 selected samples and compared IVDMD obtained values with reference ones. IVDMD values obtained with protocol 1 varied from 67 to 77%DM which faithfully reflects the values obtained by the LANO laboratory. Moreover, between these two protocols, the determination coefficient (R2) is 0.8 and the slope is 1.04 (Fig. 2), indicating that IVDMD values obtained with protocol 1 are similar to the ones obtained by the reference protocol [6].

Gelatinization step is crucial to obtain high similarity of IVDMD with the reference protocol

The impact of 7 parameters was tested by applying on the 6 selected samples the 21 protocols presented in Table 2 which vary for the presence of pepsin, the application of gelatinization, Rinsing or Neutralization, the concentration of cellulase, the presence of amyloglucosidase, the duration of incubation, the temperature of incubation. Determination coefficient, slope and repeatability (cf. “Methods” section) were criteria used to discriminate protocols according to their similarity with the Aufrère protocol [6]. The slope evolved from 0.2 to 1.2 depending protocols. Repeatability is a check of protocol quality and covers a range of 1 to 4% of error except protocol 21 that presents a lower repeatability of with a 6% error.

The dendogram clustered the 21 tested protocols in three distinct groups (Fig. 3). The groups to which Protocol 1 belongs presented highly significant correlations (higher than 0.7) with reference values than the 2 other groups (lower than 0.5) (Figs. 3 and 4). Height of the 9 tested protocols (1, 2, 10, 16, 17, 18, 19 and 20) clustered in this first group are highly correlated together (coefficient correlation R higher than 0.8) (Fig. 4). This group also presented slopes closest to 1. Moreover, protocol 21 belonging to the third group, which showed the highest mean % error (Fig. 3) and was eliminated for further statistical analyses. All other protocols showed good repeatability with less than 5% of mean % error.

Fig. 3
figure 3

Clustering of the 21 protocols according to their R2 and a coefficients. R2 is the determination coefficient between Lano reference values and corresponding values obtained by each tested protocols, a is the slope of the correlation between Lano references value and corresponding values obtained by each tested protocols and mean % error is the mean of the values obtained when divided the standard deviation by the mean

Fig. 4
figure 4

Matrix of correlation between all tested protocols

The regression analysis performed including only main effect parameters yielded in a model with the lowest Akaike Information Criterion (AIC) including the gelatinization step and the temperature of cellulasic hydrolysis incubation (Table 3). In order to clearly depict the impact of these two parameters, R2 of protocols were separated either according to their temperature of incubation (40 or 50 °C; Fig. 5A) or according if a gelatinization step was included in the protocol (Fig. 5B). R2 of protocols were not significantly different if temperature of incubation was either 40 or 50 °C. On the other hand, R2 was significantly increased (pValue = 3.6e10−3) when a gelatinization was included in the protocol, suggesting that gelatinization step is clearly a parameter that needs to be included to increase similarity with the Aufrère protocol on maize silage samples.

Table 3 Multiple regression analyses with the tested model, selected parameters in the tested model, AIC criteria of the model and p-value of each selected parameters
Fig. 5
figure 5

Plot of R2 values between reference values and values obtained for each tested protocols sorting according to temperature of incubation (A), presence or absence of gelatinization step (B), combination of gelatinization and rinsing/neutralization (C) and combination of gelatinization and temperature of incubation (D). p-value between group of protocols: horizontal lines, ***< 0.001, **< 0.01 and *< 0.05

Interactions between gelatinization and rinsing/neutralization and between gelatinization and incubation temperature also impacted IVDMD

Multiple regressions were also performed to pinpoint interactions between parameters that impacted the variation of IVDMD protocol compared to the reference protocol. This analysis identified two interactions that significantly impacted the IVDMD estimation (Table 3).

Firstly, the interaction between gelatinization and rinsing/neutralization was involved significantly in the variation of R2 of the protocols. Indeed, protocols including a gelatinization step and a rinsing step after gelatinization showed significantly higher R2 than protocols without gelatinization and/or with neutralization (Fig. 5C). It is important to note that protocols without gelatinization step and with rinsing step showed similar R2 than model with gelatinization and with neutralization step (Fig. 5C).

Secondly, the interaction between gelatinization step and the incubation temperature also significantly impacted the estimation of IVDMD. In order to better understand the impact of this interaction, Fig. 5D showed that the incubation temperature did not impacted the R2 of the protocols when a gelatinization step was included. However, without gelatinization, the protocols with incubation temperature set up at 50 °C showed higher R2 than protocols with incubation temperature set up at 40 °C (Fig. 5D).


Digestibility protocols have greatly evolved since the middle of the twentieth century. Protocols presented in the literature are very variable. There is no standard in these digestibility estimation protocols and the number of steps, the type or quantity of enzymes or even the incubation times can be very different. In particular, this can make it difficult to compare results from one study to another. Important evolution of these protocols lied in first time by the introduction of two successive steps in the protocols [5, 14] and in second time by the standardization of steps [6, 8, 9, 15]. Some protocols were proposed with a supplemental α-Amylase enzymatic step [16], or McDougall's solution and ruminal fluid (39 °C 48 h) followed by acid pepsin treatment (0.1 N 6.6 g L−1 39 °C 48 h [17]) or only one step of enzymatic hydrolysis by a cellulasic solution [18]. These different protocols gave digestibility values respectively correlated with in vivo digestibility to 0.25, 0.58 and 0.53. However, these correlation with in vivo digestibility data were always weaker than correlation obtained with Aufrère 1982 protocol [6]. This was largely described in the literature and coefficient of correlation between in vitro digestibility values and in vivo digestibility evaluation ranged from 0.9 and 0.6. For example this coefficient was of 0.91 in Aufrère [6], 0.54 in Barrière [15], 0.61 in Andrieu and in Dardenne [8, 9], 0.78 in Aufrère [11] and 0.73 Peyrat [12]. The Aufrère protocol is then considered as a robust and reliable protocol but it remains heavy and bulky to implement.

The in vitro digestibility protocol established by Aufrère [6, 10] is the reference one in digestibility studies of maize silage [19]. Moreover in France, combined with protein content, IVDMD values obtained with the Aufrère protocol are used to estimate UFL (Unité Fourrage Laitier) value with model M4 since 1995 [20] and since 2016 with model M4.2 slightly adapted from M4 model [21]. UFL value is essential for the registration process of maize silage varieties in French maize silage catalog [22]. However, the Aufrère protocol does not allow a highthroughput characterization of dry matter samples. Some attempts to adapt the Aufrère protocol to achieve a high characterization rate have been made. The use of Daisy technology have been performed [23]. The Daisy II method can be used on feeds to increase labor efficiency [24] but from our point of view, it still presents major drawbacks. Indeed, this technology always necessitates high quantity of dry matter samples and large volumes of enzymatic solutions. Moreover, the filter bags used in daisy technology can present variations in their porosity depending on the batch used which can cause undigested particle leaks [25]. Furthermore, all the sample is bathed in the same solution and we can always fear interactions between sample and inhibitions on enzymatic activities related to a particular sample. We thus proposed herein, a protocol with smallest dry matter sampling (30 mg) than Aufrère protocol with 500 mg of DM [10].

We demonstrated that few protocols were highly correlated with LANO reference values. With this new format we avoid some problems like variable bag porosity or digestive samples interactions, and it’s broadband, repeatable and comparable to Aufrère protocol. We have clearly underlined the importance of two major parameters in the implementation of these digestibility estimations.

First, we highlighted the need to perform a gelatinization step. This additional step was proposed by Aufrère between 1982 and 1983 [6, 7]. This gelatinization step allowed starch gelatinization i.e. the breaking down of intermolecular bonds of starch granules and more precisely of amylopectin the more crystalline part of starch. Gelatinization step thus allowed a better accessibility of starch molecules for subsequent hydrolysis.

The second important step we revealed is the rinsing step between acidic pretreatment and Onozuka enzymatic hydrolysis step. The combination of these two step (starch gelatinization and rinsing) was the best one to obtain values similar to the reference values. Strictly speaking, we should have retained protocol 16 or 17 as the best one. However, practical aspects allowing highthroughput analyzes also guided our choice. Thus a neutralization step was preferred to rinsing steps even if the results were a little less faithful to the references. This is why protocol 2 was chosen rather than protocol 16 or 17. Protocol 2 combined good statistical performances and presented a high level of practicality for its highthroughput implementation in the laboratory. This combination aimed reducing the arduousness and the time of manipulation while giving very reliable results. Therefore we retained the protocol 2 as the more suitable to easily and accurately measure dry matter digestibility.


Protocol 2 retained in this study first proceeds to a pretreatment step of 30 mg of dry matter with 0.1 N HCl solution with 2 g L−1 pepsin at 40 °C for 24 h. Afterwards a gelatinization step (80 °C for 30 min) is applied. Finally, after a neutralization step, an enzymatic hydrolysis step is carried out at 40 °C with an Onozuka R10 cellulase solution (1 mg mL−1) for 24 h (Fig. 1b).

This protocol is innovative on several points compared to the reference protocol of Aufrère and Doreau [7]: the reduction of the dry matter sampling necessary at the start (30 mg versus 500 mg) as well as the neutralization step between the first and the second step, which reduces handling time and improves ease of handling (Fig. 1a, b). This protocol 2 is neither more nor less than the protocol developed by Aufrère and Doreau in 1983 [7] in which the dry matter sampling is reduced to 30 mg, the rinsing replaced by a neutralization step and the glass filtering crucibles replaced by tubes. This allows us to perform experiments with a very large number of samples and reduce volumes of enzymatic solution without hampering the accuracy of the IVDMD quantification.

This protocol is reliable and can be implemented in any laboratory without any specific equipment. It makes it possible to simultaneously characterize a large number of samples, with a small amount of dry matter, with low volumes of solution and without tedious rinsing and filtration steps.


Plant materials

Plants of maize hybrids were cultivated in open field trial in summer 2018. 80 genotypes were grown in two 7 m rows with 0.80 m between rows for 11.2 m2 by genotypes and planting density of 90,000 plants/ha. At the silage stage, plants were harvested, dried in an oven (55 °C—72 h) and grind with a hammer mill (1 mm grid). From the 80 harvested samples (one sample per genotype), we selected 6 contrasting samples that cover the range of dry matter digestibility based on results obtained using protocol 1 (Table 2). This set of 6 samples has been used for the analyses carried out in this paper.

Reference values for IVDMD estimation

Estimation of the IVDMD of the 6 selected samples was performed by the reference laboratory LANO (Laboratoire Agronomique de Normandie following the initial reference protocol [6, 10]. Briefly, 50 mL of HCl 0.1 N solution with 2% of pepsin was applied on 500 mg of maize dry matter during 24 h at 40 °C. After that, samples were placed 30 min at 80 °C (gelatinization step). Rinsing was performed before adding 50 mL of enzymatic solution of Cellulase Onozuka («Onozuka R-10 from Trichoderma viride»; Serva; 10 mg ref 16,419.03) at 1 mg mL−1 during 24 h at 40 °C. A final rinsing was then performed before 12 h drying at 103 °C.

Estimation of IVDM using modified Aufrère protocol

We adapted the Aufrère protocol [10] to our lab conditions and to smaller dry matter sample quantity. This adapted protocol is called Protocol 1 (Table 2) for the rest of the article. Briefly, for this Protocol 1, 30 mg of maize dry matter powder were placed in 5 mL tubes before adding 2 mL of HCl 0.1 N solution with 2% of pepsin during 24 h at 40 °C followed. Then, gelatinization step was carried out by placing the tubes during 30 min at 80 °C in a water bath. The tubes were cooled in ice before being centrifuged at room temperature 10 min at 5000 rpm. The supernatant was discarded. Rinsing was then performed by centrifuging twice with 4 mL of water (5000 rmp—10 min at room temperature) the obtained pellet. After the first rinse centrifugation, the supernatant was discarded and the pellet was rinsed a second time with 4 mL of water. Similarly, the supernatant obtained after the second rinse centrifugation was discarded. 4 mL of Enzymatic solution of Cellulase («Onozuka R-10 from Trichoderma viride»; Serva; 10 mg ref 16,419.03) at 1 mg L−1 were then added to the pellet, vortexed and agitated during 24 h at 40 °C. The pellet recovered after centrifugation (5000 rmp—10 min at room temperature) was then freezed at − 80 °C and lyophilized 48 h before final weighing.

Protocols establishment to estimate the impact of different parameters on IVDMD estimation

We tested several protocols to study the influence of 7 parameters on the quantification of dry matter digestibility (Table 2). These 7 tested parameters are as follows:

  • presence (2 g L−1) or absence of pepsin in HCl solution (0.1 N)

  • temperature of incubation (40 or 50 °C) during cellulasic hydrolysis

  • presence (G: 30 min at 80 °C) or absence (WG) of gelatinization step

  • washing/rinsing (R) or neutralization (N) step. In the case of washing/rinsing alternative, a succession of 3 centrifugation step (5000 rmp, 10 min at room temperature) was performed. The supernatants were discarded. Afterwards, 4 mL of cellulasic solution (concentration according to the protocol) at 4.6 pH were added to the pellet and vortexed. In the case of neutralization alternative, basic solution (90 µL of NaOH solution 2 N) was added in the tube and the tubes were then vortexed.

  • presence (1.5 mL L−1) or absence of α-amyl glucosidase in cellulasic solution

  • duration (24 or 72 h) of cellulasic solution incubation

  • concentration (1 or 8 mg mL−1) of the cellulasic solution.

All these protocol variants led to the establishment of 21 protocols presented in Table 2. The six selected dry matter samples were digested using these 21 protocols with 3 replicates per protocol.

Statistical analyses

Means of IVDMD for each sample and for each protocol have been calculated. The repeatability of each protocol was estimated by the mean of the percentage of error (where error = standard deviation/mean) of each sample by protocol.

We analyzed the correlation between each protocol by matrix correlation [26]. Moreover, means of IVDMD values per sample of each protocol were correlated to IVDMD estimation from the LANO laboratory (references values). This allowed us, for each protocol, to obtain a determination coefficient (R2) and a slope (“a” coefficient from the equation Y = aX + b, where Y is IVDMD (references values) value and X is IVDMD measured with one of the 21 established protocols (Table 3). R2 and slope are reliable indicators to discriminate protocols according to their similarity to the IVDMD reference values from the LANO Laboratory that used the original Aufrère (1982) protocol. The more similar the protocol are and the more the R2 values are high and the slope closer to 1. We use the ComplexHeatmap package Rstudio to cluster and visualize the tested protocols according to their R2 slope and repeatability values (Fig. 3).

Two multiple regressions analyses were performed using R package “stats” with the “step” functions to.

  • analyze main effects of each parameter using the following model:

    $${\text{Ylmnopqr}} =\, \upmu + {\text{Cl}} + {\text{Dm}} + {\text{Fn}} + {\text{Go}} + {\text{Hp}} + {\text{Iq}} + {\text{Jr}} + {\text{Elmnopqr}}$$
  • analyze interactions between parameters using the following model:

    $${\text{Ylmnopqr}} =\, \upmu + {\text{Cl}}*{\text{Dm}}*{\text{Fn}}*{\text{Go}}*{\text{Hp}}*{\text{Iq}}*{\text{Jr}} + {\text{Elmnopqr}}$$

    where Ylmnopqr is the R2 value of the pepsin concentration “l”, the presence of the gelatinization step “m”, the rinsing or neutralization modality “n”, the cellulasic concentration “o”, the amyloglucosidase concentration “p”, the incubation times “q” and the incubation temperatures “r”. Cl is the main effect of the pepsin concentration “l”; Dm is the main effect of the gelatinization step “m”; Fn is the main effect of the Rinsing or Neutralization “n”; Go is the main effect of the cellulasic concentration “o”; Hp is the main effect of the amyloglucosidase concentration “p”; Iq is the main effect of the incubation times “q”; Jr is the main effect of the incubation temperatures “r”; Elmnopqr is the random residual term.

Wilcoxon test were carried out using R package ggplot2 fonction stat_summary [27] to highlight significant differences between R2 values of protocols.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



In vitro dry matter digestibility


Dry matter


Unité Fourrage Laitier


  1. Goffart A. Manuel de la culture et de l’ensilage des maïs et autres fourrages verts. 1877.

  2. Tilley JMA, Terry RA. A two-stage technique for the in vitro digestion of forage crops. Grass Forage Sci. 1963;18:104–11.

    Article  CAS  Google Scholar 

  3. Van Soest PJ, Wine RH, Moore LA. Estimation of the true digestibility of forages by the in vitro digestion of cell walls. 1966. Accessed 24 Nov 2020.

  4. Marten GC, Barnes RF. Prediction of energy digestibility of forages with in vitro rumen fermentation and fungal enzyme systems. IDRC, Ottawa, ON, CA; 1979. Accessed 5 May 2020.

  5. Jones DIH, Hayward MV. The effect of pepsin pretreatment of herbage on the prediction of dry matter digestibility from solubility in fungal cellulase solutions. J Sci Food Agric. 1975;26:711–8.

    Article  CAS  Google Scholar 

  6. Aufrère J. Etude de la prévision de la digestibilité des fourrages par une méthode enzymatique. Ann Zootech. 1982;31:111–30.

    Article  Google Scholar 

  7. Aufrère J, Michalet-Doreau B. In vivo digestibility and prediction of digestibility of some by-products. In: EEC Seminar. 1983.

  8. Andrieu J, Demarquilly C, Dardenne P, Barrière Y, Lila M, Maupetit P, et al. Composition and nutritive value of whole maize plants fed fresh to sheep. I: factors of variation. Compos Nutr Value Whole Maize Plants Fed Fresh Sheep Factors Var. 1993;42:221–49.

    Google Scholar 

  9. Dardenne P, Andrieu J, Barrière Y, Biston R, Demarquilly C, Femenias N, et al. Composition and nutritive value of whole maize plants fed fresh to sheep. II. Prediction of the in vivo organic matter digestibility. Ann Zootech. 1993;42:251–70.

    Article  Google Scholar 

  10. Aufrère J, Baumont R, Delaby L, Peccatte JR, Andrieu J, Andrieu JP, et al. Prévision de la digestibilité des fourrages par la méthode pepsin-cellulase. Le point sur les équations proposées. INRAE Prod Anim. 2007;20:129–36.

    Article  Google Scholar 

  11. Aufrere J. Inra Productions Animales - Prévision de la digestibilité des fourrages. INRA Prod Anim. 2007;20:129–36.

    Article  CAS  Google Scholar 

  12. Peyrat J. Digestion de l’amidon et des parois végétales du maïs fourrage chez les ruminants: conséquences sur l’évaluation de sa valeur nutritive. France: Université Blaise Pascal: Clermont-Ferrand II; 2014.

  13. Virlouvet L, El Hage F, Griveau Y, Jacquemot M-P, Gineau E, Baldy A, et al. Water deficit-responsive QTLs for cell wall degradability and composition in maize at silage stage. Front Plant Sci. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Alexander RH, McGowan M. The routine determination of in vitro digestibility of organic matter in forages—an investigation of the problems associated with continuous large-scale operation. Grass Forage Sci. 1966;21:140–7.

    Article  CAS  Google Scholar 

  15. Barrière Y, Demarquilly C, Hebert Y, Dardenne P, Andrieu J, Maupetit P, et al. Influences de la variabilité génétique et environnementale sur la digestibilité in vitro ou in vivo du maïs fourrage. Agronomie. 1991;11:151–67.

    Article  Google Scholar 

  16. Lila M, Yves B, Rémi T. Mise au point et étude d’un test enzymatique de la digestibilité de fourrages pauvres ou riches en amidon. Agronomie. 1986;6:285–91.

    Article  Google Scholar 

  17. Damiran D, DelCurto T, Bohnert DW, Findholt SL. Comparison of techniques and grinding size to estimate digestibility of forage based ruminant diets. Anim Feed Sci Technol. 2008;141:15–35.

    Article  Google Scholar 

  18. Kowalski ZM, Ludwin J, Górka P, Rinne M, Weisbjerg MR, Jagusiak W. The use of cellulase and filter bag technique to predict digestibility of forages. Anim Feed Sci Technol. 2014;198:49–56.

    Article  CAS  Google Scholar 

  19. Mabjeesh SJ, Cohen M, Arieli A. In vitro methods for measuring the dry matter digestibility of ruminant feedstuffs: comparison of methods and inoculum source. J Dairy Sci. 2000;83:2289–94.

    Article  CAS  Google Scholar 

  20. Andrieu J. Prévision de la digestibilité et de la valeur énergétique du maïs fourrage à l’état frais. INRA Prod Anim. 1995;8:273–4.

    Article  Google Scholar 

  21. Peyrat J, Nozière P, Férard A, Le Morvan A, Meslier E, Protin PV, et al. Digestibilité du maïs fourrage: équations et critères de prédiction. ARVALIS - Inst Végétal - INRA. 2016. Accessed 21 Dec 2018.

  22. VATE. Evaluation de la Valeur Agronomique Technologique et Environnementale (VATE) des nouvelles variétés à l’inscription au Catalogue Français. 2018.

  23. Tassone S, Fortina R, Peiretti PG. In vitro techniques using the DaisyII incubator for the assessment of digestibility: a review. Animals. 2020;10:775.

    Article  Google Scholar 

  24. Holden LA. Comparison of methods of in vitro dry matter digestibility for ten feeds. J Dairy Sci. 1999;82:1791–4.

    Article  CAS  Google Scholar 

  25. Adesogan AT. Effect of bag type on the apparent digestibility of feeds in ANKOM DaisyII incubators. Anim Feed Sci Technol. 2005;119:333–44.

    Article  Google Scholar 

  26. Wei T, Simko V, Levy M, Xie Y, Jin Y, Zemla J, et al. Visualization of a correlation matrix. 2021. Accessed 26 May 2021.

  27. Wickham H, Chang W, Henry L, Lin Pedersen T, Takahashi K, Wilke C, et al. Create elegant data visualisations using the grammar of graphics. 2020., Accessed 26 May 2021.

  28. Thurman RL, Wehunt EJ. Laboratory method for determining digestible nutrients. Agron J. 1955;47:302–3.

    Article  CAS  Google Scholar 

  29. Clark KW. The adaptation of an artificial rumen technique to the estimation of the gross digestible energy of forages. Theses Diss Available ProQuest. 1958;1–110.

  30. Donefer E, Niemann PJ, Crampton EW, Lloyd LE. Dry matter disappearance by enzyme and aqueous solutions to predict the nutritive value of forages. J Dairy Sci. 1963;46:965–70.

    Article  CAS  Google Scholar 

  31. Deinum B, Van Soest PJ. Prediction of forage digestibility from some laboratory procedures. NJAS Wagening J Life Sci. 1969;17:119–27.

    CAS  Google Scholar 

  32. Jones DIH, Hayward MV. A cellulase digestion technique for predicting the dry matter digestibility of grasses. J Sci Food Agric. 1973;24:1419–26.

    Article  CAS  Google Scholar 

  33. Hartley RD, Jones EC, Fenlon JS. Prediction of the digestibility of forages by treatment of their cell walls with cellulolytic enzymes. J Sci Food Agric. 1974;25:947–54.

    Article  CAS  Google Scholar 

  34. McQueen R, Van Soest PJ. Fungal cellulase and hemicellulase prediction of forage digestibility. J Dairy Sci. 1975;58:1482–91.

    Article  CAS  Google Scholar 

  35. Adegbola AA, Paladines O. Prediction of the digestibility of the dry matter of tropical forages from their solubility in fungal cellulase solutions. J Sci Food Agric. 1977;28:775–85.

    Article  CAS  Google Scholar 

  36. Clark J, Beard J. Prediction of the digestibility of ruminant feeds from their solubility in enzyme solutions. Anim Feed Sci Technol. 1977;2:153–9.

    Article  CAS  Google Scholar 

  37. Dowman MG, Collins FC. The prediction of the digestibility of silages using cellulase. J Sci Food Agric. 1977;28:1071–4.

    Article  CAS  Google Scholar 

  38. Goto I, Minson DJ. Prediction of the day matter digestibility of tropical grasses using a pepsin-cellulase assay. Anim Feed Sci Technol. 1977;2:247–53.

    Article  Google Scholar 

  39. Roughan PG, Holland R. Predicting in-vivo digestibilities of herbages by exhaustive enzymic hydrolysis of cell walls. J Sci Food Agric. 1977;28:1057–64.

    Article  CAS  Google Scholar 

  40. Kellner RJ, Kirchgessner M. Estimation of forage digestibility by a cellulase method. Z Für Tierphysiol Tierernähr Futtermittelkunde. 1977;39:9–15.

    Article  CAS  Google Scholar 

  41. Kirchgessner M, Kellner RJ. Estimation of digestibility, metabolizable energy and net energy of forage by a cellulase method. Livest Prod Sci. 1978;5:373–7.

    Article  CAS  Google Scholar 

  42. McLeod MN, Minson DJ. The accuracy of the pepsin-cellulase technique for estimating the dry matter digestibility in vivo of grasses and legumes. Anim Feed Sci Technol. 1978;3:277–87.

    Article  CAS  Google Scholar 

  43. Adamson AH, Terry GR. The relationship between the in vivo digestibility of hay and its solubility in pepsin–hydrochloric acid and fungal cellulase solutions. J Sci Food Agric. 1980;31:854–6.

    Article  CAS  Google Scholar 

  44. Aufrere J, Michalet-Doreau B. Comparison of methods for predicting digestibility of feeds. Anim Feed Sci Technol. 1988;20:203–18.

    Article  Google Scholar 

  45. De Boever JL, Cottyn BG, Andries JI, Buysse FX, Vanacker JM. The use of a cellulase technique to predict digestibility, metabolizable and net energy of forages. Anim Feed Sci Technol. 1988;19:247–60.

    Article  Google Scholar 

  46. Ferret A, Gasa J, Plaixats J, Casañas F, Bosch L, Nuez F. Prediction of voluntary intake and digestibility of maize silages given to sheep from morphological and chemical composition, in vitro digestibility or rumen degradation characteristics. Anim Sci. 1997;64:493–501.

    Article  CAS  Google Scholar 

Download references


Maize silage samples grinding was performed by Dominique Kermarrec INRAE—Unité Expérimentale sur les Ressources Génétiques Végétales en Conditions Océaniques.


The IJPB benefits from the support of Saclay Plant Sciences-SPS (ANR-17-EUR-0007). This work was financially supported by Pro-maïs and is part of the DECLIC project ( Field trail was also carried out in the frame of DECLIC. This study was supported by the French Government Grant ANR-11-BTBR-0006 BIOMASS FOR THE FUTURE, managed by the French National Research Agency under an “Investments for the Future” program (Grant ANR-11-IDEX-0003-02).

Author information

Authors and Affiliations



VM and MR piloted the study and supervised the plant sampling; PLLM, MPJ, MR and VM sampled the maize, PLLM, SG and MPJ performed the biochemical analyses; PLLM and MR performed statistical analyses on the dataset; PLLM, MR and VM interpreted the results; PLLM, MR and VM wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to V. Méchin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lopez-Marnet, PL., Guillaume, S., Jacquemot, MP. et al. High throughput accurate method for estimating in vitro dry matter digestibility of maize silage. Plant Methods 17, 89 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: