|
|
||||||||
Department of Animal Sciences, The Ohio State University, Columbus 43210
2 Corresponding author: eastridge.1{at}osu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: alfalfa kinetics of biohydrogenation fatty acid multicompartmental model
| INTRODUCTION |
|---|
|
|
|---|
Grazing cows produce milk with higher concentration of CLA than those fed preserved forages (Jahreis et al., 1997; Kelly et al., 1998; French et al., 2000). Kay et al. (2004) reported that the majority of CLA in milk from cows fed pasture was derived from desaturation of vaccenic acid (trans-11 18:1; VA) by stearoyl-CoA desaturase (EC 1.14.19.1) in the mammary gland. Whereas desaturation of VA in the mammary gland contributes the majority of CLA in milk fat, less is understood concerning factors that regulate the amount of VA produced in the rumen. The higher concentration of CLA in milk fat from grazing cows has been proposed to be a result of higher levels of soluble sugars in fresh plants compared with preserved forages (Kelly et al., 1998; French et al., 2000). Concentration of sugars in grasses average 17% of DM; fructosans occur in major amounts, followed by sucrose, fructose, and glucose (Woolford, 1984). However, ensiling and drying decrease the amount of FA (Chilliard et al., 2000) and sugars in forages (Van Soest, 1994).
Kolver and de Veth (2002) reported ruminal pH to be lower than 6.2 for cows fed pasture-based diets. A lower pH may inhibit the final step of BH, favoring formation of monoenoic trans 18:1 FA (Qiu et al., 2004). Therefore, interactions among pH and concentrations of 18:3 and sucrose may play a role in the overall BH process (Ribeiro et al., 2005).
We are not aware of any published data reporting the BH rates of individual trans 18:1 FA and CLA isomers. Understanding how ruminal BH is regulated, even in in vitro conditions, may further advance our ability to develop methods for estimating the profile of FA reaching the duodenum of ruminants. A multicompartmental model approach may allow us to characterize the influence of dietary and ruminal factors on the BH rates of FA. Because all pathways of ruminal BH have not yet been defined, such an approach may also allow for improvement in understanding the BH pathways proposed in the literature by identifying segments of the model for which more data are needed, and therefore, stimulate further research in this area.
We postulated that the BH rate of UFA in fresh alfalfa is greater than that in alfalfa hay and that decreasing pH will decrease the BH rate independent of the forage preservation method. We also postulated that an interaction between pH and sucrose concentration also influences BH of 18:3 from alfalfa hay. To quantify the effect of pH on BH rates, we generated a multicompartmental model to estimate the BH rates, FA pool size over time, and flux of all FA during in vitro incubation of forages.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Treatments were similar to those from a previous experiment in which the effect of sucrose was also tested using continuous culture fermenters (Ribeiro et al., 2005); the treatments were 1) fresh alfalfa, 2) alfalfa hay, 3) alfalfa hay plus 4% sucrose, and 4) alfalfa hay plus 8% sucrose. The concentrations of nutrients from each forage source are shown in Table 1
.
|
Ruminal fluid was obtained from 2 ruminally cannulated dairy cows per set of incubation, one consuming an alfalfa hay diet and the other consuming a TMR diet consisting of 60% (wt/wt) forage (alfalfa hay plus corn silage). These were mixed to increase diversity of micro-organisms present during incubation. Ruminal contents were strained through 2 layers of cheesecloth to retain small particles, which are the reservoir for fibrolytic bacteria that contribute to BH. The ruminal fluids from 2 cows were mixed together and taken to the laboratory under anaerobic conditions at 39°C. The ruminal fluid was blended to dislodge particle-associated bacteria and strained through 8 layers of cheesecloth. We have previously observed that this procedure allows small particles to be present in the inoculum. The inoculum was divided into 2 equal volumes, and either SB or WB was added at 20% of total volume. Inocula were gassed with CO2 and blended while 30 mL was transferred to each of the centrifuge tubes containing the various substrates (Piwonka and Firkins, 1993, 1996).
Each set of 8 treatments was incubated for 0, 1, 2, 3, 4, 6, 9, and 12 h. Times were chosen to identify peak concentrations of CLA and VA in early hours and concomitantly with a lower standard error of means previously reported (Ribeiro and Eastridge, 2004). After incubation, pH was measured and the tubes were placed in ice to stop fermentation. The tube contents were transferred to aluminum pans, frozen, freeze-dried, and kept at 80°C until analyses were performed.
Analysis of FA
The FA of substrates and residues of incubation were methylated with 2 mL of 0.5 mol/L sodium methoxide (10 min at 50°C), followed by 3 mL of 5% methanolic HCl (10 min at 80°C) as described by Kramer et al. (1997). Methyl esters were separated by gas chromatography using a HP 5890 Series II gas chromatograph (Hewlett Packard Co., Palo Alto, CA). The column was a fused silica capillary (SP-2560, 100 m x 0.25 mm i.d. x 0.2 µm film thickness; Supelco, Inc., Bellefonte, PA). Helium was used as carrier gas. Detector and injector temperatures were 260°C, and the split ratio was 80:1. Oven temperature was 166°C for 39 min, increased by 10.0°C/min to 240°C, held for 10 min, increased by 3.0°C/min to 245°C, and held for 10 min. The temperature program used was optimized to separate most of the 18:1 FA in the first isothermal range as described by Molkentin and Precht (1995). Nonadecanoic acid was used as an internal standard. Retention times and response factors were determined with methyl ester standards purchased from Nu-Chek Prep (Elysian, MN; cat. no. GLC-60) and Matreya, Inc. (Pleasant Gap, PA; FIM-FAME-7). The 18:1 FA that were not available commercially (trans-6/8, trans-9, trans-12, trans-13, trans-15, cis-12, and cis-15) were identified by order of elution as shown by Molkentin and Precht (1995).
Statistical Analyses
Data were first analyzed as a complete randomized block design using PROC MIXED of SAS (SAS Institute, 2004) according to the following model:
![]() | [1] |
where Yijk = dependent variable for treatment i on block j and time k; µ = overall mean; Ai = fixed effect of treatment i; i = 1, 2, 3, 4, 5, 6, 7, 8; cj = random effect of block j; j = 1, 2; Tk = fixed effect of time k; k = 1, 2, 3, 4, 5, 6, 7, 8; ATik = interaction of treatment i with time k; and
ijk = residual error associated with the ijkth observation.
In model [1], there is no functional form imposed on the effect of time and time by treatment interaction.
Single Available Pool, First-Order Kinetic Model
The first-order kinetic model of Ørskov and McDonald (1979) was fitted with the PROC NLIN procedure of SAS (SAS Institute, 2004) using the least squares subclass means from [1] as observations. The model for estimating BH rates of 18:3 and 18:2 was as follows:
![]() | [2] |
where Y = amount (mg/tube) of FA at time t; C = pool of unavailable FA; P = pool of potentially available FA; k = fractional rate of FA disappearance (h1); t = incubation time (h); l = lag (h); and
i = residual error associated with the ith observation.
Model [2] was fit independently of each of the 8 treatments. The k and l parameter estimates for each treatment were compared by orthogonal contrasts to test the effects of buffer, level of sucrose (linear and quadratic), forage source, and interactions between buffer and sucrose, and buffer and forage source. Significance was declared at P < 0.05 and trends at P < 0.15.
Multiple Pool, First-Order Kinetic Model
The time effect can be modeled alternatively by using the more complex kinetic diagram with multiple pools. To do so, the least squares subclass means of model [1] were used as observations to estimate fractional rates of transfer among FA pools using SAAM II software (SAAM, 1997). Observations (least squares subclass means of [1]) were weighted by the reciprocals of their standard errors. The concentrations of the intermediates of BH from 0 to 12 h were used also in SAAM II to estimate transfer rates among the pools. The decision to specify direct transfer among specific pools was made using 2 sequential steps: 1) if there were published data showing a direct link for the pools, and 2) the addition of a pool or transfer produced a smaller value of the Akaike information criterion (AIC; Cobelli and Foster, 1998):
![]() |
where WRSS = weighted residual sum of squares, and P = number of unknown model parameters.
There is evidence that cis-9, trans-11 18:2 (one of the CLA isomers) is the first product of the BH of 18:2 (Kepler et al., 1966; Noble et al., 1974); therefore, the pool of 18:2 has a transfer rate to the CLA pool. The Rosenbrock integrator method (SAAM, 1997) was used and the optimization was performed using a variance model based on relative data and a forward derivative. Because there are no data on the transfer rate values for the intermediates of BH, the initial parameter estimates for the BH rates of 18:2 and 18:3 were the same as the ones generated by the equation [1]. Means of BH rates were compared using the paired difference test (Snedecor and Cochran, 1980).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Effects of pH, Sucrose, and Forage Source on BH
Because the WB treatment resulted in a lower pH than the SB, we considered that the main buffer effect was a result of changes in pH. Lower pH is associated with an impairment of lipolysis (Van Nevel and Demeyer, 1996), BH, and lower fiber digestibility (Qiu et al., 2004). The decrease in BH with lower pH may be associated with a decrease in the numbers of cellulolytic bacteria (Wales et al., 2004), because cellulolytic bacteria play a major role in ruminal BH and these are known to be decreased at lower pH (Harfoot and Hazlewood, 1997). A pH of 6.2 is considered to be the lower threshold for optimum fiber digestion (Grant and Weidner, 1992), although the period of time that pH is less than 6.2 may be more significant in inhibiting growth of cellulolytic bacteria (Calsamiglia et al., 2002). Moreover, low pH inhibits lipolysis (Van Nevel and Demeyer, 1996), decreasing BH by diminishing the availability of free UFA. Therefore, the lower BH rates for 18:2 and 18:3 with the WB treatment were expected. Qiu et al. (2004) observed reduced cellulolytic bacterial numbers with reduced BH of 18:2 and a 50% decrease in total BH when the pH changed from 6.5 to 5.8 in continuous culture. Kolver and de Veth (2002) reported that the ruminal pH of grazing cows may be as low as 5.8. Because low pH may increase the proportion of intermediates of BH (Qiu et al., 2004), the diurnal variation in pH of grazing ruminants and the high 18:3 concentration of pasture could increase the flow of VA out of the rumen. Others have shown that lower pH decreases lipolysis (Van Nevel and Demeyer, 1996) and the extent of BH, and increases the proportion of monoenoic FA in the medium (Van Nevel and Demeyer, 1996; Griinari et al., 1998; Troegeler-Meynadier et al., 2003). The BH rates of 18:3 was almost always higher than the rates of 18:2 and it was 54% higher when fresh alfalfa, compared with hay, was incubated with the WB. Thus, fresh alfalfa may be able to maintain the number of some cellulolytic (biohydrogenating) species even at decreased pH.
We did not observe any lag time during BH of 18:2 and 18:3. Zero lag times have been reported previously with incubation of forages (Boufaied et al., 2003) and oil seeds (Enjalbert et al., 2003). Conversely, Troegeler-Meynadier et al. (2003) reported lag times between 1 and 2 h for soybean oil. Lag time is related to the time needed to release the glyceride from its matrix and is dependent on fat source (Beam et al., 2000) and might be associated with DM digestibility and time for the microbes to adapt to the substrate and incubation conditions. The UFA of fresh plants are located mainly in the chloroplasts as mono- and digalactosyldiglycerides (Butler and Bailey, 1973; Hawke, 1973), and hydrolysis of plant galactolipids has been reported to be 78 to 95% complete after 4 h of incubation (Dawson and Hemington, 1974). In fact, rapid hydrolysis of galactoglycerides and triglycerides (Beam et al., 2000) may explain the lack of lag time observed. Additionally, unesterified FA at 0 h in the incubation tubes would have contributed to the differences in lag time observed among experiments.
Boufaied et al. (2003) reported higher rates of 18:3 BH for fresh grass compared with hay, supporting our findings. Those authors argued that, because dry hay may have a lower DM digestibility, the FA could have been trapped physically within the organic matrix, decreasing rate of release, and therefore, decreasing BH of 18:2 and 18:3. Because the conservation process decreases the absolute and relative amounts of UFA of fresh forages (Chilliard et al., 2000), the percentage of FA in fresh alfalfa was 1.6 times higher than that for alfalfa hay (Table 1
). Therefore, one expected outcome in nongrazing cows is reduced concentration and out-flow of 18:3 FA from the rumen.
Sucrose had no effect on the rates of BH of 18:2 and 18:3; this observation did not support our hypothesis that sucrose would inhibit BH. We postulated that sucrose incubated with WB would result in greater impairment of BH than when incubated with SB, but we observed no interaction (P > 0.05) between sucrose and buffer source on the BH rates of 18:2 and 18:3. Conversely, continuous input of sucrose has been shown to impair the BH of 18:2 and 18:3 from alfalfa hay in continuous culture by a pH-independent process (Ribeiro et al., 2005). Perhaps, the lack of an adaptation period under in vitro conditions as occurs with continuous culture may have caused the different responses in rates of BH between the 2 methods of incubation.
Intermediates of BH
Because the proportion of 18:2 in forages is low compared with traditional diets used in nongrazing conditions, the amount of CLA in the incubation tubes was too low to be measured accurately; therefore, CLA concentrations are not reported. The concentration of VA was highest with the fresh alfalfa treatment (Figure 2
), which was expected because there was 79% more 18:3 in fresh alfalfa than in alfalfa hay. In addition, the majority of the VA may originate from BH of 18:3 when forage is the only source of FA (Sudarshan and Hawke, 1979).
|
The trans-13/14 18:1 FA had the second highest concentration of all 18:1 intermediates. Loor et al. (2003) reported that VA was the primary intermediate of BH, followed by trans-13 18:1 FA during continuous culture fermentation of grasses. The procedure to identify FA used in this trial did not allow the separation of trans-13 from trans-14 FA. The synthesis of both FA from 18:3 BH has been reported previously (Ward et al., 1964; White et al., 1970); however, the exact pathway and microorganisms involved are not known.
Kinetics of BH
Most of the research with in vitro BH used nonlinear regression and simple first-order kinetics to estimate the BH rates of 18:2 and 18:3, whereas no studies have reported BH rates of VA and other intermediates of BH (Beam et al., 2000; Enjalbert et al., 2003; Troegeler-Meynadier et al., 2003). The factors that control the concentration of VA during ruminal BH are significant because VA plays a central role in BH (Harfoot and Hazlewood, 1997), and it is the major source of CLA in milk of grazing cows (Kay et al., 2004). Modeling ruminal BH is a useful tool to characterize the influence of diet, microbial ecology, and ruminal factors on the BH process, and to postulate potential BH pathways. However, there are limited data in the literature in which the BH rates of the intermediates of BH have been estimated. Therefore, we developed a multicompartmental kinetic model to estimate the fractional rates of VA appearance and disappearance to characterize the influence of dietary and ruminal factors on pool size and flux of VA during in vitro BH.
First-order kinetics has been used previously to estimate BH rates of UFA from different sources (Beam et al., 2000; Boufaied et al., 2003; Troegeler-Meynadier et al., 2003). A second-order kinetic model would be a much more difficult model to integrate and develop. We used first-order kinetics in our model because we expected that the UFA concentration would be below saturates because of the low percentage of FA in the alfalfa samples (2%) such that pseudo first-order kinetics would apply. On pools that we evaluated, the natural logarithm concentration vs. time resulted in a linear plot (first order), and the inverse of concentration vs. time was not linear (i.e., not second order). Nonetheless, future research should evaluate if a multicompartmental, first-order kinetic model would be ideal when other substrates are utilized.
The WB treatment slowed the BH of the 18:3-Int pool (Figure 3
), and after 4 h, the 18:3-Int pool size was greater for the WB treatment, consistent with the lower BH rate of 18:3 and may be consistent also with the later peak of VA observed for the WB (Table 2
and Figure 2
). The temporal changes in the CLA and 18:3-Int pool sizes are illustrated in Figure 5
. The 18:3-Int pool increased more sharply when fresh alfalfa was incubated with the SB than with WB. Also, the 18:3-Int pool decreased faster than the CLA pool, regardless of the buffer used, because the BH rates of 18:3 was faster than for 18:2. Additionally, the CLA and 18:3-Int pools had higher BH rates compared with the 18:1 isomers, which is consistent with previous reports that those FA do not accumulate during ruminal BH as much as 18:1 isomers (Noble et al., 1974; Body, 1976); they are also associated with the same group of biohydrogenating bacteria (Harfoot and Hazlewood, 1997), reflecting the faster BH rates. The variation of the CLA concentration over time followed the same pattern as reported previously (Noble et al., 1974; Kellens et al., 1986) and was similar to the VA curve, but peaked earlier, consistent with precursor and product kinetics.
|
|
The flux of appearance of VA was higher during the first 6 h of incubation (Figure 4
). This flux was estimated by summing the fluxes from the 18:3-Int, CLA, and cis-18:1 pools to the VA pool as illustrated in our model (Figure 3
). Although the transfer rates are similar between appearance and disappearance of VA, the higher flux of appearance resulted from the higher mass of the 18:3-Int pool during the early time points. The high standard error for the BH rate of CLA and 18:3-Int to VA (Figure 3
) is a reflection of their high correlation (r = 0.70; data not shown). Additionally, because we were unable to set a direct link between 18:2-Int and CLA pools (due to overparameterization of the model), the standard error may have also been inflated.
The last step of BH is assumed to be most sensitive to pH because trans-18:1 FA accumulate when high-grain diets are fed (Griinari et al., 1998). However, the WB decreased the rates of BH of 18:2 and 18:3 without affecting the mean concentration and BH rate of VA. A shift to alternate pathways of BH that do not include VA could explain the lack of response of VA concentration to lower pH. For instance, WB decreased (P < 0.05; Figure 3
) the BH of other 18:1 FA pools (trans-9, trans-10, trans-12, and trans-13) to 18:0 and the hydrogenation of the 18:3-Int pool to the other 18:1 FA pool. Consequently, the lower rate of disappearance of 18:2 and 18:3 caused a more pronounced decrease in the flux of the UFA to other trans FA while maintaining almost similar fluxes to VA.
Limitations of the Modeling Approach
Although using a multicompartmental model to estimate BH is superior to a single-compartmental approach to model the dynamics of FA metabolism in the rumen, there are constraints and assumptions that limit interpretation. One constraint is that all rate and pool estimates are correlated and any error in estimating one rate will affect others. For instance, the 18:0 pool size may be underestimated; because 18:0 is the final product of BH, any error in estimating the precedent pools will be reflected in the 18:0 pool. Moreover, any intermediate of BH and microbial FA synthesis not taken into account that contributes to the synthesis of 18:0 will increase the error of estimating the pool size of 18:0.
Although some intermediates of BH of oleic acid have been reported (Ward et al., 1964; Mosley et al., 2002), the pathways and microbial species for BH of oleic acid have not been defined. Therefore, defining direct links between pools may be compromised. For example, we assumed that cis- and trans-15 18:1 were hydrogenated to 18:0; however, Harfoot and Hazlewood (1997) considered these to be true end products of 18:3 BH. However, there is still the possibility that cis- and trans-15 18:1 are in fact hydrogenated directly to 18:0, because we can not rule out this step based on the research done by Body (1976) and White et al. (1970). Including hydrogenation of cis- and trans-15 18:1 to 18:0 in the model improved the fit of the model (AIC value) because the concentration of both decreased after 9 h of incubation. Alternatively, these may have been isomerized to other 18:1 FA (Mosley et al., 2002; Proell et al., 2002) and then hydrogenated to 18:0; however, the model under this option did not solve.
Compartmental analysis of BH of FA in fresh alfalfa showed that the mean concentration and BH rate of VA were not affected by buffer, whereas the last step of BH of the WB treatment was slower for other trans-18:1 FA, and the mean concentration of trans-10 18:1 was higher. The concentration of trans- and cis-15 18:1 FA declined, most likely due to hydrogenation to 18:0, an observation not previously reported. Our data demonstrated the need to define even more the intermediates of BH (such as individual CLA isomers), as well the fractional rates of transfer among pools as dietary management is changed and microbial populations shift. Such knowledge would greatly improve understanding of the complex reactions that occur in the rumen and that affect the profile of FA in ruminant products.
| CONCLUSIONS |
|---|
|
|
|---|
We quantified changes in the rates of BH of trans-18:1 FA and VA, which is the key regulatory step determining the amount of VA leaving the rumen and affecting concentration of milk CLA. The multicompartmental model analysis allowed factors affecting the BH rates of 18:2 and 18:3 to be defined, as well as the synthesis and hydrogenation rates of other BH intermediates. We demonstrated the magnitude of BH pathways, other than through VA, that may explain variations in VA concentration as incubation conditions change. Because we could estimate fluxes, as well as mass of the VA pools, more information was generated from the data when using multiple pools compared with a single pool, first-order kinetic model. Future research should focus on increasing the number of individual FA pools that can be estimated by the model and quantify changes on BH rates with dietary and ruminal factors.
| FOOTNOTES |
|---|
Received for publication July 22, 2006. Accepted for publication October 27, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. L. Firkins, B. S. Oldick, J. Pantoja, C. Reveneau, L. E. Gilligan, and L. Carver Efficacy of Liquid Feeds Varying in Concentration and Composition of Fat, Nonprotein Nitrogen, and Nonfiber Carbohydrates for Lactating Dairy Cows J Dairy Sci, May 1, 2008; 91(5): 1969 - 1984. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Moate, R. C. Boston, T. C. Jenkins, and I. J. Lean Kinetics of Ruminal Lipolysis of Triacylglycerol and Biohydrogenation of Long-Chain Fatty Acids: New Insights from Old Data J Dairy Sci, February 1, 2008; 91(2): 731 - 742. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. C. Jenkins, R. J. Wallace, P. J. Moate, and E. E. Mosley BOARD-INVITED REVIEW: Recent advances in biohydrogenation of unsaturated fatty acids within the rumen microbial ecosystem J Anim Sci, February 1, 2008; 86(2): 397 - 412. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |