|
|
||||||||
Institute of Animal Breeding and Genetics, University of Goettingen, Albrecht-Thaer-Weg 3, D-37075 Goettingen, Germany
1 Corresponding author: ariek{at}gwdg.de
| ABSTRACT |
|---|
|
|
|---|
Key Words: llama colostrum milk composition lactation
| INTRODUCTION |
|---|
|
|
|---|
During lactation, milk composition undergoes specific changes. In many ungulates rising fat and protein concentrations are accompanied by increasing DM and declining sugar levels (Oftedal, 1984). The purpose of this study was to examine llama milk composition under a controlled feeding regimen during the course of a 27-wk lactation period and to describe the courses of milk constituents by suitable regression models.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
The samples were filled in test tubes (50 mL) containing 0.05% Bronopol for preservation. The samples were sent immediately in a cooler box at 2°C to the MKU (Milk Control and Research Federation, Uelzen, Germany) and analyzed for the constituents. The maximum time elapsed between sample taking and analyzing was 20 h.
Fat, protein, lactose, fat-free DM (FFDM; %), and milk urea concentration (mg/L) were determined by infrared absorption using an infrared spectral-photometer (Milkoscan FT 6500, Foss Electric, Hillerød, Denmark) following the guidelines of good laboratory practice from Lower Saxony (GLP, 2004). For calibration of the infrared spectral-photometer, cow milk was used. Correction factors were calculated by analyzing 10 llama milk samples analytically for fat (Röse-Gottlieb method as outlined in Marshall, 1992), protein (Kjeldahl method as outlined in Marshall, 1992), lactose (enzymatic method; Kleyn, 1985), and urea (Urea test-kit Reflotron, Boehringer Mannheim GmbH, Mannheim, Germany). Somatic cell count was determined using the flow cytometry method with an automated somatic cell counter (Fossomatic 5000, Foss Electric), and pH was measured using a pH meter (inoLab, WTW GmbH, Weilheim, Germany). Gross energy (GE) was estimated using the equation after Perrin (1958): GE (MJ/100 g) = 39.8 (fat %) + 23.9 (protein %) + 16.7 (lactose %). The DM concentration of milk was calculated by adding the fat concentration of the milk to the FFDM, and the ash concentration was calculated by subtracting protein and lactose from the FFDM concentration. Water concentration of milk was estimated calculating 100 DM. Milk urea was converted to MUN (mg/dL) to allow comparison with data from the literature.
For comparing values of milk constituents among species, the present data and results for major milk constituents from the literature were used to calculate the GE concentration in milk applying the equation of Perrin (1958), and values were expressed on a whole milk or GE basis.
Statistical Analyses
All statistical analyses were performed using SAS, version 8.01 (SAS Institute, 1999). There were no significant differences between the 2 trials, and data from both trials were pooled. Variations in the milk composition during the course of lactation were analyzed by the MIXED procedure using a mixed linear model. Because of the repeated measurement, female identity was included as a random effect, and parity, week of lactation, and the trial number were fitted as fixed effects. The SCC were transformed to logarithms to achieve normal distribution (logSCC).
In a preevaluation, 21 different prediction models were fitted to the weekly LSM of constituents by a nonlinear regression procedure (NLIN procedure). On the basis of best fit the following 2 models were chosen that were originally developed to describe lactation curves in cattle:
![]() |
where Yw is the constituent at week w, w is the week pp, e is the Eularian number, and a, b, c, d, and e are the parameters to be estimated.
The goodness of fit (R2) and the error mean square (EMS) were used for assessing model adequacy.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
Protein decreased in the first month from 4.98 ± 0.11% to 3.80 ± 0.12% but then steadily increased until the end of the trials to 4.42 ± 0.12%. On a GE basis, protein decreased in the first month of lactation from 31.7 ± 0.8% to 26.7 ± 0.9% and then remained constant around the average of 26.3%. The fat to protein ratio sharply increased within the first 4 wk from 80.4 ± 7.0% to 106.0 ± 7.3% and then reached 111.8 ± 6.9% until wk 27 pp.
Contrary to fat and protein, lactose increased from 5.78 ± 0.07% to 6.05 ± 0.08% in the first month of lactation and then slowly decreased to 5.57 ± 0.08% (wk 27 pp). On a GE basis, lactose showed a similar course as on a whole-milk basis and averaged 25.7 ± 0.8%. Water concentration followed a similar course as lactose and had slightly increasing values in the first 4 wk of lactation (84.4 ± 0.3% to 85.5 ± 0.3%) and decreasing ones for the rest of the lactation period (84.3 ± 0.4% in wk 27 pp). Complementarily, DM showed an opposite course and had decreasing values in the first 4 wk pp (15.61 ± 0.32% to 14.45 ± 0.36%) and increasing values for the rest of the lactation (15.68 ± 0.31% in wk 27 pp).
Gross energy concentration of the milk increased with the stage of lactation from 3.47 ± 0.13 kJ/g in wk 4 to 4.02 ± 0.11 kJ/g in wk 27 pp. The ash concentration had little variation.
The MUN reached the lowest value in wk 5 pp (16.2 ± 1.7 mg/dL) and the highest in wk 26 pp (26.7 ± 1.8 mg/dL). The 2 first-lactation animals had lower (P < 0.05) values than animals with more than 1 lactation (19.4 ± 1.4 vs. 22.6 ± 0.6 mg/dL). The mean weekly pH gradually decreased during the lactation and reached its lowest value of 6.60 ± 0.05 at wk 27 pp. Average SCC was 37 ± 2 x 1,000/mL, which is surprisingly low, even though SCC was highly variable (Table 2
). The highest weekly average value was observed in the first week with 70 ± 10 cells x 1,000/mL and the lowest in wk 25 pp with 31 ± 9 cells x 1,000/mL. Results were similar for transformed data (logSCC).
Prediction Curves
The 2 prediction models for the individual constituents shown in Figures 1
, 2
, and 3
fitted significantly to the data (P < 0.01). The estimated parameters for the 2 models and the corresponding R2 and EMS are given in Table 3
. Both regression models resulted in identical curves for fat, ash, and pH. For all other traits, the G model with more inflection points than the W model was superior in regard to R2 and EMS. However, the improvement in terms of R2 was negligible, except for protein and DM.
|
| DISCUSSION |
|---|
|
|
|---|
The high protein concentrations in llama milk from 4 to 12 h after parturition as found in the present study are consistent with values reported by Johnson (1994). Similarly, protein values reported for the Bactrian camel and the dromedary are close to 15% (Abu-Lehia, 1989; Zhang et al., 2005). Compared with domestic ruminants, protein concentration in llama colostrum is not markedly different (Hadjipanayiotou, 1995).
The low fat concentration found in the present study compares with reports for the Bactrian camel (Zhang et al., 2005) and the dromedary (Abu-Lehia, 1989). The higher fat concentrations found for ruminants are assumed to serve as a source of energy for the newborn calf (Merin et al., 2001). In contrast to cattle, suckling in llama crias takes place nearly every hour (Poullion, 2001) so that no long-term storage of nutritional energy is necessary.
Gross Composition of Milk
For ruminants it is well established that fat is the milk constituent affected most by influences such as feeding, lactational stage, health, etc. (Oftedal, 1984). High variations in the milk fat concentration were also reported for llamas (Morin et al., 1995), but the average value of 4.7% found in the present study is higher than that of 2.7% obtained by the former authors. However, the samples analyzed by Morin et al. (1995) were taken at various lactational stages and from several different locations with feeding practices not reported. The value reported by Fernandez and Oliver (1988) for llamas is in agreement with the present results, but it is unclear at which lactational stages the samples were taken and which diets were fed. Johnson (1994) reported a fat value for llama milk of 5.6%, but the methodology used was not described. For the dromedary, average fat concentrations of 3.6% have been reported (Sawaya et al., 1984), whereas Bactrian camels seem to have somewhat higher average fat concentrations (Zhang et al., 2005). Compared with other domestic ruminants the average fat concentration in llama milk is higher than in goat or cow milk but lower than in sheep milk (Table 4
).
|
The rather high protein concentration in wk 1 (Figure 1
) is probably still attributable to the transition from colostral to normal milk. Similar values for the first week pp were reported by Fernandez and Oliver (1988) and Johnson (1994) for llamas, whereas the value given by Morin et al. (1995) is somewhat lower. The reported results for the alpaca and the vicuña fall within the range of the observed protein concentration (Table 4
). Compared with domestic ruminants, only sheep reach similar average protein concentrations, whereas cows and goats have lower values (Table 4
).
The observed increase in milk fat and protein concentrations during lactation after the first month of lactation is in agreement with observations for many wild and domestic ruminants (Oftedal, 1984). Contrary to expectation, dromedaries and Bactrian camels are reported to have decreasing fat and protein concentrations over the lactation period (Merin et al., 2001; Zhang et al., 2005).
The development of the fat to protein ratio (Figure 1
) revealed that the sudden increase in the first 4 wk pp is attributable to the sharp decline of protein after transition to normal milk. Similar trends were observed in wild and domestic ruminants (Oftedal, 1984).
The high values for lactose found in the present study are consistent with other findings for llamas (Fernandez and Oliver, 1988; Morin et al., 1995) and their close relatives the alpaca (Medina and Bustinza, 1985) and the vicuña (Fernandez et al., 1997). Only Johnson (1994) reported a rather low lactose value of 3.30% for llamas, but the methodology applied was not described. Because lactose is the constituent most important to the secretion of the aqueous phase of the milk (Peaker, 1978), milk water tends to increase with increasing lactose values. Similarly, DM concentration decreases with increasing lactose values (Figure 1
). Those findings are in agreement with the general observation that species with high milk sugar concentration; for example, in the Perissodactyla, have low DM concentrations, whereas species with very low sugar concentrations such as some carnivores or lagomorphs have rather high DM concentrations that may surpass 50%, as in some seals (Oftedal, 1984).
The average DM concentration found is close to the value of 15.6% reported for llamas by Fernandez and Oliver (1988). Morin et al. (1995) observed a somewhat lower value of 13.1%. Similar results of 16.3% were also reported for the guanaco (Clavel et al., 2003).
The average ash concentration was in the range of the already published data on llamas and guanacos (Fernandez and Oliver, 1988; Johnson, 1994; Clavel et al., 2003) as well as for camels (Sawaya et al., 1984; Merin et al., 2001; Zhang et al., 2005) and domestic ruminants (Oftedal, 1984).
By calculating the GE (Table 4
) concentration from published milk composition data for llamas using the equation after Perrin (1958), the value found in the present study is in close agreement with Fernandez and Oliver (1988) and Johnson (1998). Only Morin et al. (1995) had a rather low GE value of 2.9 kJ/g, which can be attributed to the low fat concentration found in their study. The average GE concentration in llamas is similar to the milk energy of the other camelids, except the dromedary, which has lower energy concentrations (Table 4
). Because fat is the constituent in the milk with the highest caloric value (39.8 MJ/100 g), it is obvious that animals with high fat concentrations have a high GE concentration in their milk (e.g., sheep), whereas animals with low milk fat concentrations such as horses show low GE concentrations (Table 4
).
Considering the 3 energy contributors, fat, protein, and lactose, on a GE basis, it is interesting to note that fat contributed on average nearly half (48%) of the GE in the llama milk, whereas protein and lactose comprised the other half almost equally. Compared with domestic ruminants on a GE basis, llama milk seems to be nearly identical to cow milk and close to goat milk, whereas in sheep milk, fat contributes over 60% of the energy. Those findings should be considered when recommendations are given for milk replacers for llamas.
MUN, Milk pH, and SCC of Milk
Urea is the primary form of excretory N in mammals, and BUN is used to evaluate the efficiency of utilization of dietary CP by ruminants (Lewis, 1957). Because urea equilibrates rapidly throughout body fluids, including milk, MUN reflects BUN (Broderick and Clayton, 1997). To the authors knowledge, no data on MUN have been published for New World camelids, whereas publications on BUN values for llamas range from 13 to 32 mg/dL (Garry et al., 1994). Johnson (1994) reported BUN values of 19.3 and 29.6 mg/dL for llamas consuming diets containing 10 and 16% CP, respectively. The average weekly MUN values of the present study lie within the range of the published BUN data for llamas. Compared with optimum MUN levels from 7 to 14 mg/ dL for dairy cattle (Kirchgessner et al., 1986), the reference BUN for llamas or the MUN values in the present study are considerably higher. The observed high MUN values of 27 mg/dL at the end of lactation would suggest that the animals were then in a slight protein or energy surplus, or both, resulting from a decreasing milk production and a lower CP utilization. This finding could be explained by the fact that animals received 1 kg of concentrate, containing 16% CP, daily throughout the lactation regardless of the lactational stage.
Mean milk pH varies among species from 6.2 to 7.3 (Anderson, 1992). The mean milk pH in the llama (Table 2
) seems to be similar to the milk pH in goats (Baldi et al., 2002), although Rowan et al. (1996) reported a somewhat higher value of 6.93 for llamas. The decline of milk pH during lactation in the present study was also observed in the study of Anderson (1992). One explanation for the nearly linear decline may be a reduction of the ability of the mammary gland to produce bicarbonate. Anderson (1992) suggested that milk proteins could be responsible for the pH change in milk because, like blood proteins, they have the ability to control pH as a result of their buffering capacity. But because proteins are more anionic than cationic, an increase in protein concentration should result in an increase in pH. This was not observed in the present study. Another possible explanation is that an increase of organic acids such as citric, acetic, pyruvic, and lactic occurred in the milk.
The observed SCC values in the present study are within the range of the results of an extensive study conducted by Rowan et al. (1996) on udder health in llamas in North America. Compared with animals producing milk for human consumption (e.g., cows, goats, and sheep), llamas have a considerably lower SCC. The major determinant of the SCC in milk is the infection status of the gland (Rowan et al., 1996) and, as observed in the present study, the stage of lactation. Because there were no cases of clinical mastitis in both trials and also no significant effect of parity on SCC, the sole cause for the change in SCC in this study is the stage of lactation.
Prediction Curves
The curves obtained by the 2 nonlinear regression models were similar in shape for protein, GE, and DM. Comparable shapes of curves for protein percentage were also found for cattle (Stanton et al., 1992). The shape of the curves for lactose (Figure 1
) and water (Figure 2
) percentages are similar to the shape of lactation yield curves, whereas lactose had the best fit aside from protein. The prediction curves for MUN had the lowest fit of all constituents, except for ash. The reason for that is the high variation between the weekly means, with some outliers especially at wk 5, 8, 11, and 16 pp (Figure 3
).
| CONCLUSIONS |
|---|
|
|
|---|
The study showed that milk constituents in llama milk changed considerably over the lactation and were not affected by parity except for MUN.
Both regression models (W and G) applied, although developed to describe cattle milk yield lactation, were suitable to describe the course of the constituents of llama milk over the entire lactation. The differences between the 2 models in terms of model adequacy (R2) were small, so that for simplicity, the W model, although with less inflection points than the G model, seems adequate to describe the course of the constituents. These estimates may serve as a useful reference to establish standard values for the formulation of milk replacers at different stages of lactation for llama crias whose dams died or have agalactia under European housing conditions with low to medium feeding intensity. But because information on llama milk constituents is limited, especially over the entire lactation period, more systematic studies under different climatic and feeding conditions are needed.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication February 9, 2006. Accepted for publication April 20, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. S. Santos and A. M. Silvestre A Study of Lusitano Mare Lactation Curve with Wood's Model J Dairy Sci, February 1, 2008; 91(2): 760 - 766. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Riek and M. Gerken Measurements of the bodyweight and other physical characteristics of 11 llamas (Lama glama) from birth to weaning Vet Rec., October 13, 2007; 161(15): 520 - 523. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Riek, M. Gerken, and E. Moors Measurement of Milk Intake in Suckling Llamas (Lama glama) Using Deuterium Oxide Dilution J Dairy Sci, February 1, 2007; 90(2): 867 - 875. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |