Short communication: Principal components and factor analytic models for test-day milk yield in Brazilian Holstein cattle

A total of 46,089 individual monthly test-day (TD) milk yields (10 test-days), from 7,331 complete first lactations of Holstein cattle were analyzed. A standard multivariate analysis (MV), reduced rank analyses fitting the first 2, 3, and 4 genetic principal components (PC2, PC3, PC4), and analyses that fitted a factor analytic structure considering 2, 3, and 4 factors (FAS2, FAS3, FAS4), were carried out. The models included the random animal genetic effect and fixed effects of the contemporary groups (herd-year-month of test-day), age of cow (linear and quadratic effects), and days in milk (linear effect). The residual covariance matrix was assumed to have full rank. Moreover, 2 random regression models were applied. Variance components were estimated by restricted maximum likelihood method. The heritability estimates ranged from 0.11 to 0.24. The genetic correlation estimates between TD obtained with the PC2 model were higher than those obtained with the MV model, especially on adjacent test-days at the end of lactation close to unity. The results indicate that for the data considered in this study, only 2 principal components are required to summarize the bulk of genetic variation among the 10 traits.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

Validation of models for predicting milk urea nitrogen concentrations, estimating dry matter intake by the NRC (2001)

The objective of this study was to validate three different models for predicting milk urea nitrogen using field conditions, attempting to evaluate the nutritional adequacy diets for dairy cows and prediction of nitrogen excreted to the environment. Observations (4,749) from 855 cows were used. Milk yield, body weight (BW), days in milk and parity were recorded on the milk sampling days. Milk was sampled monthly, for analysis of milk urea nitrogen (MUN), fat, protein, lactose and total solids concentration and somatic cells count. Individual dry matter intake was estimated using the NRC (2001). The three models studied were derived from a first one to predict urinary nitrogen (UN). Model 1 was MUN = UN/12.54, model 2 was MUN = UN/17.6 and model 3 was MUN = UN/(0.0259 × BW), adjusted by body weight effect. To evaluate models, they were tested for accuracy, precision and robustness. Despite being more accurate (mean bias = 0.94 mg/dL), model 2 was less precise (residual error = 4.50 mg/dL) than model 3 (mean bias = 1.41 and residual error = 4.11 mg/dL), while model 1 was the least accurate (mean bias = 6.94 mg/dL) and the least precise (residual error = 5.40 mg/dL). They were not robust, because they were influenced by almost all the variables studied. The three models for predicting milk urea nitrogen were different with respect to accuracy, precision and robustness.

Effects of long-term feeding of genetically modified corn (event MON810) on the performance of lactating dairy cows

A long-term study over 25 months was conducted to evaluate the effects of genetically modified corn on performance of lactating dairy cows. Thirty-six dairy cows were assigned to two feeding groups and fed with diets based on whole-crop silage, kernels and whole-crop cobs from Bt-corn (Bt-MON810) or its isogenic not genetically modified counterpart (CON) as main components. The study included two consecutive lactations. There were no differences in the chemical composition and estimated net energy content of Bt-MON810 and CON corn components and diets. CON feed samples were negative for the presence of Cry1Ab protein, while in Bt-MON810 feed samples the Cry1Ab protein was detected. Cows fed Bt-MON810 corn had a daily Cry1Ab protein intake of 6.0 mg in the first lactation and 6.1 mg in the second lactation of the trial. Dry matter intake (DMI) was 18.8 and 20.7 kg/cow per day in the first and the second lactation of the trial, with no treatment differences. Similarly, milk yield (23.8 and 29.0 kg/cow per day in the first and the second lactation of the trial) was not affected by dietary treatment. There were no consistent effects of feeding MON810 or its isogenic CON on milk composition or body condition. Thus, the present long-term study demonstrated the compositional and nutritional equivalence of Bt-MON810 and its isogenic CON.

Identification of ABCA1 and ABCG1 in milk fat globules and mammary cells--implications for milk cholesterol secretion

The ATP-binding cassette (ABC) transporters ABCA1 and ABCG1 play an important role in cellular cholesterol homeostasis, but their function in mammary gland (MG) tissue remains elusive. A bovine MG model that allows repeated MG sampling in identical animals at different functional stages was used to test whether 1) ABCA1 and ABCG1 protein expression and subcellular localization in mammary epithelial cells (MEC) change during the pregnancy-lactation cycle, and 2) these 2 proteins were present in milk fat globules (MFG). Expression and localization in MEC were investigated in bovine MG tissues at the end of lactation, during the dry period (DP), and early lactation using immunohistochemical and immunofluorescence approaches. The presence of ABCA1 and ABCG1 in MFG isolated from fresh milk was determined by immunofluorescence. The ABCA1 protein expression in MEC, expressed as arbitrary units, was higher during the end of lactation (12.2±0.24) and the DP (12.5±0.22) as compared with during early lactation (10.2±0.65). In contrast, no significant change in ABCG1 expression existed between the stages. Throughout the cycle, ABCA1 and ABCG1 were detected in the apical (41.9±24.8 and 49.0±4.96% of cows, respectively), basal (56.2±28.1 and 54.6±7.78% of cows, respectively), or entire cytoplasm (56.8±13.4 and 61.6±14.4% of cows, respectively) of MEC, or showed combined localization. Unlike ABCG1, ABCA1 was absent at the apical aspect of MEC during early lactation. Immunolabeling experiments revealed the presence of ABCA1 and ABCG1 in MFG membranes. Findings suggest a differential, functional stage-dependent role of ABCA1 and ABCG1 in cholesterol homeostasis of the MG epithelium. The presence of ABCA1 and ABCG1 in MFG membranes suggests that these proteins are involved in cholesterol exchange between MEC and alveolar milk.

Comparison of the effect of a CIDR-Select synch versus a long-term CIDR based AI protocol on reproductive performance in multiparous dairy cows in Swiss dairy farms

Background Synchronization programs have become standard in the dairy industry in many countries. In Switzerland, these programs are not routinely used for groups of cows, but predominantly as a therapy for individual problem cows. The objective of this study was to compare the effect of a CIDR-Select Synch and a 12-d CIDR protocol on the pregnancy rate in healthy, multiparous dairy cows in Swiss dairy farms. Methods Cows (N = 508) were randomly assigned to CIDR-Select Synch (N = 262) or 12-d CIDR (N = 246) protocols. Cows in the CIDR-Select Synch group received a CIDR and 2.5 ml of buserelin i.m. on d 0. On d 7, the CIDR insert was removed and 5 ml of dinoprost was administered i.m.. Cows in the 12-d CIDR group received the CIDR on d 0 and it was removed on d 12 (the routine CIDR protocol in Swiss dairies). On d 0 a milk sample for progesterone analysis was taken. Cows were inseminated upon observed estrus. Pregnancy was determined at or more than 35 days after artificial insemination. As a first step, the two groups were compared as to indication for treatment, breed, stud book, stall, pasture, and farmer's business using chi square tests or Fisher's exact test. Furthermore, groups were compared as to age, DIM, number of AI's, number of cows per farm, and yearly milk yield per cow using nonparametric ANOVA. A multiple logistic model was used to relate the success of the protocols to all of the available factors; in particular treatment (CIDR-Select Synch/12-d CIDR), milk progesterone value, age, DIM, previous treatment of the uterus, previous gynecological treatment, and number of preceding inseminations. Results The pregnancy rate was higher in cows following the CIDR-Select Synch compared to the 12-d CIDR protocol (50.4% vs. 22.4%; P < 0.0001). Conclusion The CIDR-Select Synch protocol may be highly recommended for multiparous dairy cows. The reduced time span of the progesterone insert decreased the number of days open, improved the pregnancy rate compared to the 12-d CIDR protocol and the cows did not to have to be handled more often.

Effects of nutrient restriction on mammary cell turnover and mammary gland remodeling in lactating dairy cows

The aim of this study was to investigate the effects of a severe nutrient restriction on mammary tissue morphology and remodeling, mammary epithelial cell (MEC) turnover and activity, and hormonal status in lactating dairy cows. We used 16 Holstein x Normande crossbred dairy cows, divided into 2 groups submitted to different feeding levels (basal and restricted) from 2 wk before calving to wk 11 postpartum. Restricted-diet cows had lower 11-wk average daily milk yield from calving to slaughter than did basal-diet cows (20.5 vs. 33.5 kg/d). Feed restriction decreased milk fat, protein, and lactose yields. Restriction also led to lower plasma insulin-like growth factor 1 and higher growth hormone concentrations. Restricted-diet cows had lighter mammary glands than did basal-diet cows. The total amount of DNA in the mammary gland and the size of the mammary acini were smaller in the restricted-diet group. Feed restriction had no significant effect on MEC proliferation at the time of slaughter but led to a higher level of apoptosis in the mammary gland. Gelatin zymography highlighted remodeling of the mammary extracellular matrix in restricted-diet cows. Udders from restricted-diet cows showed lower transcript expression of alpha-lactalbumin and kappa-casein. In conclusion, nutrient restriction resulted in lower milk yield in lactating dairy cows, partly due to modulation of MEC activity and a lower number of mammary cells. An association was found between feed restriction-induced changes in the growth hormone-insulin-like growth factor-1 axis and mammary epithelial cell dynamics.

Performance and metabolic profile of dairy cows during a lactational and deliberately induced negative energy balance with subsequent realimentation

Homeorhetic and homeostatic controls in dairy cows are essential for adapting to alterations in physiological and environmental conditions. To study the different mechanisms during adaptation processes, effects of a deliberately induced negative energy balance (NEB) by feed restriction near 100 d in milk (DIM) on performance and metabolic measures were compared with lactation energy deficiency after parturition. Fifty multiparous cows were studied in 3 periods (1=early lactation up to 12 wk postpartum; 2=feed restriction for 3 wk beginning at 98+/-7 DIM with a feed-restricted and control group; and 3=a subsequent realimentation period for the feed-restricted group for 8 wk). In period 1, despite NEB in early lactation [-42 MJ of net energy for lactation (NE(L))/d, wk 1 to 3] up to wk 9, milk yield increased from 27.5+/-0.7 kg to a maximum of 39.5+/-0.8 kg (wk 6). For period 2, the NEB was induced by individual limitation of feed quantity and reduction of dietary energy density. Feed-restricted cows experienced a greater NEB (-63 MJ of NEL/d) than did cows in early lactation. Feed-restricted cows in period 2 showed only a small decline in milk yield of -3.1+/-1.1 kg and milk protein content of -0.2+/-0.1% compared with control cows (30.5+/-1.1 kg and 3.8+/-0.1%, respectively). In feed-restricted cows (period 2), plasma glucose was lower (-0.2+/-0.0 mmol/L) and nonesterified fatty acids higher (+0.1+/-0.1 mmol/L) compared with control cows. Compared with the NEB in period 1, the decreases in body weight due to the deliberately induced NEB (period 2) were greater (56+/-4 vs. 23+/-3 kg), but decreases in body condition score (0.16+/-0.03 vs. 0.34+/-0.04) and muscle diameter (2.0+/-0.4 vs. 3.5+/-0.4 mm) were lesser. The changes in metabolic measures in period 2 were marginal compared with the adjustments directly after parturition in period 1. Despite the greater induced energy deficiency at 100 DIM than the early lactation NEB, the metabolic load experienced by the dairy cows was not as high as that observed in early lactation. The different effects of energy deficiency at the 2 stages in lactation show that metabolic problems in early lactating dairy cows are not due only to the NEB, but mainly to the specific metabolic regulation during this period.

Endocrine changes and liver mRNA abundance of somatotropic axis and insulin system constituents during negative energy balance at different stages of lactation in dairy cows

The liver has an important role in metabolic regulation and control of the somatotropic axis to adapt successfully to physiological and environmental changes in dairy cows. The aim of this study was to investigate the adaptation to negative energy balance (NEB) at parturition and to a deliberately induced NEB by feed restriction at 100 days in milk. The hepatic gene expression and the endocrine system of the somatotropic axis and related parameters were compared between the early and late NEB period. Fifty multiparous cows were subjected to 3 periods (1=early lactation up to 12 wk postpartum, 2=feed restriction for 3 wk beginning at around 100 days in milk with a feed-restricted and a control group, and 3=subsequent realimentation period for the feed-restricted group for 8 wk). In period 1, plasma growth hormone reached a maximum in early lactation, whereas insulin-like growth factor-I (IGF-I), leptin, the thyroid hormones, insulin, and the revised quantitative insulin sensitivity check index increased gradually after a nadir in early lactation. Three days after parturition, hepatic mRNA abundance of growth hormone receptor 1A, IGF-I, IGF-I receptor and IGF-binding protein-3 (IGFBP-3) were decreased, whereas mRNA of IGFBP-1 and -2 and insulin receptor were upregulated as compared with wk 3 antepartum. During period 2, feed-restricted cows showed decreased plasma concentrations of IGF-I and leptin compared with those of control cows. The revised quantitative insulin sensitivity check index was lower for feed-restricted cows (period 2) than for control cows. Compared with the NEB in period 1, the changes due to the deliberately induced NEB (period 2) in hormones were less pronounced. At the end of the 3-wk feed restriction, the mRNA abundance of IGF-I, IGFBP-1, -2, -3, and insulin receptor was increased as compared with the control group. The different effects of energy deficiency at the 2 stages in lactation show that the endocrine regulation changes qualitatively and quantitatively during the course of lactation.

First report about the mode of action of combined butafosfan and cyanocobalamin on hepatic metabolism in nonketotic early lactating cows

The primary aim was to investigate the effect of combined butafosfan and cyanocobalamin on liver metabolism in early lactating cows through mRNA expression measurements of genes encoding 31 enzymes and transport proteins of major metabolic processes in the liver using 16 multiparous early lactating dairy cows. The treatments included i.v. injection of 10 mL/100 kg of body weight combined butafosfan and cyanocobalamin (TG, n = 8) on 3 d consecutively at 25 +/- 3 d in milk or injection with physiological saline solution similarly applied (CG, n = 8). Results include a higher daily milk production for TG cows (41.1 +/- 0.9 kg, mean +/- SEM) compared with CG cows (39.5 +/- 0.7 kg). In plasma, the concentration of inorganic phosphorus was lower in the TG cows (1.25 +/- 0.08 mmol/L) after the treatment than in the CG cows (1.33 +/- 0.07 mmol/L). The plasma beta-hydroxybutyrate concentration was 0.65 +/- 0.13 mmol/L for all cows before the treatment, and remained unaffected post treatment. The unique result was that in the liver, the mRNA abundance of acyl-coenzyme A synthetase long-chain family member 1, involved in fatty acid oxidation and biosynthesis, was lower across time points after the treatment for TG compared with CG cows (17.5 +/- 0.15 versus 18.1 +/- 0.24 cycle threshold, log(2), respectively). In conclusion, certain effects of combined butafosfan and cyanocobalamin were observed on mRNA abundance of a gene in the liver of nonketotic early lactating cows.

Induced hypoglycemia for 48 hours indicates differential glucose and insulin effects on liver metabolism in dairy cows

Hypoglycemia is a characteristic condition of early lactation dairy cows and is subsequently dependent on, and may affect, metabolism in the liver. The objective of the present study was to investigate the effects of induced hypoglycemia, maintained for 48 h, on metabolic parameters in plasma and liver of mid-lactation dairy cows. The experiment involved 3 treatments, including a hyperinsulinemic hypoglycemic clamp (HypoG, n=6) to obtain a glucose concentration of 2.5 mmol/L, a hyperinsulinemic euglycemic clamp (EuG, n=6) in which the effect of insulin was studied, and a control treatment with a 0.9% saline solution (NaCl, n=6). Blood samples for measurements of insulin, metabolites, and enzymes were taken at least once per hour. Milk yield was recorded and milk samples were collected before and after treatment. Liver biopsies were obtained before and after treatment to measure mRNA abundance by real-time, quantitative reverse transcription-PCR of 12 candidate genes involved in the main metabolic pathways. Milk yield decreased in HypoG and NaCl cows, whereas it remained unaffected in EuG cows. Energy-corrected milk yield (kg/d) was only decreased in HypoG cows. In plasma, concentration of beta-hydroxybutyrate decreased in response to treatment in EuG cows and was lower (0.41+/-0.04 mmol/L) on d 2 of the treatment compared with that in HypoG and NaCl cows (on average 0.61+/-0.03 mmol/L, respectively). Nonesterified fatty acids remained unaffected in all treatments. In the liver, differences between treatments for their effects were only observed in case of mitochondrial phosphoenolpyruvate carboxykinase (PEPCKm) and glucose-6-phosphatase (G6PC). In HypoG, mRNA abundance of PEPCKm was upregulated, whereas in EuG and NaCl cows, it was downregulated. The EuG treatment downregulated mRNA expression of G6PC, a marked effect compared with the unchanged transcript expression in NaCl. The mRNA abundance of the insulin receptor remained unaffected in all treatments, and no significant treatment differences were observed for genes related to lipid metabolism. In conclusion, low glucose concentrations in dairy cows affect liver metabolism at a molecular level through upregulation of PEPCKm mRNA abundance. Metabolic regulatory events in the liver are directed, apart from hormones, by the level of metabolites, either in excess (e.g., free fatty acids) or in shortage (e.g., glucose).

Effect of the prolactin-release inhibitor quinagolide on lactating dairy cows

In most mammals, prolactin (PRL) is essential for maintaining lactation, and yet the short-term suppression of PRL during established lactation by bromocriptine has produced inconsistent effects on milk yield in cows and goats. To assess the effect of the long-term inhibition of PRL release in lactating dairy cows, 5 Holstein cows in early lactation received daily intramuscular injections of 1mg of the PRL-release inhibitor quinagolide for 9 wk. Four control cows received the vehicle (water) only. During the last week of the treatments, one udder half was milked once a day (1x) and the other twice a day (2x). Blood samples were harvested at milking in wk -1, 1, 4, and 8. The daily injections of quinagolide reduced milking-induced PRL release but not the basal PRL concentration. Quinagolide induced a faster decline in milk production, which was about 5.3 kg/d lower in the quinagolide-treated cows during the last 4 wk of treatment. During wk 9, the inhibition of milk production by quinagolide was maintained in the udder half that was milked 2x but not in the half milked 1x. Milk production was significantly correlated with the quantity of PRL released at milking. Quinagolide did not affect the release of oxytocin at milking. Serum concentration of insulin-like growth factor-1 was not affected by treatment or correlated with milk production. Serum concentrations of leptin and the calciotropic hormone stanniocalcin were not affected by the treatment. In conclusion, the chronic administration of the PRL-release inhibitor quinagolide decreases milk production in dairy cows. The effect is likely the result of the reduced release of milking-induced PRL and is modulated at the level of the gland by milking frequency.

Coordinated gene expression in adipose tissue and liver differs between cows with high or low NEFA concentrations in early lactation

Dairy cows with high and low plasma non-esterified fatty acid (NEFA) concentrations in early lactation were compared for plasma parameters and mRNA expression of genes in liver and subcutaneous adipose tissue. The study involved 16 multiparous dairy cows with a plasma NEFA concentration of >500 mumol/l [n = 8, high NEFA (HNEFA)] and <140 mumol/l [n = 8, low NEFA (LNEFA)] in the first week post-partum (pp). Blood samples, adipose and liver tissues were collected on day 1 (+1d) and at week 3 pp (+3wk). Blood plasma was assayed for concentrations of metabolites and hormones. Subcutaneous adipose and liver tissues were analysed for mRNA abundance by real-time qRT-PCR encoding parameters related to lipid metabolism. Results showed that mean daily milk yield and milk fat quantity were higher in HNEFA than in LNEFA cows (p < 0.01), and the NEB was more negative in HNEFA than in LNEFA in +3wk too (p < 0.05). HNEFA cows had slightly lower (p < 0.1) insulin concentrations than LNEFA cows across the study period, and the body condition score decreased more from +1d to +3wk in HNEFA than in LNEFA (p = 0.09). The mRNA abundance of genes in the liver related to fatty acid oxidation (carnitine palmitoyltransferase 2 and very long chain acyl-coenzyme A dehydrogenase) and ketogenesis (3-hydroxy-3-methylglutaryl-coenzyme A synthase 2) were lower in HNEFA than in LNEFA cows. No differences between the two groups were observed for mRNA expression of genes in adipose tissue. The number of calculated significant correlation coefficients (moderately strong) between parameters in the liver and in adipose tissue was nearly similar on +1d, and higher for HNEFA compared with LNEFA cows in +3wk. In conclusion, dairy cows with high compared with low plasma NEFA concentrations in early lactation show differentially synchronized mRNA expression of genes in adipose tissue and liver in +3wk that suggests a different orchestrated homeorhetic regulation of lipid metabolism.

Milk flow-dependent vacuum loss in high-line milking systems: effects on milking characteristics and teat tissue condition

To study the effects of a milking system that partially compensates for milk flow-dependent vacuum loss compared with a standard (high-line) milking unit in a tie-stall barn, milk flow and vacuum patterns were recorded in 10 cows during machine milking with 2 milking systems in a crossover design for 7 d each. Before and after each treatment period postmilking teat condition was recorded by ultrasound cross-sectioning. Additionally, 2 methods to measure teat tissue condition were compared: longitudinal teat ultrasound cross-sectioning and teat tissue density measurements with the spring-loaded caliper (cutimeter method). The partial compensation of milk flow-dependent vacuum loss caused an elevation of the peak flow rate (4.74+/-0.08 vs. 4.29+/-0.07 kg/min) and a shorter duration of plateau (1.57+/-0.06 vs. 1.96+/-0.07 min) compared with the standard milking system. Total milk yield, duration of incline and decline of milk flow, average milk flow, time until peak flow rate, main milking time, and total milking time did not differ between treatments (overall means: 13.75+/-0.17 kg; 0.65+/-0.01 min; 2.88+/-0.09 min; 2.82+/-0.05 kg/min; 1.65+/-0.03 min; 5.23+/-0.09 min, and 5.30+/-0.10 min, respectively). The vacuum drop in the short milk tube during periods of high milk flow was less in the compensating vacuum than in the standard milking system (11+/-1.1 vs. 15+/-0.7 kPa). Teat measures as determined by ultrasound remained unchanged over the entire experimental period with both milking systems. Postmilking teat tissue measures including their recovery within 20 min after the end of milking show a correlation (0.85 and 0.71, respectively) between the methods used (ultrasound and cutimeter method). In conclusion, a more constant vacuum at the teat tip (within the short milk tube) during periods of high milk flow affected milk flow patterns, mainly increasing peak flow rate. However, the reduced vacuum loss did not increase the overall speed of milking. In addition, effects of higher vacuum stability on teat condition and udder health were not obvious.

A field study on characteristics and diversity of gene expression in the liver of dairy cows during the transition period

Metabolic and endocrine adaptations to support milk production during the transition period vary between individual cows. This variation between cows to adapt to lactation may have a genetic basis. The present field study was carried out to determine hepatic adaptations occurring from late pregnancy through early lactation by measuring mRNA abundance of candidate genes in dairy cows on-farm. Additionally, the objective was to observe the diversity in inter-individual variation for the candidate genes that may give indications where individual adaptations at a molecular level can be found. This study was carried out on-farm including 232 dairy cows (parity >3) from 64 farms in Switzerland. Blood and liver samples were collected on d 20+/-7 before parturition, on d 24+/-2, and on d 89+/-4 after parturition. Blood plasma was assayed for concentrations of glucose, nonesterified fatty acids, beta-hydroxybutyrate, cholesterol, triglycerides, urea, albumin, protein, insulin, insulin-like growth factor-1, leptin, 3,5,3'-triiodothyronine, and thyroxine. Liver samples were obtained at the same time points and were measured for mRNA abundance of 26 candidate genes encoding enzymes and nuclear receptors involved in gluconeogenesis, fatty acid beta-oxidation, fatty acid and triglyceride synthesis, ketogenesis, citric acid cycle, cholesterol synthesis, and the urea cycle. The cows in the present study experienced a marked metabolic load in early lactation, as presented by changes in plasma metabolites and hormones, and responded accordingly with upregulation and downregulation of almost all candidate genes involved in metabolic processes in the liver. The observed inter-individual variation for the candidate genes, which was highest for acetyl-CoA-carboxylase and glycerol-3-phosphate dehydrogenase 2, should be further investigated to unravel the regulation at molecular level for optimal adaptive performance in dairy cows.