Genetic parameter estimation for milk [beta]-hydroxybutyrate and acetone in early lactation and its association with fat to protein ratio and energy balance in Korean Holstein cattle.
Subclinical ketosis in dairy cattle can be defined as the presence of elevated levels of circulating ketone bodies in blood, milk, and urine, without expression of typical clinical signs. It has been suggested that early lactation energy status and test-day fat and protein percentages should be considered as the potential indicators for the monitoring of sub clinical ketosis risk (SCK) [1,2].
Milk fat protein ratio (FPR) is a potent, easily quantifiable trait which can be utilized to differentiate cows with higher adoptability to early lactation challenges. The early lactation energy deficit leads to lipolysis and body fat mobilization resulting in an increased fat synthesis in the udder. Simultaneously, inadequate intake of fermentable carbohydrates can cause an insufficient protein synthesis by ruminal bacteria and the flow of amino acids to the udder is compromised and milk protein content decreases. Both of these processes result in an increased FPR.
An FPR >1.5 indicates abnormally high lipolysis and has proven to be a good predictor of SCK .
Most dairy cows experience an negative energy balance (NEB) in early lactation, because their dietary intake cannot support the required energy, due to metabolic priority for milk production. The continuous mobilization of body fat reserves to balance the deficit between energy intake and energy production leads to ketogenesis, which ultimately exhausts the cow . So excess level of ketone bodies are excreted to blood, which can be detected even in milk and urine.
Serum ketone body analysis is the major reference test to diagnose hyperketonemia. However due to practical limitations and close relationship between blood and milk ketone body concentrations, milk composition records become a convenient approach for routine analysis. The concentration of milk acetone (MAc) is similar to that of blood acetone (BAc), with mean MAc/BAc ratio of 1.04 . Furthermore, a low correlation coefficient (r = 0.45) between blood AcAc and milk AcAc have been recorded . The correlation between blood [beta]-hydroxybutyrate (BHBA) and milk BHBA ranged from 0.00 to 0.87 . The mean milk AcAc to BAcAc ratio and the mean milk BHBA to blood BHBA ratio were reported as 0.40 and 0.12, respectively . Accordingly the use of milk ketone body concentrations to discriminate ketotic cows proved to be almost as reliable as blood ketone body concentrations.
For the use and application of milk ketone bodies in effective genetic evaluation and improvement, the knowledge of genetic parameters for the economically important traits is essential. Dairy cows display individual variation in metabolic adaptation during NEB . This variation can be used for the selection of high producing cows with a lower hyperketonemia risk. Considerable genetic variation in heritability estimates of milk BHBA concentration ranging from 0.14 to 0.29 has been previously reported .
Therefore, this study was conducted to estimate the genetic parameters for milk BHBA, acetone (Ac), FPR, and energy balance (EB) and to further investigate the effect of early lactation FPR and EB on raising the milk ketone body levels in Korean Holsteins by utilizing milk test day records.
MATERIALS AND METHODS
Fourier transform infrared (FTIR) measurements for milk Ac and milk BHBA levels along with routine test day records were collected from April 2012 to August 2016 by the Korea Animal Improvement Association (KAIA). Test-day milk sampling and analysis were performed according to the Korean milk-recording procedures . The dairy population was managed under recommendations of KAIA. Herds are usually monitored for hyperketonemia by regular assessment of the percentage of fresh cows having ketone body concentrations above a defined threshold value. The test day milk records included milk yield, fat and protein percentages of cows that were 1 to 100 days in milk at sampling. The cows were milked twice daily at 0500 and 1600 h. Milk composition was analyzed weekly based on samples collected from 2 consecutive milkings. Season of calving was defined as summer (May to October) and winter (November to April). The age at first calving ranged from [less than or equal to]23 months, 24 to 25 months, 26 to 28 months, and [greater than or equal to]29, respectively. The morning and evening records of registered cows with pedigree information of test days from cows that were 5 to 305 days in milk (DIM) at sampling were included in the analysis. Samples with missing sire or dam information were excluded from the analysis. Outliers over [+ or -]3 standard deviations (SD) for milk production measurements were considered to be incorrect recordings and excluded from the analysis.
Test-day milk samples were analyzed by FTIR spectroscopy using a CombiFoss FT+ system (Foss Analytical A/S, Hilleroed, Denmark) with previously developed calibration equations for milk BHBA and milk Ac. The original data set consisted of 262,940 test-day records. Fat to protein ratio data were obtained from test day records.
This EB equation was based on parity, lactation week, and milk composition volume as follows [11,12]:
eEB = 217.8 - wk2x31.9 - wk3x20.6 - wk4x15.6 - wk5x11.5 - wk6x8.0 - wk7 8 10.6 - wk8x7.2 - wk9x5.3 - wk10x4.0 - wk11x2.7 - wk12x0 - par1x34.9 - par2x7.2 - par3x6.7 - par4x0 - milkx2.11 - protx15.36 - FP x49.24 (MJ nel/d)
Where; wk2, wk3 ... to wk12 = lactation wk 2 to 12; par1 - par4 = parity categories 1 to 4; milk = milk yield (kg/d); prot = % milk protein; FP = ratio of % fat to % protein in milk.
For descriptive statistics of all parameters SAS 9.2 package (SAS Institute Inc. Cary, NC, USA) was used. The test day records of milk BHBA, Ac, EB, and FPR were analyzed for each parity. Genetic parameters, including genetic (co)variance components, were estimated by the restricted maximum likelihood procedure based on repeatability model using the Wombat program . Heritability was analyzed within each parity and across all parities using two models.
The linear model for across all parities and for individual parities were illustrated as follows (Model 1 and Model 2);
Model 1: [Y.sub.ijklmn] = [mu] + DIM + [Age.sub.1] + [Season.sub.j] + [ampm.sub.k] + [parity.sub.1] + [a.sub.m] + [p.sub.m] + [e.sub.ijklm]
Model 2: [Y.sub.ijkl] = [mu] + DIM + [Age.sub.j] + [Season.sub.k] + [ampm.sub.l] + [a.sub.m] + [p.sub.m] + [e.sub.ijkm]
Where, [Y.sub.ijklmn] is BHBA, Ac, FPR, or EB observation; [mu] is the overall mean; DIM is the covariate describing the effect of days in milk, Aget is the fixed effect of calving age i, Season is the fixed effect of calving season j, [ampm.sub.k] is the fixed effect of milk collecting time k. The [parity.sub.l] is the fixed effect of parity l. The [a.sub.m] is the additive genetic effect of cow m, [p.sub.m] is the permanent environmental effect of cow m, e is the random residual effect associated with each record.
Heritability was calculated using this equation:
[h.sup.2] = [[sigma].sup.2.sub.A]/[[sigma].sup.2.sub.P]
Where [h.sup.2] is heritability, [[sigma].sup.2.sub.A] is additive genetic variance and [[sigma].sup.2.sub.P] is phenotypic variance.
A total of262,940 milk samples were evaluated for BHBA, Ac, protein, fat and milk yield at KAIA. The mean BHBA level for above data set was 43.48 [micro]mol/L with the range of 0 to 4,280 [micro]mol/L while the mean Ac level was 127.49 [micro]mol/L with the range of 0 to 3,570 [micro]mol/L. Average milk yield was 16.11 kg/d while mean fat and protein percentages were 3.27 and 3.84 respectively (Table 1).
This study evaluated the relationship between milk ketone bodies (BHBA and Ac) with the indicator traits (FPR and EB). Mean milk BHBA and Ac level declined with lactation week. In the early lactation (up to 6th week) both traits declined at a higher rate. During that period, Ac shows higher slope than BHBA. But after that the mean BHBA and Ac concentration become static. The FPR ranged from maximum value of 1.30 and minimum 1.15. DIM were split into 7 day intervals. The mean FPR was >1.2 up to 5th week of lactation. Cows with a FPR [greater than or equal to] 1.2 had greater milk BHBA and Ac up to 6th week of lactation, then the milk BHBA and Ac level became static. Highest mean BHBA and AC concentrations were 73.78 [micro]mol/L and 158.05 [micro]mol/L while lowest level was 33.48 [micro]mol/L and 85.01 [micro]mol/L respectively (Figure 1).
Figure 2 shows the relationship between mean EB and ketone body concentrations. In the 2nd week of lactation EB was -13.18 MJ nel/d. The negativity was reduced with lactation week and became positive at 5th week of the lactation (35 DIM). In the 5th week of lactation EB was 1.81 MJ nel/d and it increased up to 15.1 MJ nel/d by the 12th week of lactation. During the phase of NEB, elevated levels of milk BHBA and milk Ac concentrations can be observed. Afterward mean EB becomes positive and ketone body concentrations stabilize and even out. In this study, the estimated heritability for milk BHBA and milk Ac for all parities was 0.19 and 0.28 respectively. Heritability estimates for EB and FPR were 0.50 and 0.46 respectively. Genetic and phenotypic relationships of milk BHBA, Ac, EB, and FPR were studied. Phenotypic correlations among the traits considered were lower than genetic correlations with the exception of BHBA and Ac. Genetic and phenotypic correlations between Milk BHBA and EB, Ac and FPR, EB and FPR were negative while other correlations were positive. The highest negative correlation was observed between EB and FPR (Genetic correlation of -0.86 and phenotypic correlation of -0.81) (Table 2).
The heritability estimates for milk BHBA, EB, and FPR declined with the parity while heritability for Ac deviated from that pattern. The highest heritability for milk BHBA (0.14) was observed at parity 1 and it declined to 0.09 at the 4th parity. The observed 0.52 heritability for EB at parity 1 declined to 0.19 at parity 4 and 0.43 heritability for FPR at parity 1 declined to 0.16 at parity 4. Milk Ac showed the lowest heritability at the 1st parity (0.31) then it ranged from 0.29 to 0.31 during the subsequent parities (Table 3, Figure 3).
Variance components (additive genetic variance and phenotypic variance) and heritability estimates for different parities are shown in Table 3. The proportion of phenotypic variance was larger than additive genetic variance.
There was an apparent trend of increasing heritability with increasing parity.
The energy metabolism of dairy cows along with the milk composition including ketone bodies was the concern of this study. Dairy herds are monitored for the prevalence of hyper ketonemia by regular assessment milk composition volumes above a defined threshold value.
Various thresholds of BHBA and Ac in milk have been used to define SCK. Consequently, the concentrations rather than presence or absence of hyperketonemia in cows were used to investigate genetic variations in ketosis susceptibility.
Data from the first week of lactation were excluded in estimating FPR and EB, because of possible effect of the colostrum component of milk. FPR was considered as an indicator of lipo-mobilization and NEB status of postpartum cows. Several studies have indicated different values for the optimal FPR, whereas 1.33 was set as a cut off value  while FPR values higher than 1.5 suggest a great energy deficiency and SCK if ketone bodies are present . Animals with FPR greater than 2.0 in early lactation showed an increase in postpartum diseases . A positive correlation of BHBA levels with the fat: protein ratio ([R.sup.2] = 0.42; p<0.001) was evident . Values recorded in this study suggest that FPR of 1.3 indicates the likelihood of elevated levels of ketone bodies with the probable risk of SCK.
Dairy cows exhibit variation in the adaptation to NEB with respect to individual and herd factors. This study shows, elevated levels of ketone body concentrations during the period of NEB (up to 4 week). It suggests that heavy fat mobilization and ketogenesis in this high metabolism period causes a reduction of energy reserves. Buttchereit et al  suggested that animals change to a state of positive EB around DIM 40 to 80, which coincides with findings of de Vries and Veerkamp .
The genetic parameter estimation focused on the challenge of using the existing genetic variation for genetic selection of high-performing dairy cattle with a lower hyperketonemia susceptibility. This study estimated heritabilities for milk BHBA, milk Ac, FPR, and EB in early lactation to assess the susceptibility to hyperketonemia and investigate the possibility of using test day records in milk for genetic improvements. The previous studies suggested that the prevalence of hyperketonemia varies across parity, stage of lactation and season in cows [1,8], Therefore these factors were included in the models and heritability was analyzed within each parity (j = 1, 2, 3, 4) individually and across all parities.
The heritability estimates were somewhat higher when compared with previous studies on genetic variation. So, we can suggest that a considerable part of the observed variation could be allocated to fixed and random effects included in the model. Residual variance could be low. But our heritability estimates are also within the range of heritabilities in several studies on clinical ketosis. The heritability estimates of a previous study for milk BHBA and milk Ac were 0.04 to 0.17 and 0.22 to 0.29 respectively . Koeck et al  found that heritabilities for milk BHBA at different stages of early lactation from 5 to 100 DIM were between 0.14 and 0.29. High genetic merit cows experience more severe NEB in early lactation . Heritability estimates of EB ranged from, 0.12 to 0.49 over monthly intervals . In this study high heritability was observed in EB (0.50). Thus, it suggests that a favorable metabolic adaptation to the NEB in early lactation could be used for genetic selection.
Generally, tests for fat and protein content of milk as a basis for regular dairy herd improvement have been considered as alternate measures for determining the early lactation nutritional status of dairy cows . According to the breed and herd factors, milk fat and protein levels can vary tremendously. Holsteins have the lowest fat and protein content when compared with other breeds (USDA- AIPL summary of herds DHI test during 2004). Heritability estimate for FPR in this study was 0.46. Such a high heritability estimate indicates inheritance accounts for considerable difference between cows in protein and fat content in milk. German Holstein cattle showed e heritability estimates of 0.06 and 0.30 for EB and FPR respectively .
The highest negative phenotypic correlation -0.81 was observed between EB and FPR while highest negative genetic correlation -0.86 also observed between EB and FPR. So, variation between cows exists for fat protein mobilization during the metabolic adaptation to NEB in early lactation, which increases the hyperketonemia risk.
Animals with higher FPR have less body reserves and NEB. Variation may exist in cows for dietary intake, fat mobilization, metabolic gene expression which all influence the metabolic adaptation to NEB in early lactation.
Milk BHBA, Ac, EB, and FPR can be used as indicator traits for hyperketonimia. Milk test day records, provide practical routine data collection for breeding programs. FTIR predictions for milk BHBA and Ac can be used for screening cows for hyperketonemia. Milk BHBA, Ac, EB, and FPR showed a considerable genetic variation with heritability along with parity. Milk test day records provide useful data for selective breeding programs. Selective breeding may contribute to maintaining the milk ketone bodies at optimum level during early lactation.
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
This work was conducted with the support of the Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01268002) Rural Development Administration, Jeonju-si, Jeollabuk-do, Republic of Korea.
[1.] Duffield TF, Kelton DF, Leslie KE, Lissemore KD, Lumsden JH. Use of test day milk fat and milk protein to detect subclinical ketosis in dairy cattle in Ontario. Can Vet J 1997;38:713-8.
[2.] Toni F, Vincenti L, Grigoletto L, Ricci A, Schukken YH. Early lactation ratio of fat and protein percentage in milk is associated with health, milk production, and survival. J Dairy Sci 2011; 94:1772-83.
[3.] Cejna V, Chladek G. The importance of monitoring changes in milk fat to milk protein ratio in Holstein cows during lactation. J Cent Eur Agric 2006;6:539-46.
[4.] Collard BL, Boettcher PJ, Dekkers JC, Petitclerc D, Schaeffer LR. Relationships between energy balance and health traits of dairy cattle in early lactation. J Dairy Sci 2000;83:2683-90.
[5.] Enjalbert F, Nicot MC, Bayourthe C, Moncoulon R. Ketone bodies in milk and blood of dairy cows: relationship between concentrations and utilization for detection of subclinical ketosis. J Dairy Sci 2001;84:583-9.
[6.] Andersson L. Concentrations ofblood and milk ketone bodies, blood isopropanol and plasma glucose in dairy cows in relation to the degree of hyperketonaemia and clinical signs. Transbound Emerg Dis 1984;31:683-93.
[7.] Geishauser T, Leslie K, Kelton D, Duffield T. Evaluation of five cowside tests for use with milk to detect subclinical ketosis in dairy cows. J Dairy Sci 1998;81:438-43.
[8.] Van der Drift SG, Van Hulzen KJ, Teweldemedhn TG, et al. Genetic and nongenetic variation in plasma and milk beta-hydroxybutyrate and milk acetone concentrations of earlylactation dairy cows. J Dairy Sci 2012;95:6781-7.
[9.] Koeck A, Jamrozik J, Schenkel FS, et al. Genetic analysis of milk beta-hydroxybutyrate and its association with fat-to-protein ratio, body condition score, clinical ketosis, and displaced abomasum in early first lactation of Canadian Holsteins. J Dairy Sci 2014;97:7286-92.
[10.] Cho KH, Park B, Choi J, et al. Development of international genetic evaluation models for dairy cattle. J Anim Sci Technol 2013;55:1-6.
[11.] Heuer C. The use of test day information to predict energy intake of dairy cows in early lactation. J Dairy Sci 2004;87:593-601.
[12.] Heuer C, Van Straalen W, Schukken Y, Dirkzwager A, Noordhuizen J. Prediction of energy balance in a high yielding dairy herd in early lactation: model development and precision. Livest Prod Sci 2000;65:91-105.
[13.] Meyer K, WOMBAT-A program for mixed model analyses by restricted maximum likelihood. User notes. Version 1.0. Animal Genetic and Breeding Unit. Armidale, Australia: University of New England; 2010.
[14.] Garcia C, Montiel R, Borderas T, Girard V Relationship between [beta]-hydroxybutyrate and the fat: protein ratio of milk during early lactation in dairy cows. Arch Med Vet 2015;47:21-5.
[15.] Buttchereit N, Stamer E, Junge W, Thaller G. Short communication: Genetic relationships among daily energy balance, feed intake, body condition score, and fat to protein ratio of milk in dairy cows. J Dairy Sci 2011;94:1586-91.
[16.] De Vries M, Veerkamp R. Energy balance of dairy cattle in relation to milk production variables and fertility. J Dairy Sci 2000;83:62-9.
[17.] Lee S, Cho K, Park M, et al. Genetic parameters of milk beta- hydroxybutyric acid and acetone and their genetic association with milk production traits of Holstein cattle. Asian-Australas J Anim Sci 2016;29:1530-40.
[18.] Veerkamp R, Oldenbroek J, Van Der Gaast H, Van Der Werf J. Genetic correlation between days until start of luteal activity and milk yield, energy balance, and live weights. J Dairy Sci 2000;83:577-83.
[19.] Spurlock DM, Dekkers JC, Fernando R, Koltes DA, Wolc A. Genetic parameters for energy balance, feed efficiency, and related traits in Holstein cattle. J Dairy Sci 2012;95:5393-402.
Umanthi Ranaraja (1), KwangHyun Cho (2), MiNa Park (2), SiDong Kim (2), SeokHyun Lee (1), and ChangHee Do (1) *
* Corresponding Author: ChangHee Do Tel: +82-42-821-5784, Fax: +82-42-823-2766, E-mail: firstname.lastname@example.org
(1) Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea
(2) National Institute of Animal Science, RDA, Cheonan 31000, Korea
Submitted Jun 12, 2017; Revised Sept 16, 2017; Accepted Jan 10, 2018
Caption: Figure 1. The relationship between milk BHBA and Ac concentration with FPR along with lactation week (FPR, fat to protein ratio; Ac, acetone; BHBA, [beta]-hydroxybutyrate). The bar graph represents the FPR values and line graphs represent the ketone body concentrations.
Caption: Figure 2. The relationship between milk BHBA and acetone concentration with energy lactation week (EB, energy balance; Ac, acetone; BHBA, [beta]-hydroxybutyrate). The bar values and line graphs represent the ketone body concentrations.
Caption: Figure 3. Heritability estimates of milk BHBA, milk Ac, EB, and FPR according to parities from parity 1-4 (FPR, fat to protein ratio; EB, energy balance; Ac, Acetone, BHBA-[beta]-hydroxybutyrate).
Table 1. Descriptive statistics for the traits including milk [beta]-hydroxybutyrate (BHBA) and acetone Traits No of Mean SD Minimum Maximum records Protein (%) 262,940 3.27 0.35 0.4 15.63 Fat (%) 262,940 3.84 0.92 0.19 18.89 BHBA 262,940 43.48 59.48 0 4,280 ([micro]mol/L) Acetone 262,940 127.49 91.76 0 3,570 ([micro]mol/L) Milk yield 262,940 16.11 5.15 11.8 98 (kg/d) SD, standard deviation. Table 2. Heritabilities (in bold in the diagonal) and genetic correlations (above the diagonal) and phenotypic correlations (below diagonal) Item BHBA Ac EB FPR BHBA 0.19 0.52 -0.28 0.40 Ac 0.57 0.28 0.22 -0.15 EB -0.24 0.09 0.50 -0.86 FPR 0.35 -0.01 -0.81 0.46 BHBA, [beta]-hydroxybutyrate; Ac, acetone; EB, energy balance; FPR, fat protein ratio. Table 3. Additive genetic variance ([[sigma].sup.2.sub.A]), phenotypic variances ([[sigma].sup.2.sub.P]), and heritability ([h.sup.2]) estimates of different parities Item Parity 1 Parity 2 Parity 3 BHBA ([[sigma].sup.2.sub.A]) 178.68 145.38 140.58 ([[sigma].sup.2.sub.P]) 1,203 1,319.8 1,559.1 ([h.sup.2]) 0.14 0.11 0.09 Ac ([[sigma].sup.2.sub.A]) 1,085.4 1,320.7 1,650.9 ([[sigma].sup.2.sub.P]) 4,710.2 4,446.7 5,164.7 ([h.sup.2]) 0.23 0.29 0.31 EB ([[sigma].sup.2.sub.A]) 117.88 116.44 86.343 ([[sigma].sup.2.sub.P]) 224.82 277.15 291 ([h.sup.2]) 0.52 0.42 0.29 FPR ([[sigma].sup.2.sub.A]) 0.0338 0.0352 0.0203 ([[sigma].sup.2.sub.P]) 0.0788 0.0938 0.0958 ([h.sup.2]) 0.42 0.37 0.21 Item Parity 4 All BHBA ([[sigma].sup.2.sub.A]) 155.78 109.66 ([[sigma].sup.2.sub.P]) 1,626.8 832.14 ([h.sup.2]) 0.09 0.19 Ac ([[sigma].sup.2.sub.A]) 1,537.1 1,726.86 ([[sigma].sup.2.sub.P]) 5,226.3 5,780 ([h.sup.2]) 0.29 0.29 EB ([[sigma].sup.2.sub.A]) 54.839 115.88 ([[sigma].sup.2.sub.P]) 277.38 222.28 ([h.sup.2]) 0.19 0.52 FPR ([[sigma].sup.2.sub.A]) 0.0154 0.0438 ([[sigma].sup.2.sub.P]) 0.0956 0.0983 ([h.sup.2]) 0.16 0.44 BHBA, [beta]-hydroxybutyrate; Ac, acetone; EB, energy balance; FPR, fat protein ratio.
|Printer friendly Cite/link Email Feedback|
|Author:||Ranaraja, Umanthi; Cho, KwangHyun; Park, MiNa; Kim, SiDong; Lee, SeokHyun; Do, ChangHee|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Jun 1, 2018|
|Previous Article:||A demonstration of the H3 trimethylation ChIP-seq analysis of galline follicular mesenchymal cells and male germ cells.|
|Next Article:||Genetic diversity analysis of Thai indigenous chickens based on complete sequences of mitochondrial DNA D-loop region.|