J. Anim. Sci. 2003. 81:520-528
© 2003 American Society of Animal Science
Simplified in situ method for estimating ruminal dry matter and protein degradability of concentrates1
V. Olaisen*,
T. Mejdell
,
H. Volden
,1 and
N. Nesse*,2
* Department of Chemical Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim;
and
SINTEF Applied Chemistry, Chemical Engineering, NO-7465 Trondheim, Norway; and
and
Department of Animal Science, Agricultural University of Norway, NO-1432 As
2 Correspondence:
P.O. Box 5025 (phone: +47-64-94-80; fax: +47-64-94-79-60, E-mail:
harald.volden{at}ihf.nlh.no).
 |
Abstract
|
|---|
In this study, dry matter and crude protein in situ degradation data from different concentrate feeds were used to test the accuracy of effective degradability (ED) measures when using reduced ruminal incubation times compared with models based on seven or eight incubation times. The ED was estimated both with and without correction for nylon bag particle loss. The crude protein ED corrected for particle loss of the calibration data set was widely distributed in a range from 16 to 90% with an overall mean value of 60.4%, and the dry matter ED was distributed in the range from 22.7 to 80.7%, with a mean value of 56.9%. The simplified method was developed based on bilinear regression models where all combinations of one to three disappearance values were tested to find the optimal time point combinations to estimate ED. Bilinear regression models based on two and three ruminal incubation times gave similar estimates to a standard in situ method over a wide range of passage rates both for the data set used to parameterize the models and the independent data set used to evaluate the models. Using two incubation times, the bilinear model based on 4 and 24 h gave the most accurate estimates, and the models based on 2, 8, and 24 h for uncorrected data and 4, 8, and 24 h for corrected data were most accurate of the three time points bilinear models. The number of nylon bags used by these models was reduced by 58 to 78% compared with the standard in situ method, and the total incubation time needed was substantially reduced.
Key Words: Cattle Degradation Dry Matter Mathematical Model Protein Rumen
 |
Introduction
|
|---|
Several methods are used to estimate ruminal degradation of DM and CP (Hvelplund and Weisbjerg, 1998). The in situ method has achieved the widest use. In spite of a number of limitations, no in vitro method has been generally accepted as a satisfactory alternative (Stern et al., 1997; Hvelplund and Weisbjerg, 1998). Therefore, the in situ technique has been chosen as the reference method to measure rumen CP degradation in several protein evaluation systems for ruminants (Madsen, 1985; Vérité and Peyraud, 1989; Tamminga et al., 1994). However, the method requires a large number of nylon bags to be ruminally incubated for each feed sample and, in turn, a substantial amount of human work. Therefore, a need exists for simpler in situ techniquesin particular, for routine analysis (Vanzant et al., 1996). Some simplified in situ techniques have been proposed for degradability measurements. For example, Wilkerson et al. (1995) used only a 16-h ruminal incubation time to estimate the ruminal CP degradability in roughages by assuming that the 16-h sample directly estimated the escape CP, whereas Broderick (1994) and Calsamiglia et al. (1994) calculated effective degradability (ED) by estimating the ruminal degradation rate from double-point incubations. However, questions remain as to whether it is possible to develop a simplified in situ method for concentrate feeds that can be used to estimate ED directly without losing accuracy. Therefore, the objective of this study was to evaluate a novel mathematical approach to calculate ED of CP and DM in concentrate feeds from in situ studies using a minimal number of incubation times. We evaluated these models for in situ data both corrected and not corrected for small particle loss of the wash fraction.
 |
Materials and Methods
|
|---|
Feeds
Data from an earlier study by Volden and Harstad (1995) and from different unpublished in situ studies at the Agricultural University of Norway were used to develop the models. Descriptions of the feeds used in the present study to develop the simplified models are shown in Tables 1 and 2
. The material consisted of 81 individual feed ingredients and 54 concentrate mixtures. The concentrate mixtures consisted mainly of barley and oats treated with an annular gap expander (OE expander; Amandus Kahl, Reinbek, Germany) in the range of 110 to 130°C.
View this table:
[in this window]
[in a new window]
|
Table 1. Water solubility of dry matter, in situ dry matter wash fraction, and uncorrected and corrected effective degradability (ED) for different concentrates
|
|
View this table:
[in this window]
[in a new window]
|
Table 2. Water solubility of crude protein, in situ crude protein wash fraction, and uncorrected and corrected effective degradability (ED) for different concentrates
|
|
To evaluate the simplified models, we used an independent data set (H. Volden, unpublished data) consisting of 32 individual feed ingredients (Table 3
). In these feeds, there was no correction for small particle loss of the wash fraction.
View this table:
[in this window]
[in a new window]
|
Table 3. Dry matter and crude protein in situ effective degradability of the feeds used to evaluate the simplified models
|
|
In situ Measurements
Three nonlactating cows fed a diet consisting of two-thirds hay and one-third concentrate mixture on DM basis at maintenance level (Madsen et al., 1995) were used to obtain the in situ disappearance profiles of DM and CP. Feed samples were ground through a 1.5-mm screen. Approximately 2 g of feed was placed into nylon bags (6 x 12 cm, 13 to 15 mg/cm2) with a pore size of 36 µm (ZBF AG, CH 8803, Rüschlicon, Switzerland). Nylon bags were incubated in the rumen for 2, 4, 8, 16, 24, and 48 h. For 44 feed samples, nylon bags were additionally incubated for 72 h. The number of nylon bags used varied by incubation time and cow: three (0 and 2 h), four (4, 8, 16, and 24 h), six (48 h), and eight (72 h). The rationale for increasing the number of bags with increased incubation time was to ensure adequate residue for subsequent analysis. Both DM and CP disappearance values were measured for each cow. The in situ 0-h wash fraction of DM and CP was measured by washing nylon bags containing samples in a domestic washing machine for 3 x 10 min in cold water. Water solubility of DM and CP was obtained by soaking a 1-g sample of feed prepared for in situ measurements in 40 mL of distilled water in a 100-mL centrifuge tube at 20°C for 1 h with shaking every 10 min. The sample was then filtered through a Schleicher & Schuell 5891 black ribbon filter (Ref. No. 300-011) and further washed with 3 x 40 mL of distilled water. After washing, the filter was dried at 45°C for 48 h and then weighed and analyzed for N in the feed residue. The ruminal incubation program and washing and drying procedures were completed as described by Volden and Harstad (1995).
All cows were cared for according to laws and regulations controlling experiments on live animals in Norway (i.e., the Animal Protection Act of December 20, 1974, and the Animal Protection Ordinance concerning Experiments on Animals of January 15, 1996).
Estimating Effective Degradability by the Reference Method
Effective degradability was estimated with or without correction for particle loss of the wash fraction. The wash fraction was assumed to be totally and instantaneously degraded when not corrected for particulate loss. The particulate loss was measured as the difference between the in situ 0-h wash fraction and the water-soluble fraction. Three assumptions were made for the properties of the particle loss fraction (Dhanoa et al., 1999): 1) it consists of a degradable and undegradable fraction of the same ratio as the particles remaining in the bag after washing; 2) the degradation rate was the same for the particle loss fraction as for the particles remaining in the bag; and 3) the passage rate was the same for the particle loss fraction as for the particles remaining in the bag.
These three assumptions give the following equation to adjust the experimental disappearance values when correcting for particle loss in the wash fraction:
 | [1] |
where D(ti) and Dcorr (ti) are the measured and corrected percentage disappearance at time ti, s is the water solubility (%), and w is the wash fraction (%).
The in situ residues were not corrected for microbial contamination. Although this will affect the absolute degradation values, it will not change the principles in development of simplified in situ methods and their evaluation.
In the reference method, with a full-time incubation program, the kinetics of the in situ DM and CP disappearance were calculated by a generalized Von Bertalanaffy model (López et al., 1999) expressed as:
 | [2] |
where D(t) is the percentage of disappearance from the bags at time t, A and B are the soluble and potentially degradable fractions, and c is the first-order degradation rate. The v parameter will adjust the shape of the curve where low values result in a sigmoid shape (Figure 1
, SoyPass), high values result in an initially very steep curve (Figure 1
, barley), and a v = 1 gives a first-order kinetics. This model was used both for the uncorrected and corrected experimental disappearance values. The nonlinear parameters in both models were estimated from a Levenberg-Marquardt algorithm implemented in MATLAB (The MathWorks, Inc., Natick, MA).

View larger version (19K):
[in this window]
[in a new window]
|
Figure 1. Dry matter ( ) and crude protein () disappearance values (not corrected for particle loss of the wash fraction) and their respective curves fitted with the equation D(t) = A + B(1 - e-ct)1/v.
|
|
The ED was calculated by Fogler (1992) as:
 | [3] |
where the r subscript denotes the reference method and E (t) is the rumen residence time distribution function (RTD) of the feed particles. The RTD function is a normalized, dimensionless ruminal outflow concentration profile, where the total area under the curve is equal to 1. Equation [3]
may be used for any combination of rumen degradation (D[t]) and RTD model (E[t]). We used the RTD for an ideally mixing compartment, E(t) = kpe-kpt, where kp is the passage rate (h-1) of feed particles out of the rumen. In the calculations of ED, the passage rate was set at 0.06 h-1 (if no other value was specified), which is used for concentrates in several protein-evaluating systems (Vérité and Peyraud, 1989; Tamminga et al., 1994). The disappearance and RTD models were inserted into Eq. [3]
and integrated numerically to estimate effective degradability.

View larger version (13K):
[in this window]
[in a new window]
|
Figure 3. Correlation plot for effective degradability calculated by the reference method (EDr) and effective degradability calculated using the three incubation times (2, 8, and 24 h) and the bilinear regression method (EDb) at a passage rate of 0.06 h-1 for the data set used to derive the bilinear models. Data were not corrected for particle loss of wash fraction.
|
|
Simplified Method: Estimating Effective Degradability Using a Bilinear Regression Model
Effective Degradability was estimated by simple bilinear models, generally described by:
 | [4] |
where the b subscript denotes the bilinear model, a0 and ai are the parameters, D(ti) is the percentage of DM or CP disappearance from the bags at the times ti, and n is the number of incubation times used in the model. For uncorrected data, this general model was used to test all one- (7), two- (21), and three-timepoint (35) model combinations from 0 (wash sample) to 48 h to determine the timepoints that were optimal for estimating ED. Both the wash and the water-soluble fractions had to be measured when using corrected disappearance values (Eq. [1]
), therefore, the water solubility was incorporated in all bilinear models using disappearance values corrected for particle loss of the wash fraction:
 | [5] |
where Dcorr (ti) is the corrected percentage of DM or CP disappearance, a0, a1, and ai are the parameters, and s is the water solubility (%). The bilinear models were parameterized separately for the DM and CP D(t) data sets and for a data set containing both the DM and CP D(t) values.
The accuracy relative to the reference model of the different bilinear models was evaluated by calculating the root mean square error (RMSE):
 | [6] |
 |
Results and Discussion
|
|---|
Effective Degradability Estimated by the Reference Method
Tables 1 and 2
show the mean values for DM and CP water solubility, in situ wash fraction, and estimated ED for the feeds used to parameterize the bilinear models, respectively. There was a large difference between water solubility and wash fraction for many of the feedsin particular, for oats where the wash machine loss for DM and CP was more than 45 percentage units higher than the water solubility. These large differences also resulted in large differences between corrected and uncorrected degradability estimates, in particular, for oats, concentrate mixtures, and peas. The mean ED varied largely between the different feeds. Calculated with a passage rate of 0.06 h-1, the uncorrected CP ED values in the data set used to parameterize the models ranged from 23 to 95%, with a mean value of 70.2%. The corrected CP ED values ranged from 16 to 90%, with a mean value of 60.4%. Also, the shape of the disappearance curves varied widely, and Figure 1
illustrates some of the disappearance curve variation in the data set. The v parameter in the generalized Van Bertalanaffy model (Eq. [2]) made this model very flexible, and curves with both slow and steep initial degradation rate were fitted well by this model (Figure 1
).
Table 3
shows DM and CP uncorrected mean ED in the independent data set used to test the bilinear models. The total mean CP ED of this data set was 70.4% and ranged from 31.0 to 89.8%, which was similar to the data set used to parameterize the models.
Effective Degradability Estimated by Bilinear Regression Models
Table 4
shows RMSE and r2 values for the one, two, and three ruminal incubation time models with the lowest RMSE values compared to the reference method. Using a passage rate of 0.06 h-1, more than 97% of the variation in ED was explained when the one- incubation time (8 h) bilinear model was used. For passage rates lower than 0.05 h-1, a model based on the 16-h incubation time showed the lowest RMSE value of the one-timepoint models (data not shown). However, with 17% of the residuals larger than 4 percentage units (uncorrected model) in absolute value, the errors seemed to be too large for using only a single timepoint model in routine analyses. Using a two-timepoint model, only 2% of the residuals were larger than 4 percentage units in absolute value, and the model based on 4- and 24-h incubation times generally had the lowest RMSE of the two-point model across the passage rates tested (Figure 2
). Using uncorrected data, the model based on 2-, 8-, and 24-h incubation times had the lowest RMSE of the three-timepoint models for passage rates higher than 0.047 h- (Figures 2 and 3
), and the model based on 2, 16, and 48 h had lowest RMSE for slower passage rates. Using data corrected for particle loss, the model based on 4-, 8-, and 24-h incubation times had the lowest RMSE of the three-timepoint models. The parameters in the two- and three-timepoint models varied with the passage rate and are presented in Tables 5 and 6
. These parameters can be combined with in situ disappearance values to estimate ED directly if using the same incubation times as in Tables 5 and 6
. The parameters can be estimated by linear interpolation if the same incubation times and passage rates between the values in Tables 5 and 6
are used. These parameters can also be used for data corrected or not corrected for microbial contamination. However, new parameters must be estimated if different incubation times than those presented in Table 5 and 6
are used. These parameters were estimated using both DM and CP data. The parameters and RMSE values were only slightly different when the parameters were estimated using DM and CP data separately. As we discuss later, the bilinear regression method may be considered a simple numerical integration of the degradation Eq. [3]. Since these parameters represent the area sections under the E(t) curve, it is reasonable to assume that the parameters will be similar for different feed fractions, such as DM and CP. The parameters presented in this paper can probably also be applied to other feed fractions, such as starch, NDF, and individual amino acids, although the models presented in this study were developed for DM and CP data only. The current study deals strictly with concentrate feeds. Forages have slower passage rates than concentrates (0.02 to 0.05 h-1; Poncet et al., 1995). Figure 2
shows that the three-point bilinear model based on 2, 16, and 48 h gave low RMSE values at low passage rates. We therefore believe that a similar approach will work well with forages too.
View this table:
[in this window]
[in a new window]
|
Table 4. Root mean square error (RMSE), coefficient of determination, and percentage reduction in the number of nylon bags by different bilinear regression models
|
|
View this table:
[in this window]
[in a new window]
|
Table 5. Parameters in two and three ruminal incubation times for different passage rates of bilinear regression models not corrected for particle loss of wash fractiona
|
|
View this table:
[in this window]
[in a new window]
|
Table 6. Parameters in two and three ruminal incubation times for different passage rates of bilinear regression models corrected for particle loss of wash fractiona
|
|
Using the parameters in Table 5
on the independent data set gave ED estimates very similar to those estimated by the reference method (Figure 4
). The accuracy relative to the reference method for the two- and three-timepoint bilinear models was nearly equal to the accuracy found in the data set used to parameterized the models. For the three-timepoint model, the r2 values were even higher for the evaluation data set than the data set used to parameterize the models. The predictability of ED using the bilinear models was therefore good.
Poncet et al. (1995) reported that the passage rate of labeled concentrate particles is usually in the range of 0.04 to 0.06 h-1 in dairy cows. The French (Vérité and Peyraud, 1989) and Dutch (Tamminga et al., 1994) protein-evaluating system used a passage rate of 0.06 h-1 for concentrates, whereas the Nordic AAT/PBV system used passage rates of 0.04 to 0.08 h-1, depending on the country (Madsen et al., 1995). Both two- and three-timepoint bilinear regression models fitted the data well in the actual passage rate range (0.04 to 0.08h-1). The estimates relative to the standard method were substantially more precise for the three-timepoint model than for the two-timepoint model (Figure 2
), and RMSE did not vary much for the three-timepoint model within the range of realistic passage rates evaluated. Therefore, the question of which model to choose is one of what to stress the most: simplicity or accuracy. The bilinear regression models substantially reduced the number of nylon bags compared to the reference in situ method (Table 4
). Therefore, the amount of manual human work and, in turn, costs, can be considerably reduced. In addition, the reduction of total incubation time to 24-h will increase the number of samples that can be analyzed for each cow. Additional measurements of the wash value and water solubility have to be carried out if the degradability is to be corrected for small particle loss. Still, there will be a substantial reduction of human work compared to the reference method.
The bilinear regression method may be considered a simple numerical integration of the degradation equation (Eq. [3]). Integrating this equation using a mathematical model for E(t) and assuming four stepwise constant D(t) values gives:
 | [7] |
where the integration limits, ti1, ti2, and ti3, are unknown. If we further assume that D(t4) is constant for all feeds, then the last term of Eq. [7] becomes a constant (a0). This assumption is of course crude, but if ti3 is high, then the contribution to ED from this term would be small. Writing the values of the integrals as parameters results in the following equation:
 | [8] |
Equation [7]
can be expressed as:
 | [9] |
This equation is the three-timepoint bilinear regression model where the integration limits in Eq. [8]
can be calculated from the estimated parameters a1, a2, and a3. Figure 5
shows two residence time distributions with passage rates of 0.04 h and 0.08 h-1, and three integration limits. The areas under the graphs from 0 to ti1 and from ti1 to ti2 are larger for the RTD curve with the highest passage rate (Table 5
). For this curve, the area under the curve from ti3 to infinity is also very small; therefore, we do not make an error by assuming a constant value for the last term in Eq. [7]. However, for the slowest RTD curve, this area is rather big, and the same assumption will cause much larger errors. Figure 2
shows that the error increases as the passage rate decreases for models with 24 h as the longest incubation time.

View larger version (11K):
[in this window]
[in a new window]
|
Figure 5. The residence-time distributions for passage rates of 0.04 h-1 () and 0.08 h-1 (). Three integration limits are indicated (ti1, ti2, and ti3).
|
|
Some simplified in situ techniques have been used for degradability measurements. Wilkerson et al. (1995) used only a 16-h ruminal incubation time to estimate the ruminal protein degradability in roughages by assuming that the 16-h sample directly estimated the escape protein. Broderick (1994) proposed a model that was based on estimating the degradation rate from the wash sample (0 h) and one sample incubated in the rumen (16 h) and by using acid detergent insoluble protein as an estimate of the total undegradable fraction. Vanzant et al. (1996) used this simplified technique to calculate the protein degradability in prairie hay and alfalfa. Using only the 0- and 16-h points to construct degradation rates resulted in similar estimates of protein degradation when compared with full time-series calculations. Calsamiglia et al. (1994) observed that degradation rate and degradability, based on double-point incubations (0 and 24 h or 2 and 24 h), explained at least 99% of the variation in measurements determined from curves based on seven timepoints in the range of 0 to 24 h. Nevertheless, in the present article, ED is calculated directly from the experimental disappearance values, which is in contrast to estimated degradation parameters, and where these are used in further calculations of degradability. Degradation parameters can be misleading when estimated from too few experimental values. Using two or three incubation times, one wrong disappearance value can result in an entire degradation curve being wrong. The approach used in the present paper is not that vulnerable for one wrong disappearance value since each value only contributes to a fraction of the calculated ED values. This makes the bilinear model approach more robust than the methods described by Broderick (1994) and Calsamiglia et al. (1994).
In this study, we have only tested models based on the discrete incubation times already used in the standard in situ method (0, 2, 4, 8, 16, 24, and 48 h). A more rigorous approach would be an optimization based on all possible incubation times. However, although somewhat limited in its approach, this study clearly demonstrated the potential of reducing the number of incubation times.
 |
Implications
|
|---|
The objective of this study was to test accuracy of in situ effective degradability measures in concentrates when using reduced numbers of ruminal incubation times. The simplified method was based on bilinear regression models, where one to three disappearance values were used to estimate effective degradability directly. Bilinear regression models, based on two and three incubation times, gave similar effective degradability estimates as calculated from seven or eight incubation times. Replacing the full-timescale nylon bag procedure with a two or three ruminal incubation times model may reduce the amount of human work considerably, and this reduction may be obtained without much loss in the accuracy of the method relative to the standard in situ method. However, this model is only valid for estimation of effective degradability or escape values and not for estimation of degradation rates.
 |
Footnotes
|
|---|
1 This study was financially supported by the Res. Council of Norway, SINTEF Group, and Felleskjøpet fôrutvikling. 
Received for publication December 6, 2001.
Accepted for publication September 13, 2002.
 |
Literature Cited
|
|---|
Broderick, G. A. 1994. Quantifying forage protein quality. Page 200 in Forage Quality, Evaluation, and Utilization. G. C. Fahey Jr., ed. Am. Soc. Agron., Madison, WI.
Calsamiglia, S., I. K. Yoon, and M. D. Stern. 1994. Effect of various incubation times on in situ estimation of ruminal crude protein degradation. J. Anim. Sci. 72(Suppl. 1):171.
Dhanoa, M. S., J. France, S. López, J. Dijkstra, S. J. Lister, D. R. Davies, and A. Bannik. 1999. Correcting the calculation of extent of degradation to account for particulate matter loss at zero time when applying the polyester bag method. J. Anim. Sci. 77:33853391.[Abstract/Free Full Text]
Fogler, H. S., 1992. Pages 738741 in Elements of Chemical Reactor Engineering. Prentice-Hall Int., Inc., London.
Hvelplund, T., and M. R. Weisbjerg. 1998. In vitro techniques to replace in vivo methods for estimating amino acid supply. Page 131 in In Vitro Techniques for Measuring Nutrient Supply to Ruminants. E. R. Deaville, E. Owen, A. T. Adesogan, C. Rymer, J. A. Huntington, and T. L. J. Lawrence, ed. Occasional Pub. No. 22, Br. Soc. of Anim. Sci. Proc., July 1997, Univ. of Reading, U.K.
López, S., J. France, J. Dijkstra, and M. S. France. 1999. Comparison of mathematical models to describe disappearance curves obtained using the polyester bag technique for incubating feeds in the rumen. J. Anim. Sci. 77:18751888.[Abstract/Free Full Text]
Madsen, J., 1985. The basis for the proposed Nordic protein evaluation system for ruminants. The AAT-PBV system. Acta Agric. Scand. 25(Suppl.):920.
Madsen, J., T. Hvelplund, M. R. Weisbjerg, J. Bertilsson, I. Olsson, R. Sprndly, O. M. Harstad, H. Volden, M. Tuori, T. Varvikko, P. Huhtanen, and B. L. Olafsson. 1995. The AAT/PBV protein evaluating system for ruminants. A revision. Norw. J. Agric. Sci. 19(Suppl.):137.
Poncet, C., B. Michalet-Doreau, T. McAllister, and D. Rémond. 1995. Dietary compounds escaping rumen digestion. Page 167 in Recent Development in the Nutrition of Herbivores. M. Journet, E. Grenet, M.-H. Farce, M. Thériez, and C. Demarquilly, ed. Proc. IVth Int. Symp. Nutr. Herbivores, INRA Editions, Paris.
Stern, M. D., A. Bach, and S. Calsamiglia. 1997. Alternative techniques for measuring nutrient digestion in ruminants. J. Anim. Sci. 75:22562276.[Abstract/Free Full Text]
Tamminga, S., W. M. Van Straalen, A. P. J. Subnel, R. G. M. Meijer, A. Steg, C. J. G. Wever, and M. C. Blok. 1994. The Dutch protein evaluation system: The DVE/OEB-system. Livest. Prod. Sci. 40:139155.
Vanzant, E. S., R. S. Cochran, E. C. Titgemeyer, S. D. Stafford, K. C. Olson, D. E. Johnson, and G. St. Jean. 1996. In vivo and in situ measurements of forage protein degradation in beef cattle. J. Anim. Sci. 74:27732784.[Abstract]
Vérité, R., and J. L. Peyraud. 1989. Protein: The PDI system. Page 33 in Ruminant Nutrition: Recommended Allowances and Feed Tables, R. Jarrige, ed. John Libbey Eurotext, Paris.
Volden, H., and O. M. Harstad. 1995. Effect of rumen incubation on the true indigestibility of feed protein in the digestive tract determined by nylon bag techniques. Acta Agric. Sci. Scand. 45:106115.
Wilkerson, V. A., T. J. Klopfenstein, and W. W. Stroup. 1995. A collaborative study of in situ forage protein degradation. J. Anim. Sci. 73:583588.[Abstract]