Body and milk traits as indicators of dairy cow energy status in early lactation

The inclusion of feed intake and efficiency traits in dairy cow breeding goals can lead to increased risk of metabolic stress. An easy and inexpensive way to monitor postpartum energy status (ES) of cows is therefore needed. Cows' ES can be estimated by calculating the energy balance from energy intake and output and predicted by indicator traits such as change in body weight (Delta BW), change in body condition score (Delta BCS), milk fat:protein ratio (FPR), or milk fatty acid (FA) composition. In this study, we used blood plasma nonesterified fatty acids (NEFA) concentration as a biomarker for ES. We determined associations between NEFA concentration and ES indicators and evaluated the usefulness of body and milk traits alone, or together, in predicting ES of the cow. Data were collected from 2 research herds during 2013 to 2016 and included 137 Nordic Red dairy cows, all of which had a first lactation and 59 of which also had a second lactation. The data included daily body weight, milk yield, and feed intake and monthly BCS. Plasma samples for NEFA were collected twice in lactation wk 2 and 3 and once in wk 20. Milk samples for analysis of fat, protein, lactose, and FA concentrations were taken on the blood sampling days. Plasma NEFA concentration was higher in lactation wk 2 and 3 than in wk 20 (0.56 +/- 0.30, 0.43 +/- 0.22, and 0.13 +/- 0.06 mmol/L, respectively; all means +/- standard deviation). Among individual indicators, C18:1 cis-9 and the sum of C18:1 in milk had the highest correlations (r = 0.73) with NEFA. Seven multiple linear regression models for NEFA prediction were developed using stepwise selection. Of the models that included milk traits (other than milk FA) as well as body traits, the best fit was achieved by a model with milk yield, FPR, Delta BW, Delta BCS, FPR x Delta BW, and days in milk. The model resulted in a cross-validation coefficient of determination (R(2)cv) of 0.51 and a root mean squared error (RMSE) of 0.196 mmol/L. When only milk FA concentrations were considered in the model, NEFA prediction was more accurate using measurements from evening milk than from morning milk (R(2)cv = 0.61 vs. 0.53). The best model with milk traits contained FPR, C10:0, C14:0, C18:1 cis-9, C18:1 cis-9 x C14:0, and days in milk (R(2)cv = 0.62; RMSE = 0.177 mmol/L). The most advanced model using both milk and body traits gave a slightly better fit than the model with only milk traits (R(2)cv = 0.63; RMSE = 0.176 mmol/L). Our findings indicate that ES of cows in early lactation can be monitored with moderately high accuracy by routine milk measurements.Peer reviewe

The effects of dietary resin acid inclusion on productive, physiological and rumen microbiome responses of dairy cows during early lactation

Dairy cows have intense fluctuations in digestive, metabolic and hormonal systems around calving which predispose them to various disorders and health problems. The aim of the current experiment was to investigate feed and nutrient intake, rumen fermentation, rumen bacterial communities, milk production, milk fatty acid composition and plasma biomarker profiles of dairy cows to assess the modulation of these functions by in-feed resin acid inclusion. Thirty-six Nordic Red cows were used in a continuous feeding trial starting 3 weeks prepartum and lasting for 10 weeks into the lactation. The cows were fed grass silage ad libitum and the dietary treatments were 1) control with basal concentrate (CON), 2) CON supplemented with tall oil fatty acids (TOFA; 90 % fatty acids and 9% resin acids) at 7.0 g/cow/day and 3) CON supplemented with resin acid concentrate (RAC; 37.5% resin acids) at 1.7 g/cow/day. The mixture of resin acids in TOFA and RAC, consisting mostly of abietic and dehydroabietic acids, originated from coniferous tree species Pinus sylvestris L. and Picea abies L. Feed intake and milk production were measured throughout the experimental period. Milk and blood samples were collected at weeks 2, 3, 6 and 10, and rumen fluid was sampled at weeks 2 and 10 of lactation to analyse rumen fermentation and rumen bacterial communities. The dynamics in feed intake and milk production with progressing lactation showed typical curvilinear trends (P for timePeer reviewe

Simulation Study on Heterogeneous Variance Adjustment for Observations with Different Measurement Error Variance

Abstract The Nordic Holstein yield evaluation model describes all available milk, protein and fat test-day yields from Denmark, Finland and Sweden. In its current form all variance components are estimated from observations recorded under conventional milking systems. Also the model for heterogeneity of variance correction is developed for the same observations. As automated milking systems are becoming more popular the current evaluation model needs to be enhanced to account for the different measurement error variances of observations from automated milking systems. In this simulation study different models and different approaches to account for heterogeneous variance when observations have different measurement error variances were investigated. Based on the results we propose to upgrade the currently applied models and to calibrate the heterogeneous variance adjustment method to yield same genetic variance for both milking systems

Principal component approach in variance component estimation for international sire evaluation

<p>Abstract</p> <p>Background</p> <p>The dairy cattle breeding industry is a highly globalized business, which needs internationally comparable and reliable breeding values of sires. The international Bull Evaluation Service, Interbull, was established in 1983 to respond to this need. Currently, Interbull performs multiple-trait across country evaluations (MACE) for several traits and breeds in dairy cattle and provides international breeding values to its member countries. Estimating parameters for MACE is challenging since the structure of datasets and conventional use of multiple-trait models easily result in over-parameterized genetic covariance matrices. The number of parameters to be estimated can be reduced by taking into account only the leading principal components of the traits considered. For MACE, this is readily implemented in a random regression model.</p> <p>Methods</p> <p>This article compares two principal component approaches to estimate variance components for MACE using real datasets. The methods tested were a REML approach that directly estimates the genetic principal components (direct PC) and the so-called bottom-up REML approach (bottom-up PC), in which traits are sequentially added to the analysis and the statistically significant genetic principal components are retained. Furthermore, this article evaluates the utility of the bottom-up PC approach to determine the appropriate rank of the (co)variance matrix.</p> <p>Results</p> <p>Our study demonstrates the usefulness of both approaches and shows that they can be applied to large multi-country models considering all concerned countries simultaneously. These strategies can thus replace the current practice of estimating the covariance components required through a series of analyses involving selected subsets of traits. Our results support the importance of using the appropriate rank in the genetic (co)variance matrix. Using too low a rank resulted in biased parameter estimates, whereas too high a rank did not result in bias, but increased standard errors of the estimates and notably the computing time.</p> <p>Conclusions</p> <p>In terms of estimation's accuracy, both principal component approaches performed equally well and permitted the use of more parsimonious models through random regression MACE. The advantage of the bottom-up PC approach is that it does not need any previous knowledge on the rank. However, with a predetermined rank, the direct PC approach needs less computing time than the bottom-up PC.</p

Principal component and factor analytic models in international sire evaluation

<p>Abstract</p> <p>Background</p> <p>Interbull is a non-profit organization that provides internationally comparable breeding values for globalized dairy cattle breeding programmes. Due to different trait definitions and models for genetic evaluation between countries, each biological trait is treated as a different trait in each of the participating countries. This yields a genetic covariance matrix of dimension equal to the number of countries which typically involves high genetic correlations between countries. This gives rise to several problems such as over-parameterized models and increased sampling variances, if genetic (co)variance matrices are considered to be unstructured.</p> <p>Methods</p> <p>Principal component (PC) and factor analytic (FA) models allow highly parsimonious representations of the (co)variance matrix compared to the standard multi-trait model and have, therefore, attracted considerable interest for their potential to ease the burden of the estimation process for multiple-trait across country evaluation (MACE). This study evaluated the utility of PC and FA models to estimate variance components and to predict breeding values for MACE for protein yield. This was tested using a dataset comprising Holstein bull evaluations obtained in 2007 from 25 countries.</p> <p>Results</p> <p>In total, 19 principal components or nine factors were needed to explain the genetic variation in the test dataset. Estimates of the genetic parameters under the optimal fit were almost identical for the two approaches. Furthermore, the results were in a good agreement with those obtained from the full rank model and with those provided by Interbull. The estimation time was shortest for models fitting the optimal number of parameters and prolonged when under- or over-parameterized models were applied. Correlations between estimated breeding values (EBV) from the PC19 and PC25 were unity. With few exceptions, correlations between EBV obtained using FA and PC approaches under the optimal fit were ≥ 0.99. For both approaches, EBV correlations decreased when the optimal model and models fitting too few parameters were compared.</p> <p>Conclusions</p> <p>Genetic parameters from the PC and FA approaches were very similar when the optimal number of principal components or factors was fitted. Over-fitting increased estimation time and standard errors of the estimates but did not affect the estimates of genetic correlations or the predictions of breeding values, whereas fitting too few parameters affected bull rankings in different countries.</p

Estimation of social genetic effects on feeding behaviour and production traits in pigs

Pigs are housed in groups during the test period. Social effects between penmates may affect average daily gain (ADG), backfat thickness (BF), feed conversion rate (FCR), and the feeding behaviour traits of pigs sharing the same pen. The aim of our study was to estimate the genetic parameters of feeding behaviour and production traits with statisticalmodels that include social genetic effects (SGEs). The data contained 3075 Finnish Yorkshire, 3351 Finnish Landrace, and 968 F1-crossbred pigs. Feeding behaviour traits were measured as the number of visits per day (NVD), time spent in feeding per day (TPD), daily feed intake (DFI), time spent in feeding per visit (TPV), feed intake per visit (FPV), and feed intake rate (FR). The test period was divided into five periods of 20 days. The number of pigs per pen varied from 8 to 12. Two model approaches were tested, i.e. a fixed group size model and a variable group size model. For the fixed group size model, eight random pigs per pen were included in the analysis, while all pigs in a pen were included for the variable group size model. The linear mixed-effectsmodel included sex, breed, and herd*year*season as fixed effects and group (batch*pen), litter, the animal itself (direct genetic effect (DGE)), and penmates (SGEs) as random effects. For feeding behaviour traits, estimates of the total heritable variation (T-2 +/- SE) and classical heritability (h(2) +/- SE, values given in brackets) from the variable group size model (e.g. period 1) were 0.34 +/- 0.13 (0.30 +/- 0.04) for NVD, 0.41 +/- 0.10 (0.37 +/- 0.04) for TPD, 0.40 +/- 0.15 (0.14 +/- 0.03) for DFI, 0.53 +/- 0.15 (0.28 +/- 0.04) for TPV, 0.66 +/- 0.17 (0.28 +/- 0.04) for FPV, and 0.29 +/- 0.13 (0.22 +/- 0.03) for FR. The effect of social interaction was minimal for ADG (T-2 = 0.29 +/- 0.11 and h(2) = 0.29 +/- 0.04), BF (T-2 = 0.48 +/- 0.12 and h(2) = 0.38 +/- 0.07), and FCR (T-2 = 0.37 +/- 0.12 and h(2) = 0.29 +/- 0.04) using the variable group size model. In conclusion, the results indicate that social interactions have a considerable indirect genetic effect on the feeding behaviour and FCR of pigs but not on ADG and BF. (C) 2020 The Authors. Published by Elsevier Inc. on behalf of The Animal Consortium.Peer reviewe

Estimation of parameters for heterogeneous variance adjustment on test-day data

vo

Genetic parameters for residual energy intake and energy conversion efficiency in Nordic Red dairy cattle

201

Short communication: Multivariate outlier detection for routine Nordic dairy cattle genetic evaluation in the Nordic Holstein and Red population

201